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Chapter  1 


INTRODUCTION 


It  is  well  known  that  organic  matrix  fiber-reinforced  laminated  composites  are  very 
susceptible  to  transverse  low-velocity  impact,  which  could  cause  significant  damage 
embedded  within  the  materials.  Such  damage  involves  extensive  intraply  matrix- 
cracking  and  interply  delaminations.  Depending  on  the  extent  of  the  damage,  the 
strength  and  stiffness  of  the  materials  can  be  significantly  reduced  [1-12].  Therefore, 
the  knowledge  of  impact  damage  in  laminated  composites  is  critically  important  for 
the  application  of  the  materials  in  structural  design. 

Hence,  this  subject  has  attracted  significant  attention  among  researchers  recently. 
Considerable  work  both  experimental  and  analytical  has  been  performed  [13-52].  and 
most  of  the  studies  are  focused  on  the  damage  resulting  from  a  point-nose  projec¬ 
tile.  With  point-nose  impact  the  damage  pattern  is  very  complicated  and  three- 
dimensional,  involving  multiple  delaminations  along  interfaces  and  a  considerable 
number  of  intraply  matrix  cracks.  Because  the  damage  is  mostly  embedded  inside 
the  materials,  it  is  very  difficult  to  detects.  Although  several  new  techniques  have 
been  cited  in  the  recent  literature  [53],  C-scan  and  X-ray  machines  are  still  most 
commonly  used  for  evaluating  impact  damage  in  composites.  Specimens  are  also  fre¬ 
quently  sliced  into  pieces  for  visual  inspection  of  internal  damage.  These  procedures 
are  still  very  time-consuming  and  tedious. 

Some  of  the  prominent  work  in  this  area  is  briefly  mentioned  in  the  following: 
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Sjoblom  et  al  [25]  used  a  pendulum  type  impactor  with  strain  gauges  attached  to 
obtain  load  histories.  They  found  that  permanent  damage  was  indicated  by  jumps 
in  measured  energy  loss.  Clark  et  al  [33]  proposed  a  model  which  explained  some 
characteristics  of  impact  damage,  such  as  the  relative  sizes  of  delaminations  and 
occurrence  of  prominent  delamination/cracking  features.  Gosse  et  al  [30]  proposed  a 
K-rule  which  hypothesized  that  delaminations  would  occur  at  the  interface  bounded 
by  the  matrix  cracks  between  two  adjacent  plies  with  different  orientations.  Joshi 
et  al  [31,46]  studied  experimentally  impact  damage  in  laminated  composite  plates 
for  damage  mechanism  and  characterization.  They  showed  pictures  of  a  cross  ply 
plate  sliced  longitudinally  and  transversely  for  inspection  after  impact.  Quantitative 
discussion  about  failure  mechanism  was  presented.  Stori  et  al  [40]  performed  extensive 
impact  tests  to  evaluate  the  impact  resistance  of  carbon  fiber  reinforced  composites 
with  various  matrix  materials. 

Although  some  progress  has  been  made  in  understanding  impact  damage,  the 
knowledge  of  impact  damage  is  still  limited  by  the  complexity  of  the  impact  damage 
modes  resulting  from  the  previous  experiments.  Accordingly,  information  of  the  im¬ 
pact  damage  mechanisms  and  mechanics  is  still  not  well  developed,  and  the  governing 
parameters  controlling  the  impact  damage  in  laminated  composites  are  also  not  fully 
understood. 

In  addition  to  experimental  studies,  considerable  analytical  work  has  also  been 
performed  by  several  investigators.  Sun  et  al  [20,23]  developed  an  empirical  contact 
law  for  simulating  the  contact  force  distribution  of  a  projectile  on  impacted  com¬ 
posites.  A  plate  finite  element  method  associated  with  the  contact  law  was  used  to 
evaluate  transient  dynamic  behavior  of  an  impacted  composite  plate.  Wu  et  al  [13-16] 
performed  an  experimental  and  analytical  study  on  the  transverse  impact  behavior  of 
Graphite/Epoxy  laminated  composite  plate;  a  non-dimensional  empirical  expression 
was  developed  and  used  to  predict  the  damage  induced  by  low  velocity  impact.  Gu  et 
al  [361  developed  a  model  for  estimating  the  impact  damage  size  in  SMC  composites. 
However,  due  to  the  lack  of  understanding  of  the  basic  mechanisms  and  mechanics  of 
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impact  damage,  so  far  there  are  no  analytical  models  available  for  accurately  predict¬ 
ing  impact  damage  initiation  and  the  extent  of  the  damage  in  laminated  composites. 

Therefore,  in  this  investigation,  the  major  focus  of  the  attention  is  first  con¬ 
centrated  on  understanding,  through  experiments  as  well  as  analysis,  the  basic  im¬ 
pact  damage  mechanisms  and  mechanics  of  laminated  composites  resulting  from  low- 
velocity  impact.  The  basic  parameters  governing  the  impact  damage  mechanics  and 
mechanisms  are  then  identified  and  thoroughly  studied.  The  knowledge  gained  from 
the  studies  is  then  utilized  to  develop  an  analytical  model  for  predicting  damage  in 
composites  as  a  result  of  low-velocity  impact. 

In  Chapter  2,  the  major  approach  for  achieving  the  objectives  of  the  investigation 
is  briefly  outlined.  Chapter  3  describes  the  development  of  a  new  impact  facility 
and  a  new  testing  approach  for  studying  impact  damage  in  laminated  composites. 
Simplified  impact  damage  patterns  will  be  presented  and  the  major  characteristics 
of  impact  damage  in  laminated  composites  will  be  discussed.  Chapter  4  presents  an 
analytical  model  which  consists  of  a  stress  analysis  and  a  failure  analysis  for  analyzing 
the  experimental  results  obtained  from  the  tests  in  the  previous  chapter.  In  Chapter  5. 
effects  of  impact  parameters  such  as  stacking  sequence,  thickness  of  the  laminate  and 
mass  of  impactor  on  the  impact  damage  are  studied.  Both  experimental  and  analytical 
results  are  presented.  In  Chapter  6,  a  model  will  be  presented  for  predicting  impact 
damage  initiation  and  the  extent  of  the  damage  in  laminated  composites  resulting 
from  a  point  nose  impactor.  The  comparisons  between  the  numerical  simulations  and 
the  point  nose  impact  test  data  obtained  from  the  investigation  will  also  be  shown. 
The  information  regarding  the  3-D  impact  computer  code  will  be  described  in  Chapter 
7. 
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Chapter  2 


METHOD  OF  APPROACH 


2.1  THE  OBJECTIVES 

This  investigation  has  been  performed  to  study  the  impact  damage  in  fiber-reinforced 
polymer  matrix  laminated  composites  subjected  to  low-velocity  impact  by  a  foreign 
point-nose  projectile  as  shown  in  Figure  2.1. 

The  major  objectives  of  the  investigation  are  twofold: 

1.  to  study  impact  damage  mechanisms  and  mechanics  of  fiber- reinforced  lami¬ 
nated  composites  due  to  low-velocity  impact,  and 

2.  to  develop  adequate  models  for  predicting  the  low-velocity  impact  damage  in 
the  materials. 

In  order  to  achieve  the  objectives,  this  investigation  was  carried  out  in  three 
sequential  steps  as  follows: 

1.  Generation  of  Simplified  2-Dimensional  Impact  Test  Data. 

2.  Analysis  of  Simplified  2-Dimensional  Impact  Test  Data. 

3.  Development  of  Point-nose  Impact  Damage  Model. 
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Figure  2.1:  Description  of  the  problem.  A  laminated  composite  panel  subjected  to 
transverse  impact  by  a  low-velocity  spherical  nose  projectile. 


The  major  tasks  executed  in  each  step  are  briefly  outlined  in  the  next  section. 

2.2  THE  MAJOR  TASKS 

2.2.1  Generation  of  Simplified  2- Dimensional  Impact  Test  Data. 

A  unique  impact  test  facility  was  designed  and  built  for  this  study.  In  this  investi¬ 
gation,  the  line- nose  impactors  were  selected  because  they  could  produce  uniformly 
distributed  transient  dynamic  load  across  the  specimen.  As  a  consequence,  a  uniform 
and  consistent  impact  damage  pattern  could  be  produced  in  composites,  significantly 
simplifying  the  impact  damage  pattern.  Accordingly,  extensive  tests  were  performed 
to  study  the  impact  damage  mechanisms  and  failure  modes  of  laminated  composites 
subjected  to  line-loading  impact.  The  initial  failure  mode  and  failure  mechanisms 
were  thoroughly  examined  from  the  test  data.  The  effects  of  ply  orientation,  iliick- 
ness,  and  mass  of  the  impactor  on  impact  damage  were  also  evaluated. 
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2.2.2  Analysis  of  Simplified  2-Dimensional  Impact  Test  Data. 


An  analytical  model  was  also  developed  for  predicting  damage  in  laminated  compos¬ 
ites  resulting  from  the  line-loading  impact.  Dynamic  stresses  and  strains  in  compos¬ 
ites  were  calculated  from  a  transient  dynamic  finite  element  code  based  on  a  plane 
strain  condition,  which  was  developed  during  the  work.  A  modified  Hashin  matrix 
failure  criterion  [54]  was  proposed  to  predict  the  initial  damage,  and  a  post-failure 
analysis  after  the  occurrence  of  the  initial  damage  was  also  performed  to  determine 
the  effect  of  the  initial  damage  on  the  initiation  of  delamination  and  micro-matrix 
cracks. 


2.2.3  Development  of  Point-nose  Impact  Damage  Model. 

In  this  work,  an  analytical  model  was  developed  based  on  the  previous  line-loading 
impact  results,  by  extending  the  previous  two-dimensional  model  to  predict  the  im¬ 
pact  damage  resulting  from  a  point-nose  impactor.  The  model  consists  of  a  stress 
analysis  for  determining  the  stress  distributions  inside  the  materials  during  impact, 
and  a  failure  analysis  for  predicting  the  initiation  and  the  extent  of  the  impact  dam¬ 
age.  The  information  on  the  transient  dynamic  stresses  is  required  for  accurately 
predicting  impact  damage  from  the  failure  analysis.  In  the  failure  analysis,  a  matrix 
failure  criterion  and  a  delamination  criterion  were  proposed  for  predicting  the  initial 
impact  damage  and  the  extent  of  the  damage  due  to  impact,  respectively.  Point-nose 
impact  experiments  were  also  performed  during  the  investigation. 
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Chapter  3 


2-DIMENSIONAL  LINE-LOADING  IMPACT  EXPERIMENT 


3.1  INTRODUCTION 

The  objective  of  this  experiment  is  to  gain  a  fundamental  understanding  of  damage 
mechanisms  and  mechanics  of  laminated  composites  due  to  low-velocity  impact  and 
to  determine  essential  parameters  governing  impact  damage.  It  was  observed  from 
the  previous  experiments  that  impact  damage  can  be  significantly  affected  by  the 
laminate  configuration,  including  ply  orientation  and  thickness,  and  the  energy  of  the 
impactor.  Similar  phenomena  have  also  been  recorded  in  the  literature  for  impact 
damage  resulting  from  traditional  point-nose  impactors  [14,30,31,44,46,47].  Appar¬ 
ently,  there  exist  relationships  that  relate  the  impact  damage  to  the  material  proper¬ 
ties,  geometry,  ply  orientation,  and  mass  and  velocity  of  the  impactor.  Because  the 
experimental  procedures  for  inspecting  the  embedded  damage  caused  by  a  point-nose 
impactor  are  very  tedious  and  time-consuming,  it  is  very  difficult  to  rely  simply  on  the 
point-nose  impact  test  to  establish  these  relationships  and  fundamentally  understand 
these  phenomena.  None  of  these  relationships  have  been  established  experimentally 
or  analytically  in  the  literature.  However,  it  is  believed  that  the  knowledge  of  such 
information  is  critically  important  to  design  engineers  in  selecting  proper  materials, 
layups  and  geometry  for  designing  composite  structures  with  better  impact  resistance. 

In  order  to  achieve  this  objective,  a  new  impact  tester  was  designed  and  built  for 
the  investigation.  The  impactors  were  designed  with  a  line-edge  nose,  rather  than 
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a  point  nose,  which  considerably  reduced  the  complexity  of  the  impact  damage  into 
several  basic  patterns.  The  major  advantages  of  using  the  line-loading  impactor  are 
that  not  only  can  the  damage  pattern  be  simplified,  but  most  importantly,  it  can  be 
inspected  easily  from  the  sides  of  the  specimens,  without  resorting  to  sophisticated 
nondestructive  testing  machines  such  as  a  C-scan  or  an  X-ray. 


3.2  A  NEW  IMPACT  TESTER 

The  major  apparatus  of  the  facility  consists  of  a  pressure  tank,  a  precision-made 
barrel,  a  high  precision  timer,  optical  fiber  photoelectric  sensors,  and  supporting 
fixtures,  as  shown  in  Figure  3.1.  The  essential  characteristics  of  the  design,  which 
is  different  from  any  other  available  in  the  literature,  are  the  use  of  a  barrel  with  a 
rectangular  cross  section  and  the  impactors.  The  barrel  is  made  of  four  carbon-steel 
strips,  precisely  assembled  with  the  maximum  tolerance  of  the  dimensions  of  the  inner 
cross  section  varied  along  the  axis  of  the  barrel  within  0.01  cm  The  inner  dimensions 
of  the  barrel  are  7.62  cm  in  height  and  6.35  cm  in  width,  with  a  wall  thickness  of  1.27 
cm  The  length  of  the  barrel  is  about  137  cm 

Because  the  barrel  is  rectangular,  the  impactors  can  be  designed  in  two  parts, 
a  common  rectangular  base  and  a  nose  which  can  be  made  into  different  shapes, 
which  can  be  easily  mounted  and  dismounted  from  the  base  (see  Figure  3.2).  Thus, 
different  damage  patterns  can  be  produced  by  selecting  the  appropriate  nose  shape 
for  the  impactor.  For  instance,  the  impact  damage  caused  by  a  spherical  nose  is 
considerably  different  from  that  produced  by  a  line-nose  impactor. 

During  the  test,  the  impactor  was  driven  by  compressed  air  from  the  air  tank 
through  the  rectangular  barrel.  The  velocity  of  the  impactor  was  controlled  by  se¬ 
lecting  proper  weights  for  the  base  and  the  nose  of  the  impactor,  and  by  adjusting  the 
air  pressure  from  the  air  tank.  The  setup  of  the  tester  is  designed  so  that  once  the 
impactor  strikes  the  target,  it  immediately  rebounds  back  into  the  barrel.  Therefore, 
this  facility  can  be  used  to  evaluate  the  impact  damage  in  composites  as  a  function 
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Figure  3.1:  A  schematic  of  the  impact  test  facility. 


Line  Loading  Impactor 


Figure  3.2:  Description  of  the  problems  and  impactors.  (j4):  a  line-nose  impactor. 
( B ):  a  point-nose  impactor. 
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of  the  weight  and  velocity  of  the  impactor  and  the  nose  shape  of  the  impactor.  For 
instance,  the  velocity  of  a  60  gram  impactor  can  be  achieved  as  high  as  40  m/sec. 
The  heavier  the  mass  of  an  impactor,  the  lower  the  maximum  velocity  will  be. 

For  this  investigation,  a  line-nose  impactor  was  chosen,  as  shown  in  Figure  3.2. 
The  use  of  the  line-nose  impactor  results  in  an  uniformly  distributed  transient  dy¬ 
namic  load  across  the  specimen  which  is  clamped  on  two  parallel  free  edges.  It  was 
expected  that  such  a  uniform  loading  would  produce  a  consistent  and  uniform  damage 
pattern  throughout  the  specimen  width,  hence  substantially  simplifying  the  impact 
damage  mechanisms  from  a  three-dimensional  to  a  two-dimensional  event. 

3.3  EXPERIMENTS 

Extensive  impact  tests  were  performed  to  study  impact  damage  in  laminated  compos¬ 
ites  subjected  to  line-nose  impact.  Different  ply  orientations  and  various  thicknesses 
of  the  specimens  were  selected  for  the  tests.  Table  3.1  lists  the  ply  orientations  and 
the  geometries  of  the  specimens.  During  testing,  different  weights  and  velocities  of 
the  impactors  were  used  as  additional  test  parameters. 

T300/976  Graphite/Epoxy  prepregs  were  selected  to  fabricate  specimen  panels. 
An  autoclave  was  used  to  cure  the  panels.  The  dimensions  of  the  panels  were  about 
25.4  cm  wide  by  25.4  cm  long.  After  curing,  each  panel  was  sliced  with  a  diamond- 
coated  saw  into  three  specimens  of  7.62  cm  x  15.24  cm  All  the  specimens  were  X-rayed 
before  testing  to  evaluate  any  internal,  pre-existing  damage  caused  by  either  curing 
or  cutting. 

At  first,  a  widely  used  cure  cycle,  as  shown  in  Figure  3.3  (A),  recommended  by 
the  manufacturer  was  adopted  to  cure  the  [06/  ±  45^/905],  and  [08/906/08]  panels. 
However,  when  X-rayed,  it  was  found  that  the  [0e/  ±  454/90s],  specimens  contained 
significant  amounts  of  internal  matrix  cracks  in  the  0  and  90  degree  layers.  Further¬ 
more,  the  [08/906/08|  panels  were  completely  debonded  along  the  0  and  90  interfaces 
after  being  cured.  Hence,  it  seemed  the  cure  cycle  was  the  culprit;  the  heating  and 
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Ply  Orientation 

Thickness  (h) 
(mm) 

Span  length  (L) 
(cm) 

Width  (W) 
(cm) 

[07/902/07] 

2.30 

10.0 

7.6 

[06/904/06] 

2.30 

10.0 

7.6 

[04/±454/904/±454/04] 

2.88 

10.0 

7.6 

[03/±454/906/±454/03] 

2.88 

10.0 

7.6 

[06/±454/90io/±454/06] 

4.31 

7.6 

7.6 

[04/904/04/904/04] 

2.88 

10.0 

7.6 

[903/03/903] 

1.30 

10.0 

7.6 

[03/903/03/903/03/903/03] 

3.02 

10.0 

7.6 

Clamped  Area 
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cooling  rates  were  too  high  to  allow  the  matrix  to  relax  from  a  viscous  state  to  a 
solid  state,  which  caused  excessively  high  residual  thermal  stresses  to  build  up  inside 
the  laminates.  Accordingly,  a  new  cure  cycle  with  a  slightly  lower  cure  temperature 
was  chosen  with  slower  heating  and  cooling  rates  as  shown  in  Figure  3.3  (B).  No  pre¬ 
matrix  cracks  were  found  in  any  of  the  panels  cured  under  the  new  cycle.  Overall, 
more  than  100  specimens  were  tested  during  this  investigation.  The  results  from  the 
tests  are  summarized  in  the  following  section. 

3.4  EXPERIMENTAL  RESULTS 

3.4.1  Matrix  Cracks  and  Delaminations 
[06/904/06] 

Figure  3.4  shows  a  typical  schematic  of  the  impact  damage  pattern  in  [06/904/06] 
laminated  composites  impacted  by  a  line-nose  impactor.  Basically,  damage  appeared 
in  the  three  possible  locations,  each  of  which  consisted  of  a  matrix  crack  and  a  delam¬ 
ination;  one  near  the  center  region  and  the  other  two  close  to  the  clamped  boundaries. 
Because  of  line-loading  impact,  the  resulting  impact  damage,  matrix  cracks  and  de¬ 
laminations  are  clearly  visible  from  the  sides  of  each  specimen.  The  damage  near  the 
center  of  the  specimens  is  an  instant  result  from  the  impact  and  is  the  focus  of  the 
study. 

A  side  view  photograph  of  the  damage  near  the  center  region  of  a  lifesize  specimen 
of  [06/904/06]  after  impact  is  shown  in  Figure  3.5.  Clearly,  there  is  a  matrix  crack 
aligned  about  45  degrees  from  the  impact  direction  leading  to  delaminations  along 
the  upper  and  lower  interfaces  between  the  0  degree  and  90  degree  plies.  This  crack 
was  not  located  directly  beneath  the  impacted  area,  but  at  some  distance  away  and 
it  always  appeared  near  the  same  location  in  all  the  tested  [06/904/06]  specimens. 

It  was  discovered  that  if  these  matrix  cracks  did  not  appear  in  the  impacted 
specimens,  then  there  was  no  delamination.  Once  damage  occurred,  delamination 
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(B)  Modified  Cure  Cycle 


Figure  3.3:  Cure  cycles  used  in  the  test.  (/I):  the  original  manufacures’  recommended 
cure  cycles.  ( B ):  the  modified  cure  cycle. 
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Figure  3.5:  Photographs  of  an  undamaged  and  a  damaged  [06 /904 / Oe]  composite. 
Top:  a  side-view  photograph  of  a  life-size  damaged  specimen.  Center:  X-radiograph 
of  the  specimen  before  impact.  Bottom:  X-radiographs  of  the  specimen  after  impact. 


was  always  very  extensive  and  expanded  across  the  width  of  the  specimen.  The 
delamination  along  the  upper  interface  of  the  90  degree  layers  propagated  toward 
the  center  of  the  impacted  area,  while  the  one  located  at  the  bottom  interface  grew 
significantly  away  from  the  location  of  the  critical  crack  and  toward  the  boundary. 
The  delaminations  seemed  to  initiate  from  the  matrix  crack.  The  results  strongly 
suggested  that  these  matrix  cracks  were  the  initial  failure  mode  of  impact  damage  in 
laminated  composites.  Delaminations  were  initiated  once  the  cracks  had  propagated 
to  the  interfaces  between  the  plies  containing  different  ply  orientations.  These  types 
of  matrix  cracks  are  referred  to  as  the  “critical  matrix  cracks”  hereafter  in  this  paper. 

Figure  3.5  also  shows  the  X-radiographs  of  a  [0e/904/06]  specimen  before  and 
after  impact.  Clearly,  before  impact  there  was  no  damage  in  the  specimen.  However, 
after  impact  the  X-radiograph  shows  significant  amounts  of  internal  damage  in  the 
specimen.  Not  only  was  the  delamination  the  apparent  failure  mode,  but  extensive 
micro-cracks  were  also  found.  To  clarify,  an  enlarged  view  of  a  similar  damage  pattern 
near  the  center  region  is  presented  schematically  in  Figure  3.4.  Clearly,  these  micro¬ 
cracks  were  not  seen  by  the  naked  eye,  even  with  the  aid  of  a  binocular  microscope 
with  20  times  of  magnification,  nor  were  they  revealed  with  a  C-scan  machine.  They 
were,  however,  detectable  from  a  X-radiograph,  indicated  by  straight,  parallel  white 
lines.  Apparently,  these  micro-cracks  were  generated  along  with  delamination  growth. 

Theoretically,  a  pair  of  critical  matrix  cracks  should  exist  near  the  impacted  area 
if  the  specimen  was  symmetric  with  respect  to  the  loading  and  if  the  internal  flaws 
in  the  laminate  were  uniformly  distributed.  However,  in  reality,  only  a  few  tests  were 
found  to  produce  double  critical  cracks  near  the  center  region.  The  white  double 
curve  edges  shown  in  the  X-radiograph  of  the  damaged  specimen  are  the  boundaries 
of  the  contact  area  between  the  upper  and  lower  interfaces  of  the  delamination  in  the 
laminate.  Delaminations  shown  in  the  side-view  photograph  of  the  figure  actually 
were  equal  when  the  specimen  was  viewed  on  either  side.  In  most  of  the  case  studies, 
delaminations  propagated  throughout  the  width  of  the  specimens  quite  uniformly. 
The  delamination  fronts  near  the  free  edges  were  found  to  have  propagated  slightly 
further  than  those  in  the  inner  areas.  However,  it  was  found  that  the  energy  of  the 
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impactor  required  to  produce  the  initial  impact  damage  was  nearly  linearly  propor¬ 
tional  to  the  specimen  width,  indicating  that  the  effect  of  the  free  edge  was  minimal 
for  this  study,  since  the  actual  delamination  size  was  not  the  major  concern  of  the 
study. 

[O7/9O2/O7] 

The  damage  pattern  that  appeared  in  [O6/9O4/O6]  specimens  was  seemingly  repro¬ 
duced  in  the  [O7/9O2/O7]  specimens.  However,  for  a  given  mass  of  the  impactor,  a 
higher  velocity  was  required  to  cause  damage  in  the  [07/902/07]  specimens  than  for 
the  [O6/9O4/O6]  laminates.  As  a  result,  once  damage  occurred  it  was  more  violent  in 
the  [O7/9O2/O7]  composites  than  in  the  [0e/904 /Oe]  specimens. 

Figure  3.6  shows  three  photographs  of  a  tested  specimen.  The  top  figure  shows 
a  side-view  of  a  specimen  after  impact.  Again,  there  existed  a  critical  matrix  crack 
in  the  90  degree  layer  near  the  impacted  area.  An  extensive  delamination  was  found 
between  the  0  and  90  degree  layer  interfaces  initiating  from  the  critical  crack.  The 
delamination  propagated  along  the  lower  and  upper  interfaces  from  the  tips  of  the 
critical  matrix  crack. 

The  results  of  the  X-radiographs  of  the  specimen  before  and  after  impact  are  also 
presented  in  Figure  3.6.  Clearly,  no  damage  was  found  before  testing,  but  significant 
damage,  including  micro-cracks  and  delaminations,  appeared  in  the  specimen  after 
impact.  In  addition  to  the  micro-cracks  in  the  90  degree  layers,  the  0  degree  layers 
also  contained  a  few  matrix  cracks  which  apparently  enhanced  the  extent  of  the 
delaminations.  Accordingly,  the  size  of  the  delaminations  was  much  more  non-uniform 
in  shape  than  in  the  [06/904/06]  composites. 

[04/  ±  454/902],  and  [03/  ±  454/903], 

Laminates  containing  ±45  degree  angle  plies  were  also  evaluated  in  the  study. 
Again,  all  the  damaged  specimens  contained  critical  matrix  cracks.  Figures  3.7  and 
3.8  show  photographs  of  typical  [03/  ±454/903],  and  [04 /  ±454/902],  specimens  after 
impact,  respectively.  Surprisingly,  in  addition  to  the  critical  cracks  located  near  the 
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Sideview 


Figure  3.6:  Photographs  of  an  undamaged  and  a  damaged  [O7/9O2/O7]  compos¬ 
ite.  Top:  a  side-view  photograph  of  a  life-size  damaged  specimen.  Center:  an 
X-radiograph  of  the  specimen  before  impact.  Bottom:  an  X-radiograph  of  the  speci¬ 
men  after  impact. 
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center  and  the  boundaries  of  the  specimens,  cracks  were  also  found  between  the  central 
impacted  region  and  the  boundaries,  as  shown  in  Figures  3.7  and  3.8.  All  the  critical 
matrix  cracks  appeared  in  the  90  degree  layer  and  led  to  extensive  delaminations. 
Delamination  occurred  along  the  interface  between  the  ±45  degree  and  90  degree 
layers.  By  examining  the  X-radiographs  in  Figures  3.7  and  3.8,  it  appears  that  besides 
the  micro-cracks  in  the  90  degree  layers,  there  are  also  micro-cracks  found  in  the  ±45 
degree  layers  within  the  delamination  region.  Apparently,  these  micro-cracks  in  the 
±45  degree  layers  were  generated  along  with  the  delamination  propagation  as  were 
those  in  90  degree  layers. 

[04/904/04/904/04] 

All  the  specimens  considered  above  contained  only  one  group  of  90  degree  plies  lo¬ 
cated  in  middle  plane  of  the  specimens.  Tests  were  also  performed  on  [04/904/04/904/04] 
specimens  which  contained  two  90  degree  ply  groups  located  off  the  middle  plane.  It 
was  intended  to  evaluate  the  effect  of  the  location  and  number  of  90  degree  ply  groups 
on  impact  damage.  Typical  results  generated  from  the  tests  are  shown  in  Figure  3.9. 

Near  the  centrally  impacted  region,  a  critical  matrix  was  found  in  each  90  degree 
ply  group,  leading  to  delamination  along  each  interface  between  the  0  degree  and  90 
degree  ply  groups.  Sometimes  a  double  critical  matrix  crack  symmetrically  appearing 
from  the  impacted  area  was  produced,  as  shown  in  Figure  3.10.  The  delamination 
on  the  lower  interface  propagated  much  more  extensively  than  the  one  on  the  upper 
interface.  These  two  delaminations  were  separated  by  the  0  degree  layer  group  located 
in  the  middle  plane. 

[03/903/03/903/03/903/03] 

Like  the  above  tested  specimens,  the  results  from  the  [03 / 903 /03/ 903 /03 / 903 / 03] 
specimens,  which  contained  three  90  degreee  ply  groups,  showed  multiple  critical 
matrix  cracks  in  the  90  degree  layer  groups.  A  photograph  of  a  typical  test  result  is 
shown  in  Figure  3.11.  Once  again,  these  critical  cracks  led  to  multiple  delaminations 
which  were  individually  constrained  by  the  0  degree  ply  groups.  One  interesting 
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After 


Figure  3.7:  Photographs  of  a  damaged  [O3 /  ±  454/903],  composite.  Top:  a  side-view 
photograph  of  a  life-size  damaged  specimen.  Bottom:  an  X-radiograph  of  the  speci¬ 
men  after  impact. 
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After 


Figure  3.8:  Photographs  of  a  damaged  [04/  ±454/902]s  composite.  Top:  a  side-view- 
photograph  of  a  life-size  damaged  specimen.  Bottom:  an  X-radiograph  of  the  speci¬ 
men  after  impact. 
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Figure  3.9:  A  schematic  of  a  typical  damage  pattern  of  [0n/90m/0„/90m/0„]  compos¬ 
ites. 


After 


Figure  3.10:  Photographs  of  a  damaged  [O4/9O4/O4/9O4/O4]  composite.  Top:  a 
side-view  photograph  of  a  life-size  damaged  specimen.  Bottom:  an  X-radiograph 
of  the  specimen  after  impact. 
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note  worth  mentioning  is  that  the  critical  matrix  cracks  occurring  near  the  central 
impacted  region  were  much  closer  to  the  center  line  than  those  of  the  other  tested 
specimens  with  different  ply  orientations. 

[9O3/O3/9O3] 

Tests  were  also  performed  on  [9O3/O3/9O3]  laminates  containing  90  degree  ply 
groups  on  the  outer  layers.  Intuitively,  it  was  expected  that  the  90  degree  ply  groups 
are  prone  to  produce  intraply  matrix  cracking,  because  of  excessive  out-of-plane  bend¬ 
ing  stress  due  to  impact.  The  intuition  was  confirmed  by  the  experiments.  A  photo¬ 
graph  of  a  lifesize  tested  specimen  of  [9O3/O3/9O3]  is  shown  in  Figure  3.12.  Clearly, 
a  matrix  crack  initiated  directly  beneath  the  impacted  area  on  the  bottom  90  degree 
ply  group  and  produced  a  delamination  along  the  interface  between  the  90  and  0 
degree  ply  groups.  It  is  noted  that  the  matrix  crack  is  aligned  vertically  (parallel  to 
the  loading  direction)  and  different  from  the  other  embedded  matrix  cracks,  which 
were  all  inclined  by  an  angle  with  respect  to  the  loading  direction,  found  in  previously 
tested  specimens  with  the  0  degree  ply  groups  on  the  outer  surfaces. 


3.4.2  Impact  Energy  Threshold 

One  unique  feature  of  the  2-dimensional  impact  results  is  that  once  damage  occurs, 
delamination  is  always  very  extensive.  Figure  3.13  presents  the  estimated  delami¬ 
nation  size  in  specimens  of  various  ply  orientations  as  a  function  of  impact  energy. 
There  apparently  exists  an  impact  energy  threshold  beyond  which  damage  occurs. 
For  instance,  the  energy  threshold  of  [06/902]a  composites  is  about  100  J/m  for  a 
given  impactor  of  1.142  kg/m,  beyond  which  significant  delaminations  are  produced. 
No  damage,  including  matrix  cracks  and  delaminations,  was  found  in  any  of  the  spec¬ 
imens  tested  below  that  energy  level.  Similar  results  are  also  shown  in  Figures  3.13 
for  other  ply  orientations.  It  is  noted  that  all  the  dimensions  above  are  normalized 
by  width,  based  on  the  two-dimensional  approach. 

It  was  strongly  implicated  that  the  impact  energy  threshold  was  associated  with 
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Figure  3.11:  Photographs  of  a  damaged  [O3/9O3/O3/9O3/O3/9O3/O33  composite.  Top: 
a  side-view  photograph  of  a  life-size  damaged  specimen.  Bottom:  an  X-radiograph  of 
the  specimen  after  impact. 


Figure  3.12:  Photographs  of  a  damaged  [9O3/O3/9O3]  composite.  Top:  a  side-view- 
photograph  of  a  life-size  damaged  specimen.  Bottom:  an  X-radiograph  of  the  dam¬ 
aged  specimen. 
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Normalized  Dclamination  Area  (D/LH) 


Figure  3.13:  Illustration  of  the  relationships  between  impact  damage  size  and  the 
impact  energy  for  different  ply  orientations. 
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the  energy  required  to  initiate  the  impact  damage.  The  existence  of  an  impact  en¬ 
ergy  threshold  for  point-loading  cases  also  has  been  reported  by  other  investigators 
[25,34,47].  Additionally,  it  appears  that  the  impact  energy  threshold  also  strongly 
depends  on  the  ply  orientation,  seen  in  Figure  3.13.  The  energy  threshold  may  vary 
for  different  ply  orientations.  This  information  on  the  impact  energy  threshold  could 
be  useful  in  design  as  a  guideline  for  selecting  appropriate  composites  and  for  de¬ 
termining  the  configuration  and  layups  of  composite  structures  in  order  to  sustain 
impact. 

3.4.3  Pre-existing  Micro-cracks  Due  to  Thermal  Residual  Stresses 

Due  to  the  thermal  expansion  coefficient  mismatch,  manufacturing  produces  signifi¬ 
cant  thermal  residual  stresses  in  laminated  composites  [55-58].  The  amount  of  resid¬ 
ual  stress  depends  strongly  on  the  degree  of  the  thermal  coefficient  mismatch,  the 
ply  orientation,  and  the  cure  cycle.  For  [06 /  ±454/905],  layups  cured  under  the  stan¬ 
dard  cure  cycle  given  in  Figure  3.2,  residual  stresses  exceeded  the  in  situ  transverse 
strength  of  the  90  degree  plies  in  the  laminate  and  caused  micro-matrix  cracks,  as 
shown  in  Figure  3.14.  However,  for  [04/  ±  454/902]s  and  [03/  ±  454/903],  laminates 
cured  under  the  new  cycle,  no  damage  was  found  after  curing. 

For  the  specimens  with  pre-matrix  cracks,  the  impact  energy  required  to  initiate 
damage  was  found  to  be  substantially  lower  than  for  those  laminates  with  similar 
layups  but  without  pre-matrix  cracks.  Figure  3.13  shows  the  reduction  of  the  impact 
energy  threshold  in  [06/±454/905],  laminates  compared  to  those  of  the  [03/±454/903]5 
and  [04/  ±  454/902]a  laminates  without  pre-matrix  cracks.  Although  the  actual  ply 
stacking  sequences  are  slightly  different  among  those  compared,  the  reduction  of  the 
energy  threshold  of  [06/  ±  454/905],  composites  is  very  substantial  and  can  be  pri¬ 
marily  attributed  to  the  existence  of  the  pre-matrix  cracks.  Figure  3.14  also  shows 
an  X-radiograph  of  a  [0e/  ±  454/905],  specimen  with  pre-matrix  cracks  after  impact. 
Surprisingly,  no  additional  pre-matrix  cracks  were  found  after  impact,  and  the  de¬ 
laminations  were  initiated  from  one  of  the  existing  pre-matrix  cracks  and  propagated 
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Figure  3.14:  Photographs  of  an  impacted  specimen  of  a  [06 /  ±  454 / 905]a  compos¬ 
ite.  Top:  a  side-view  photograph  of  a  life-size  damaged  specimen.  Bottom:  an 
X-radiograph  of  the  specimen  before  and  after  impact.  A  considerable  number  of 
internal  preexisting  cracks  were  found  before  impact. 
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into  interfaces  in  a  similar  pattern. 

Apparently,  the  critical  matrix  cracks  in  [0«/  ±454/905]4  composites  were  created 
from  existing  pre-matrix  cracks  and  hence  required  much  less  energy  than  normally 
required  to  produce  the  critical  crack.  Accordingly,  the  impact  energy  threshold  of 
the  pre-cracked  composites  is  considerably  lower  than  for  composites  without  flaws. 
Hence,  these  results  strongly  indicates  that  the  impact  energy  threshold  is  related  to 
the  energy  needed  to  initiate  the  critical  matrix  cracks. 

It  is  worth  noting  that  even  if  laminates  contain  no  pre- matrix  cracks,  there  always 
exists  a  considerable  amount  of  residual  stresses  in  the  materials  which  could  have 
a  significant  effect  on  impact  damage.  For  the  materials  containing  higher  residual 
stresses,  it  is  expected  that  the  materials  would  have  a  much  lower  damage  resistance 
to  impact  than  those  containing  lower  residual  stresses.  Hence,  the  energy  threshold 
of  a  composite  also  depends  strongly  upon  the  amount  of  residual  stresses  left  in  the 
material  due  to  manufacturing. 

3.5  CONCLUSION 

An  experimental  investigation  was  performed  to  study  impact  damage  in  T300/976 
Graphite/epoxy  laminated  composites.  Based  on  the  test  results,  the  following  re¬ 
marks  can  be  made: 

1.  Matrix  cracking  initiates  impact  damage. 

2.  Delamination  is  always  accompanied  by  a  “critical”  matrix  crack. 

3.  Considerable  micro-cracks  can  be  generated  along  with  delamination  growth 
during  impact. 

4.  There  exists  an  impact  energy  threshold  above  which  impact  damage  occurs. 

5.  Pre-existing  micro-cracks  induced  by  thermal  stresses  can  substantially  reduce 
impact  resistance  of  composites. 
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6.  Stacking  sequence  significantly  affects  impact  resistance  of  composites. 
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Chapter  4 


2-DIMENSIONAL  LINE-LOADING  IMPACT  ANALYSIS 


4.1  INTRODUCTION 

In  order  to  further  substantiate  the  experimental  findings  and  to  understand  the 
basic  mechanics  causing  low-velocity  impact  damage,  an  analytical  investigation  was 
subsequently  performed.  An  analytical  model  was  developed  for  predicting  the  impact 
damage  resulting  from  the  tests.  The  numerical  simulations  based  on  the  model  were 
also  carried  out  to  understand  the  physics  of  the  impact  damage  mechanics. 

The  analytical  model  described  consists  of  a  stress  analysis  for  calculating  tran¬ 
sient  dynamic  stresses,  strains  and  deformations  during  impact,  and  a  failure  analysis 
for  predicting  impact  damage.  Numerical  simulations  of  the  test  conditions  based 
on  the  analysis  are  also  presented.  The  effects  of  geometry  and  ply  orientation  of 
the  composites,  and  the  velocity  and  mass  of  the  impactor  on  impact  damage  are 
evaluated  as  well.  Based  on  the  study,  the  model  combined  with  the  test  results 
provide  essential  information  on  the  mechanisms  and  mechanics  of  impact  damage  in 
laminated  composites  due  to  low-velocity  impact. 
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4.2  STATEMENT  OF  PROBLEM 


The  problem  is  described  as  follows:  consider  a  symmetric  laminated  composite  plate 
clamped  along  two  parallel  edges  and  free  on  the  other  two  edges.  A  line-nose  impactor 
strikes  uniformly  across  the  width  on  the  center  line  of  the  specimen  (see  Figure  3.2). 
For  a  given  mass  of  the  impactor,  it  is  desired  to  determine: 

1.  The  initial  failure  mode. 

2.  The  location  of  the  initial  damage. 

3.  The  velocity  required  to  cause  the  initial  impact  damage. 

4.  The  initiation  of  delamination  and  micro-cracks. 

5.  The  direction  of  delamination  growth. 

6.  The  essential  parameters  governing  the  impact  damage. 

4.3  ANALYSIS 

During  the  investigation,  an  analytical  model  was  developed  for  predicting  the  impact 
damage  in  laminated  composites  due  to  line- loading  impact.  Since  the  load  distribu¬ 
tion  across  the  width  is  assumed  to  be  uniform  at  any  instant  time  and  the  boundary 
conditions  are  symmetric  with  respect  to  the  loading  condition,  the  impact  response 
of  the  plate  is  assumed  to  be  two-dimensional  and  independent  of  the  x2  axis,  the 
direction  along  the  width  (see  Figure  3.2).  Ignoring  the  effect  of  free  edges,  the  plane 
strain  condition  is  adopted  for  the  development  of  a  two-dimensional  transient  dy¬ 
namic  finite  element  analysis  for  the  investigation.  The  finite  element  analysis  is  used 
to  calculate  transient  dynamic  stresses,  strains,  and  deformations  inside  the  laminate 
during  impact. 


34 


The  experiments  presented  in  the  previous  Chapter  strongly  indicate  that  impact 
damage  is  initiated  by  matrix  cracking  and  that  there  exist  critical  matrix  cracks 
from  which  delaminations  and  micro-cracks  initiate.  In  order  to  substantiate  the 
findings,  the  proposed  model  focuses  on  predicting  the  initial  failure  mode  and  the 
critical  matrix  cracks,  and  determining  the  mechanics  governing  the  initiation  of 
delaminations  and  micro-cracks.  Thus,  in  the  model,  a  matrix  failure  criterion  is 
adopted  and  modified  for  predicting  the  initial  failure  mode.  A  post-failure  analysis 
after  the  occurrence  of  the  initial  damage  is  also  performed  to  determine  the  effect  of 
the  initial  damage  on  the  initiation  of  delamination  and  micro-cracks. 

4.3.1  Stress  Analysis 

Based  on  plane  strain  condition,  the  equilibrium  equation  at  an  instant  time  in  a 
variational  form  can  be  expressed  as  [59,60] 

f  pu,tt6u,6v  -f  f  (jjBtjdv—  [  T,6u,da  =  0  {  (4.1) 

Jv  Jv  J$a  [  j  =  1. 3, 5 

where  u,  and  u,,«  are  the  displacements  and  the  accelerations  (tt1>t(  =  d2ujdt 2), 
respectively.  T i  is  the  contact  force  distribution  during  impact.  a}  and  are  the 
stresses  and  strains  due  to  the  mechanical  loading  in  contracted  notations,  i.e., 


{<71,03,(75}  —  {<711,033,013}  (4.2) 

{«lif2'«5}  =  {« 11,  f33i  £13}  (4.3) 

For  laminated  composites,  the  material  properties  may  vary  from  layer  to  layer 
throughout  the  thickness,  depending  upon  the  laminate  stacking  sequence.  Hence. 
Eq.  (4.1)  can  be  rewritten  on  a  layer-by-layer  basis  as  follows: 
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Y  J  |y  npnUiM8nUi  dv  +  J  ncrjSnejdv  dh  =  J  Trfiiida  (4.4) 

where  iV  is  the  total  number  of  the  layers  of  the  laminate  considered.  The  forward 
superscript  n  for  all  the  variables  indicates  the  corresponding  layer  within  the  lami¬ 
nate. 

Thus,  the  mechanical  stresses  in  the  n-th  layer  are  related  to  the  mechanical 
strains  of  the  corresponding  layer  through  the  following  equation 


%=  n[D]k}nek  (4.5) 

where  is  the  stiffness  matrix  of  the  n-th  ply  of  the  composites  based  on  the 

plane  strain  assumption.  A  complete  expression  of  the  stiffness  matrix  "[Dj*]  for  a 
single  layer  based  on  the  plane  strain  assumption  is  given  in  the  Appendix  A. 

In  order  to  solve  Eq.  (4.1),  the  contact  load  distribution  7’,,  between  the  impactor 
and  the  composite,  must  first  be  known.  In  this  investigation,  the  Hertzian  contact 
law  was  adopted  to  simulate  the  contact  load  distribution.  Because  of  the  use  of  a 
cylindrical  line-nose  impactor,  the  Hertzian  contact  law  has  a  considerably  different 
expression  from  the  one  that  is  most  commonly  used  for  a  spherical  point-nose  im¬ 
pactor  [61].  Accordingly,  the  contact  load  distribution  T,  (=  /)  for  a  line-loading 
impact  is  related  to  the  indentation  depth  a  (=  6s  —  6c  =  the  change  in  the  dis¬ 
tance  between  the  center  of  the  impactor  nose  and  the  mid-plane  of  the  plate)  by  the 
following  expression  [4] 


o  =  fix.  +  Xp){l  -  ln  l/r(x.  +  XP)]}  when  6S  >  6C  ^ 

0  =  /  when  6s  <  6q 

where  6$  is  the  impactor  displacement,  6c  is  the  plate  displacement  measured  at  the 
center  of  the  mid-plate  of  the  laminate,  which  is  opposite  to  the  impactor’s  surface. 
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and  Xa  and  Xp  are  the  constants  defined  as  [61] 


1  -  v] 

Xa  =  — vr- 

ir  Et 

1 

XP  7 rEyy 


(4.7) 

(4-8) 


and  r,  u„,  and  E»  are  the  local  radius,  the  Poisson’s  ratio,  and  the  Young’s  modulus 
of  the  impactor,  respectively.  Evv  (=  Ezz )  is  the  modulus  of  elasticity  of  the  impacted 
composite  ply  in  the  direction  transverse  to  the  fibers. 

A  transient  dynamic  finite  element  program  was  developed  based  on  the  analysis. 
An  eight  node  isoparametric  element  was  used,  and  a  typical  finite  element  mesh 
generated  from  the  program  is  shown  in  Figure  4.1. 


4.3.2  Finite  Element  Formulation 

The  stresses  and  strains  in  the  laminate  were  calculated  by  a  two-dimensional,  tran¬ 
sient  finite  element  method.  The  equations  used  in  the  numerical  calculations  are 
presented  in  this  section. 
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Figure  4.2:  A  schematic  of  the  8-node  serendipity  quadrilateral  element. 


Finite  Element  Model 

In  the  finite  element  analysis,  an  8-node  serendipity  quadrilateral  element  was 
adopted  (Figure  4.2).  The  displacements  at  any  point  in  the  element  (ui,ti3)  can  be 
expressed  as 


8 

«,  =  NrV  9=1,3  (4.9) 

r=l 

where  tt,r  are  the  displacements  of  nodal  point  r  in  the  x,  direction.  The  shape 
functions  of  the  8-node  quadrilateral  elements  Nr(r  =  1  ~  8)  are  given  as  follows 
[59,60] 
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(4.10) 


£  and  (  are  natural  coordinates  for  each  element  (Figure  4.2)  whose  values  vary  from 
—  1  to  +1.  The  coordinates  (xl5x3)  of  any  point  inside  the  element  are  related  to  the 
natural  coordinate  through  the  shape  functions 


x,  =  £Nrx,r  g=l,3  (4.11) 

r=l 

where  xqr  are  the  coordinates  of  the  nodal  point  r. 

From  Eq.  (4.9)  the  strains  at  any  point  in  the  laminate  can  be  written  as 


{€lli  c33i2Ci3}T  —  {tin  e33i  7l3}T 

=  DBr]{uIr,U3r}T  (4‘12) 

r=l 

where  cn  and  €33  are  the  normal  strains,  £13  are  the  components  of  shear  strain  tensor, 
and  713  are  the  engineering  shear  strains.  The  symbols  {  }  and  [  ]  represent  vectors 
and  matrices,  respectively.  The  superscript  T  means  transposition  of  a  vector  or  a 
matrix.  The  [Br]  matrices  is  defined  as 
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[B,]  = 


0 


r  =  1  ~  8 


(4.13) 


‘Nr>1 
0  Nr,3 
.Nr,3  Nril. 

The  symbols  Nr,t  represent  the  derivatives  of  the  shape  functions  with  respect  to 
x,  (N r,t  =  dNr/dxi).  The  stresses  are  related  to  the  strains  by  the  Eq.  (4.5). 

Governing  Equations 

By  neglecting  damping,  the  governing  equation  can  be  written  as  [60] 


|A/)(</)  +  [/>]{</)  =  {F)  (4.14) 

where  [A/]  and  [A']  are  the  mass  and  stiffness  matrices,  and  {A},  {d}  and  {d}  are 
the  force,  displacement,  and  acceleration  vectors,  respectively.  The  mass  matrix  [A/] 
is  the  sum  of  the  element  mass  matrix  [A/e] 


N'l 

[A/]  =  £[A/1  (4.15) 

e=l 

where  Nei  is  the  number  of  total  elements.  The  element  mass  matrix  is  the  sum  of 
the  element  mass  submatrices 


[AH-ED^W  (4-16) 

p=  1  r=l 

In  general,  the  element  mass  submatrices  are  [59] 

{M'lr  =  JjBpf  p[Br}da  p, r  =  1  ~  8  (4.17) 

where  At  is  the  area  of  element,  and  p  is  the  ply  density.  For  computational  con¬ 
venience,  the  mass  matrix  lumping  procedure  is  used.  In  this  method  the  diagonal 
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terms  of  the  consistent  mass  matrix  is  scaled  so  that  the  mass  is  preserved.  In  the 
case  of  the  rectangular  8-node  serendipity  elements  used  in  our  analysis,  the  element 
mass  matrix  becomes  [59] 
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The  stiffness  matrix  [I\  ]  is  the  sum  of  element  stiffness  matrices  [A'e] 


N'l 

[A'l  =  E(A"I 


e=  1 


The  element  stiffness  matrix  is  the  sum  of  the  element  stiffness  submatrices 


8  8 


!A"I  =  E  EIA-W 


p=l  r= 1 


where 


[A'eW=  /  [Bp]T[D][Br]da  p,  r  =  1  ~  8 

JA' 

The  elasticity  matrix  [D]  in  Eq.  (4.21)  is  same  as  in  Eq.  (4.5). 


(4.18) 


(4.19) 


(4.20) 


(4.21) 
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In  the  present  problem  with  the  absence  of  the  body  forces,  the  contact  force 
caused  by  the  impactor  is  the  only  external  load  to  be  considerd.  The  contact  force 
vector  {F}  is  defined  as 


{ F }  =  f{U)  (4.22) 

where  {{7}  is  a  unit  vector  whose  component  is  —1  in  the  direction  of  the  contact 
force.  All  other  components  of  {U}  are  zero  and  /  is  the  magnitude  of  the  contact 
force  which  is  a  function  of  time. 

At  time  t  +  At,  Eq.  (4.14)  is  written  as 


[M]{d}t+At  +  [K]{d)t+At  =  {F}t+At  (4.23) 

where  the  superscript  refers  to  time.  The  Newmark  method  will  be  employed  to  obtain 
the  solutions  to  this  equation.  Accordingly,  the  velocity  and  acceleration  vectors  at 
time  t  +  At  are  written  as 


{d}'+A<  =  {dY  +  (1  -  A)A t{d}(  +  XAt{d) 


t  +  At 


and 


{d} 


t+At  _ 


8A  t2 


{d} 


t+At 


1  i  -  —  8  • 

— = — {tf}'  — —  {</}*  - 1 — -{dY 
8 At2  ‘  ]  8 A tw  8  1 


(4.24) 

(4.25) 


The  parameters  8  and  A  are  constants  whose  values  depend  on  the  finite  difference 
scheme  used  in  the  calculations.  The  constant-average-acceleration  method  was  em¬ 
ployed,  which  is  implicit  and  unconditionally  stable.  For  this  method  8  is  1/4  and  A 
is  1/2.  Although  the  velocity  vector  is  not  required  in  Eq.  (4.24),  it  is  presented  here 
because  it  will  be  needed  subsequently.  By  substituting  Eq.  (4.25)  into  Eq.  (4.23). 
we  obtain 


[/;']{</},+A'  =  {f},+A( 


(4.26) 
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where  [ K\  is  the  effective  stiffness  matrix,  and  {F}<+A(  is  the  effective  force  vector. 
These  parameters  are  defined  as 


[A'l  =  J^IM]  +  (A)  (4.27) 

and 

{F}t+At  =  {HY  +  {F}‘+A‘  (4.28) 

where  {HY  is  the  following  vector 

w  -  I"l(ss?  w  +  jsW  +  «•») 

Referring  to  Eq.  (4.26),  it  is  noted  that  the  displacements,  velocities,  and  accel¬ 
erations  at  time  t  are  known  at  every  point  inside  the  plate.  Hence,  the  unknowns 
in  Eq.  (4.26)  are  the  displacement  vector  {<f}<+A‘  and  the  force  vector  {F}‘+At.  To 
determine  these  two  variables,  an  additional  expression  is  needed.  This  expression  is 
described  below. 

Solution  Procedure 

Combining  Eq.  (4.26)  with  Eq.  (4.28),  results  in 


[I<]{dY+At  =  {HY  +  {F},+At  (4.30) 

The  displacement  vector  {d},+At  is  expressed  as  the  sum  of  the  displacements  due  to 
the  force  {H }‘  and  the  contact  force  {F},+A< 


{d}t+Al  =  {d}‘„+A‘  +  {(f}^+A< 

Eq.  (4.30)  and  Eq.  (4.31)  give 


(4-31) 
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[/i-]((<f}£4'  +  {</}?4')  =  {H}’  +  {FT** 


(4.32) 


From  Eq.  (4.32), 


=  {H)‘  (4.33) 

and 

[K]{d}£At  =  {F}t+At  (4.34) 

In  order  to  determine  the  displacement  vectors  {d}//  and  {d}/r  in  Eq.  (4.31  )-(4.32), 
the  two  force  vectors  {//}  and  {F}  at  every  nodal  point  must  be  specified.  At  each 
nodal  point  { H }  is  given  by  Eq.  (4.29)  and  {F}  equals  to  zero  except  at  the  point 
of  impact  where  {F}  =  f{U)  [see  Eq.  (4.22)].  The  solutions  of  Eq.  (4.30)-(4.31 ) 
are  obtained  in  two  steps.  First,  the  forces  {H}  and  the  displacements  {d}u  are 
calculated  from  Eq.  (4.29)  and  (4.33).  This  calculation  is  straightforward,  since  the 
mass  matrix  [M],  the  time  step  At,  as  well  as  the  displacement  {d},  the  velocity  {d}, 
and  acceleration  {d}  at  time  t  are  all  known.  Second,  the  contact  force  vector  {F} 
and  the  displacement  {d}f  are  calculated  as  follows.  At  time  t  +  At,  Eq.  (4.22)  is 
written  as 


{Fy+*‘  =  ft+At{U} 

where  ft+At  is  the  contact  force  at  time  t  +  At.  Eq.  (4.34)  and  (4.35)  yield 

{K){d}F+*'  =  f+At{V} 

For  a  unit  contact  force  ( ft+At  =  1),  Eq.  (4.3G)  becomes 


(4.35) 


(4.36) 
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[A'lW"  =  W ) 


(4.37) 


where  {d}^^  is  the  displacement  caused  by  the  unit  contact  force.  For  a  given 
effective  stiffness  matrix  the  displacement  vector  can  be  calculated  from  Eq. 

(4.37).  From  Eqs.  (4.37)  and  (4.38),  it  can  be  seen  that  and  {d}£/"A<  can  be 

related  by 


(4.38) 


Eqs.  (4.31)  and  (4.38)  give 

{d}'+A‘  =  {d}£A‘  +  /,+*‘{d}(/+A‘  (4-39) 

In  Eq.  (4.39),  the  unknowns  are  the  displacement  vector  {d}  and  the  scalar  contact 
force  /  at  time  f  +  At.  During  impact  the  contact  force  can  expressed  as 

at+A(  _  ft+At(Xa  +  xp)|  1  -  In  [ft+£itr(x,  +  x>)]|  when  Ss  >  Sc 
0  =  fl+At  when  6$  <  6C 

Xs,Xp>&Si$Ci  and  r  are  known  constants  defined  previously,  a  is  the  indentation  depth 
which  varies  with  time.  At  time  t  +  A t  this  depth  is 

af+A<  =  £+*'  -  S£At  (4.41) 

S'c**'  is  the  displacement  of  the  center  of  the  midsurface  of  the  plate  in  the  direction 
of  the  impact.  With  the  use  of  Eq.  (4.39),  6c  can  be  expressed  as 

f>c“  =  S'cT  -  /‘+A,*cif‘  (4.42) 
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where  6tQAt,  S'c^1  and  Seif1  are  the  magnitudes  of  the  displacements  at  the  center 
point  of  the  mid-surface  in  the  direction  of  impact  at  time  t  +  A t.  S s+At  is  the 
position  of  the  center  point  of  the  impactor.  At  time  t  +  At,  the  magnitude  of 
is  determined  by  Newton’s  second  law  (See  Appendix  B) 

4+A'  = 

(4.43) 

=  +  vtAt  -  A-(2f<  +  /,+a,)A(2 

where  m  is  impactor  mass  and  iq  is  the  impactor  velocity  at  time  t.  By  combining 
Eqs.  (4.40)-(4.43),  the  contact  force  can  be  expressed  as  follows: 


o  =  /,+A'(X-.  +  X>){l-ln[/'+4V(x,  +  Xp)]} 

-«!i-v1Al  +  gL(2/'  +  /'«>)Af2 

(4.44) 

+  (6'cT  +  !'+*%¥' )  when  6S  >  Sc 

0  =  ft+At  when  6$  <  Sc 

Newton-Raphson  method  was  ed  to  calculate  the  contact  force  ft+At  in  Eq. 

(4.44).  Once  the  value  of  ft+At  is  known,  the  displacement  vector  {<f}  at  time  t  +  At 
is  calculated  from  Eq.  (4.39),  and  the  velocity  and  acceleration  vector  at  time  t  +  At 
are  calculated  from  Eqs.  (4.24)  and  (4.25).  From  the  known  displacements,  the  strains 
are  calculated  from  Eq.  (4.12),  and,  from  the  known  strains,  the  stresses  are  obtained 
from  Eq.  (4.5).  The  aforementioned  procedures  are  repeated  at  each  time  step. 

Once  the  stresses  are  calculated,  initial  damage  of  the  composites  can  be  deter¬ 
mined  by  substituting  the  stresses  into  an  appropriately  selected  failure  criterion. 
Accordingly,  the  accuracy  of  the  prediction  of  damage  relies  on  the  proper  selection 
of  the  failure  criterion  and  the  accuracy  of  information  regarding  the  stresses  in  the 
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laminate  and  the  ply  strength  distributions  within  the  composites.  It  is  noted  that  the 
stresses  thus  calculated  from  the  analysis  are  associated  directly  with  the  mechanical 
loading.  There  also  exists  a  certain  amount  of  thermal  residual  stresses  in  composites 
resulting  from  fabrication  which  could  inherently  weaken  the  strength  of  the  material 
and  its  ability  to  sustain  the  impact.  The  thickness  and  ply  orientation  of  the  com¬ 
posite  could  also  affect  the  ply  strength  distributions  within  the  composite  [62-65]. 
Hence,  in  order  to  accurately  predict  the  initial  impact  damage,  the  in  situ  strength 
distribution  of  composites  as  well  as  the  transient  dynamic  stress  distributions  should 
be  considered  in  conjunction  with  the  failure  criterion. 

4.3.3  Failure  Analysis 

The  major  focus  of  the  failure  analysis  is  to  numerically  simulate  the  impact 
resulting  from  the  line-loading  impact.  Hence,  the  proposed  failure  analysis 
trates  on  predicting  the  occurrence  and  location  of  the  critical  matrix  cracks 
initiation  of  delamination  and  micro-cracks. 

Initial  Damage 

A  matrix  failure  criterion  is  utilized  to  determine  the  initial  failure  mode  of  impact 
damage,  the  critical  matrix  cracks.  In  this  study,  a  three-dimensional  matrix  failure 
criterion  originally  proposed  by  Hashin  [54]  was  adopted  and  modified  for  predicting 
the  matrix  cracks.  The  initial  impact  damage  is  predicted  whenever  the  calculated 
stresses  satisfy  the  criterion  and  the  corresponding  location  is  associated  with  the 
place  where  the  matrix  crack  occurs. 

Here,  based  on  line-loading  impact,  there  are  only  three  major  stress  components 
that  can  contribute  to  initial  matrix  cracking  in  the  90  degree  layers  in  the  plane 
strain  condition,  as  shown  schematically  in  Figure  4.3.  These  are  the  interlaminar 
shear  stress  013  (=  ayz),  in-plane  tensile  stress  <ru  (=  ayv),  and  out-of-plane  normal 
stress  033  (=  <7„).  However,  it  will  be  shown  later  that  the  out-of-plane  normal  stress 
<733  is  found  to  be  very  small  in  comparison  with  the  other  two  stress  components 


damage 
concen- 
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during  the  entire  impact  event.  Hence,  the  matrix  failure  criterion  can  be  simplified 
as  (<ryy  >  0), 


rcVy\2  (nZvA2  2  (  €m  >  1  Failure 

\  nYt  )  \  ”5,  /  M  |ew<l  No  failure 


(4.45) 


where  the  subscript  of  x  and  y  are  the  local  coordinates  of  the  n  —  th  layer  parallel  and 
normal  to  the  fiber  directions,  respectively,  and  2  is  the  out-of-plane  direction.  nYt  and 
n5,  are  the  in  situ  ply  transverse  tensile  strength  and  interlaminar  transverse  shear 
strength  within  the  n-th  ply  of  laminate  under  consideration,  respectively  [62.G3]. 
"<7yI  and  ”<fyy  are  the  averaged  interlaminar  transverse  shear  stress  and  the  averaged 
in-plane  transverse  tensile  stress,  respectively,  within  the  n-th  ply,  which  can  be 
expressed  as 

"Vyz  =  T~  f  Oyzdz  (4.46) 

Jtn- 1 

and 


avv 


dz 


(4.47) 


where  tn  and  <„_i  are  the  upper  and  lower  interfaces  of  the  n-th  ply  or  ply  group  in 
the  laminate  and  hn  is  the  thickness  of  the  ply  or  ply  group. 


Based  on  Eq.  (4.45),  the  determination  of  the  initial  impact  damage  strongly 
depends  on  the  in  situ  strength  of  the  material  within  the  laminate.  It  is  well  known 
that  the  strengths  of  a  single  ply  within  a  laminate  can  be  substantially  different 
from  the  strengths  measured  directly  from  a  unidirectional  composite  [62-65];  the 
difference  has  been  attributed  to  the  thermal  residual  stresses,  ply  orientation,  and 
the  thickness  of  the  laminate.  However,  the  actual  strength  distributions  of  a  laminate 
as  a  function  of  these  parameters  are  still  not  well  established. 

Here,  an  empirical  relationship  of  the  ply  transverse  tensile  strength  within  a 
laminate  as  a  function  of  the  laminate  thickness  and  stacking  sequence  was  adopted, 
which  was  proposed  by  Chang  and  Lessard  [62].  The  relationship  can  be  expressed 
as 


48 


Figure  4.3:  The  stress  components  based  on  the  plane  strain  condition  contributing 
to  transverse  matrix  cracking. 
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where  Yt  is  the  in  situ  ply  transverse  tensile  strength  of  the  laminate,  A 6  is  the 
maximum  ply  angle  change  between  the  ply  or  the  ply  group  under  consideration 
and  its  neighboring  plies,  and  Mc  is  the  total  number  of  consecutive  (clustered)  ply 
groups  that  are  considered.  Y°  is  the  transverse  tensile  strength  of  a  [90„]3  composite 
(n  >  6).  The  parameters  of  A  and  B  are  constants  which  can  be  determined  from 
experiment  [62]. 

Note  that  the  strength  distribution  thus  determined  corresponds  to  the  mechanical 
load  that  causes  the  first  matrix  crack  in  the  laminate.  The  distribution  depends  on 
the  manufacturing  process,  ply  orientation  and  thickness  of  the  laminate.  Hence,  the 
value  of  Yt  determined  from  Eq.  (4.48)  can  be  used  in  Eq.  (4.45)  for  predicting  initial 
impact  damage. 

However,  the  interlaminar  shear  strength  distribution  within  a  laminated  com¬ 
posite  has  not  yet  been  fully  studied.  It  is  still  not  understood  how  the  strength 
distribution  is  affected  by  the  stacking  sequence  and  the  laminate  thickness,  but  such 
distribution  has  been  frequently  assumed  to  be  in  the  same  order  of  magnitude  as  the 
ply  shear  strength  distribution  within  a  laminate.  Accordingly,  here,  it  was  assumed 
that  the  interlaminar  shear  strength  normal  to  the  fiber  direction  is  equivalent  to  in 
situ  ply  shear  strength  S  within  the  laminate.  The  distribution  of  the  in  situ  ply 
shear  strength  S  has  been  observed,  experimentally,  to  exhibit  similar  behavior  as 
the  in  situ  transverse  tensile  strength  Y(  [62].  Hence,  the  expression  for  the  in  situ 
ply  shear  strength  distribution  proposed  by  Chang  and  Lessard  [62]  was  adopted  for 
estimating  the  interlaminar  shear  strength  distribution  normal  to  the  fiber  direction 
within  a  laminate.  The  expression  similar  to  Eq.  (4.48)  can  be  written  as 


where  S  is  the  in  situ  ply  shear  strength  of  a  laminate  and  S°  is  the  ply  shear  strength 
measured  from  a  unidirectional  composite  with  more  than  eight  layers.  Again,  C  and 
D  are  the  material  constants  which  can  be  determined  from  experiment  [62], 

Delamination  Growth  and  Micro-cracks 

Once  the  critical  matrix  crack  is  predicted,  the  next  step  is  to  discover  if  such 
a  crack  could  initiate  delaminations  and  micro-cracks  in  a  composite  during  impact. 
Accordingly,  a  post-failure  analysis  was  performed  to  simulate  the  subsequent  re¬ 
sponse  of  the  composite  containing  the  initial  damage.  The  analysis  was  executed 
by  reducing  the  material  stiffnesses  within  the  damaged  element  where  the  critical 
matrix  crack  was  predicted.  For  matrix  cracking,  the  material  properties  within  the 
damaged  layer  were  reduced  as  follows  [62,66,67] 
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where  Ex  and  Ez  are  the  tensile  moduli  of  the  n-th  layer  along  the  fibers  and  normal 
to  the  middle  plane,  respectively.  Gxz  and  Gxy  are  the  out-of-plane  shear  moduli  of 
the  n-th  layer.  The  expression  for  the  undamaged  material  properties  of  n[D]  is  given 
in  Appendix  A. 

Eq.  (4.50)  indicates  that  the  damaged  element  can  not  sustain  any  additional 
transverse  tensile  stress  and  out-of-plane  shear  stress  due  to  the  presence  of  the  crack. 
This  reduction  is  reasonable  within  the  neighborhood  of  the  crack  for  a  quantitative 
study,  although  a  dynamic  fracture  analysis  may  result  in  more  qualitative  findings. 
The  stresses  and  strains  were  then  recalculated  using  Eq.  (4.5)  at  the  same  instant 
in  time  as  the  occurrence  of  the  critical  matrix  cracks.  A  Singly,  the  redistributed 
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stress  distributions  near  the  neighborhood  of  the  critical  matrix  cracks  could  be  stud¬ 
ied. 


4.4  VERIFICATION 

During  the  investigation,  a  computer  code  designated  as  “2DIMPACT”  was  devel¬ 
oped  based  on  the  proposed  analysis.  In  order  to  verify  the  calculations  of  stresses 
and  strains  from  the  code,  comparisons  were  made  between  the  existing  analytical  so¬ 
lutions  and  the  calculated  results  from  the  code.  The  stress  distributions  through  the 
thickness  of  a  composite  plate  subjected  to  quasi-static  cylindrical  bendings  were  cal¬ 
culated  and  compared  with  the  analytical  solutions  obtained  by  Pagano  [68].  Figure 
4.4  shows  the  calculated  results  from  the  code  and  the  analytical  solutions  [68].  Ex¬ 
cellent  agreements  were  obtained.  Furthermore,  the  dynamic  response  of  a  cantilever 
beam  subjected  to  a  periodic  loading  was  determined  from  the  code  and  compared 
with  the  analytical  solutions  in  Figure  4.5  [69].  Again,  the  calculations  agreed  with 
the  analysis  very  well. 

Unfortunately,  no  analytical  solutions  of  the  transient  dynamic  stress  distributions 
due  to  line-loading  impact  were  available  in  the  literature  for  comparison.  Hence, 
the  solutions  of  a  beam  subjected  to  a  line  loading  impact  based  on  the  code  were 
compared  to  the  solutions  of  the  beam  subjected  to  a  point  loading  impact.  It  is 
reasonable  to  consider  that  if  the  width  of  the  beam  is  narrow  enough,  the  response 
of  the  beam  and  the  contact  force  distribution  would  be  consistent  in  both  cases.  The 
solutions  resulting  from  the  point-loading  impact  were  calculated  from  the  ‘"IMPACT" 
code  which  was  developed  previously  by  Wu  et  al  [13,16]. 

Figure  4.6  shows  the  central  deflections  of  the  beam  for  both  loading  cases  at 
a  given  impact  energy.  The  two  curves  coincide  with  each  other  very  well.  The 
contact  force  distributions  as  a  function  of  duration  of  time  for  both  cases  are  also 
shown  in  Figure  4.6.  Again,  based  on  two  completely  different  contact  algorithms, 
the  two  calculated  contact  force  distributions  were  very  consistent  with  each  other 
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Figure  4.4:  Comparison  between  the  calculated  and  the  analytical  solutions  for  a 
laminated  composite  panel  subjected  to  cylindrical  bending.  The  analytical  solutions 
are  taken  from  [68]. 
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Figure  4.5:  Comparison  between  the  calculated  and  the  analytical  solutions  for  a 
laminated  composite  panel  subjected  to  periodic  loading.  The  analytical  solutions 
are  taken  from  [69]. 
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and  exhibited  a  very  similar  pattern.  Accordingly,  it  was  believed  that  the  code 
could  provide  accurate  information  about  the  stresses  and  deformations  inside  the 
structures  during  impact.  Hence,  numerical  solutions  were  generated  from  the  code 
to  compare  with  the  test  results,  and  the  results  of  the  comparison  are  presented  in 
the  next  section. 


4.5  COMPARISON 

Numerical  simulations  of  the  specimens  under  the  test  condition  are  performed  by 
using  the  code.  The  material  properties  of  T300/976  Graphite/Epoxy  used  in  the 
calculations  are  summarized  in  Table  4.1.  The  results  of  the  calculations  compared 
to  the  experimental  data  are  presented  as  follows: 


4.5.1  Initial  Impact  Damage 

To  predict  the  initial  matrix  cracking,  the  test  results  of  the  [06/904/06]  and  [03/  ± 
454/903]s  specimens  are  first  simulated  by  the  code.  The  weight  and  velocity  of  the 
impactor  used  in  the  calculations  are  the  same  as  the  test  conditions.  Figures  4.7  and 
4.8  show  the  peaks  of  the  values  of  ew  in  90  degree  layers  (matrix  failure  criterion) 
as  a  function  of  position  for  [06/904/06]  and  [03/  ±  454/903],  specimens  under  the 
test  condition  respectively.  For  >  1,  matrix  cracking  is  predicted.  For  instance. 
Figure  4.7  indicates  that  there  are  two  possible  locations  in  90  degree  layers  for  matrix- 
cracking  at  the  given  test  condition  (M  =  1.142  kg/m  and  V  =  14.5  m/sec);  one  is 
close  to  the  central  impacted  area  at  t  =  351  //sec,  and  the  other  is  near  the  clamping 
areas  at  t  =  423  //sec.  The  central  crack  occurs  much  earlier  than  the  others  near  the 
boundaries.  Apparently,  the  former  is  the  direct  result  of  impact  from  the  impactor 
and  the  latter  is  primarily  due  to  bending  and  stretching  of  the  plates  resulting  from 
the  constraints  of  the  boundaries. 

Surprisingly,  a  rather  different  strength  ratio  distribution  ( e a/ )  is  obtained  for  the 
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Figure  4.6:  Comparison  of  the  calculations  between  the  “2DIMPACT”  and  ‘‘IM¬ 
PACT”  [13,16]  for  a  composite  beam  subjected  to  a  line-loading  impact  or  a 
point-loading  impact,  respectively.  Top:  the  central  deflection  of  the  beam  for  both 
loading  cases  at  a  given  impact  energy.  Bottom:  the  contact  distribution  as  a  function 
of  the  duration  of  time  for  both  loading  cases. 
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Moduli 


Symbol  (unit) 


In-plane  longitudinal  modulus 

Exx  (Gpa) 

156 

In-plane  transverse  modulus 

Eyy  (Gpa) 

9.09 

In-plane  shear  modulus 

Gxy  (Gpa) 

6.96 

Out-of-plane  shear  modulus 

Gyz  (Gpa) 

3.24 

In-plane  Poisson's  ratio 

nxy 

0.228 

Out-of-plane  Poisson's  ratio 

nyz 

0.400 

Density 

r  (kg/m3) 

1540 

Strength 

Symbol  (unit) 

Longitudinal  tension 

XT(Mpa) 

1520 

Longitudinal  compression 

Xc(Mpa) 

1590 

Transverse  tension 

Yt  (Mpa) 

45 

Transverse  compression 

Yc(Mpa) 

252 

Ply  longitudinal  shear 

S°  (Mpa) 

105 

Strength  parameters  (From  Eq.  4.48) 

A* 

1.3 

Strength  parameters  (From  Eq.  4.48) 

B* 

0.7 

Strength  parameters  (From  Eq.  4.49) 

C* 

2.0 

Strength  parameters  (From  Eq.  4.49) 

D* 

1.0 

Impactor 

Symbol  (unit) 

Modulus 

Es  (Gpa) 

207 

Poisson’s  ratio 

ns 

0.3 

Nose  radius 

r  (mm) 

1.5 

*  Data  taken  from  [  62  ] 

Table  4.1:  Material  properties  of  T300/976  used  in  the  calculations. 
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Strength  Ratio 


Figure  4.7:  The  calculated  maximum  strength  ratios  of  e a/  (failure  criterion)  for  a 
[0e/904 /Oe]  specimen  occurring  at  two  separate  instants. 
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Strength  Ratio 


Figure  4.8:  The  calculated  maximum  strength  ratios  of  e^f  (failure  criterion)  for  a 
[03 /  ±  454/903],  specimen  occurring  at  three  separate  instants. 


[03/  ±  454/903],  composites.  No  matrix  cracking  is  predicted  in  ±45  degree  layers 
under  the  given  loading  condition.  Besides  the  two  possible  locations  for  matrix 
cracking  similar  to  those  of  [O6/9O4/O6]  composites,  an  additional  peak  of  the  strength 
ratio  ( e\f  =  1)  occurs  first  at  time  t  =  321  //sec  and  appears  at  a  quarter  of  the  test 
span  of  the  specimen,  measured  from  the  boundary  (about  half  the  distance  from  the 
center  of  the  specimen  to  the  boundary).  The  location  of  the  critical  matrix  cracks 
corresponds  to  the  position  where  e\f  reaches  the  unity.  Figure  4.9  shows  the  locations 
of  the  predicted  critical  cracks  compared  to  the  experimental  data  for  [06/904/06]  and 
[03 /  ±454/903],  composites.  The  experimental  results  in  Figure  4.9  were  taken  from 
an  average  of  the  measured  data  of  three  to  four  tested  specimens. 

It  is  worth  noting  that  the  critical  matrix  crack  near  the  central  region  of  the 
specimens  does  not  occur  directly  beneath  the  impacted  area,  but  a  distance  from 
it.  In  order  to  understand  the  phenomenon,  the  distributions  of  the  stresses  near 
the  central  impacted  region  were  carefully  studied.  It  was  found,  for  instance,  in 
[06/904/06]  specimens  that  the  interlaminar  shear  stress  and  the  in-plane  transverse 
tensile  stress  in  the  90  degree  layers  reach  the  peaks  at  t  =  351  //sec  near  the  location 
where  the  critical  matrix  occurs,  but  the  out-of-plane  normal  tensile  stress  cr33  is 
always  negligibly  small  compared  to  the  others  and  decreases  rapidly  once  it  is  away 
from  the  impacted  area  (see  Figure  4.10).  Apparently,  the  interlaminar  shear  stress 
and  the  in-plane  transverse  tensile  stress  are  the  dominant  stresses  causing  the  critical 
matrix  cracking  during  impact.  Thus,  the  effect  of  the  out-of-plane  normal  tensile 
stress  on  the  initiation  of  impact  damage  is  negligible.  As  a  result,  without  out- 
of-plane  normal  stress,  delamination  initiation  which  requires  the  combination  of 
out-of-plane  normal  stress  and  interlaminar  shear  stress  would  most  unlikely  occur 
earlier  than  the  matrix  failure,  indicating  that  the  critical  matrix  cracking  is  the  initial 
failure  mode. 

Furthermore,  the  distributions  of  interlaminar  shear  stress  across  the  thickness  at 
various  locations  were  also  calculated  and  studied.  Figure  4.11  shows  the  interlaminar 
shear  stress  distributions  at  six  different  cross  sections  near  the  central  impacted  areas 
at  time  t  =  351  //sec.  Interestingly,  the  distributions  varied  considerably  from  cross 
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Initial  Matrix  Crack  Locations 


Figure  4.9:  Comparison  of  the  locations  of  the  initial  matrix  cracks  between  the 
averaged  data  and  the  predictions. 
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Figure  4.10:  The  instantaneous  stress  distributions  near  the  center  impacted  area  of 
a  [06/90* /06]  specimen. 
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section  to  cross  section,  and  the  shear  distribution  in  90  degree  layer  reached  the 
maximum  at  the  location  near  the  critical  matrix  crack.  The  magnitude  of  such 
shear  distribution  is  considerably  higher  than  the  one  obtained  statically  by  a  line¬ 
loading  force  under  the  same  out-of-plane  deflection  as  shown  in  Figure  4.12.  The 
difference  can  be  attributed  to  the  transient  dynamic  effect.  Clearly,  the  interlaminar 
shear  stress  distribution  is  strongly  influenced  by  the  impact  loading  and,  hence,  is  a 
function  of  the  velocity  and  the  mass  of  the  impactor  as  well. 

The  distribution  of  the  transverse  normal  stress  at  time  t  =  351  //sec,  given  in 
Figure  4.10,  was  also  compared  with  a  static  solution  based  on  the  same  amount  of 
deflection  corresponding  to  the  dynamic  case  at  t  =  351  //sec.  The  comparison  is 
shown  in  Figure  4.13.  Apparently,  the  distribution  is  not  very  sensitive  to  dynamic 
effect  and  is  primarily  attributed  to  bending  and  stretching.  Hence,  it  can  be  con¬ 
cluded  that  the  initial  impact  damage  is  governed  by  transient  dynamic  interlaminar 
shear  stress  and  in-plane  normal  stress.  The  location  and  the  amount  of  energy  re¬ 
quired  to  cause  the  initial  impact  damage  strongly  depend  on  the  velocity  and  mass 
of  the  impactor  and  the  geometry  and  boundary  condition  of  the  specimen. 

A  very  consistent  prediction  of  the  location  of  the  critical  matrix  cracks  for  each 
possible  90  degree  layer  in  the  laminates  was  also  obtained  for  specimens  with  other 
ply  orientations  and  geometries.  The  results  of  the  comparison  are  summarized  in 
Figure  4.14.  All  of  these  critical  matrix  cracks  located  near  the  central  impacted 
region  of  the  specimens  were  embedded  inside  the  laminates  and  inclined  of  an  angle 
about  45  degrees  from  the  loading  direction.  The  inclination  of  the  matrix  cracks  is 
due  to  the  combination  of  the  interlaminar  shear  stress  and  the  transverse  normal 
stress  in  the  90  degree  layers. 

However,  numerical  simulations  were  also  performed  on  [903/03/903]  specimens  in 
which  the  critical  cracks  were  found  experimentally  in  the  outer  90  degree  layers.  At 
a  given  mass  and  velocity  of  the  impactor,  the  code  predicted  that  a  critical  matrix 
crack  would  be  located  at  the  center  of  the  specimen  in  the  bottom  90  degree  layers. 
The  prediction  also  coincided  with  the  experiment  very  well.  It  was  found  that  the 
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Figure  4.11:  The  calculated  interlaminar  shear  stress  distribution  across  the  thickness 
at  six  different  cross  sections  near  the  central  impacted  areas  of  a  [0e/904 /06]  specimen 
at  the  time  corresponding  to  the  first  critical  matrix  crack. 
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Figure  4.12:  Comparison  of  the  interlaminar  shear  stress  distribution  along  the  lower 
90/0  interface  of  a  [06/904/06]  composite  specimen  between  the  static  and  the  tran¬ 
sient  dynamic  calculations.  The  static  load  was  applied  to  generate  an  equal  amount 
of  lateral  deflection  as  that  of  the  dynamic  case  at  the  compared  instant. 
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Figure  4.13:  Comparison  of  the  in-plane  tensile  stress  distribution  in  the  90°  layers 
along  the  lower  90/0  interface  of  a  [06/9004/06]  composite  specimen  between  the  static 
and  the  transient  dynamic  calculations.  The  static  load  was  applied  to  generate  an 
equal  amount  of  lateral  deflection  as  that  of  the  dynamic  case  at  the  compared  instant. 


Figure  4.14:  Comparison  of  the  locations  of  the  initial  matrix  cracks  between  the 
averaged  data  and  the  predictions. 


67 


critical  matrix  crack  is  dominated  by  the  in-plane  transverse  tensile  stress  <7yy,  and 
hence,  the  crack  surface  was  normal  to  the  loading  direction. 

Therefore,  based  on  the  results  of  the  calculations,  it  can  be  concluded  that  the 
critical  matrix  crack  is  the  initial  failure  mode  and  the  interlaminar  shear  stress/strength 
and  in-plane  transverse  tensile  stress/strength  are  very  critical  for  initiating  impact 
damage.  The  inclined  cracks  are  dominated  by  the  interlaminar  shear  stress,  and  the 
vertical  cracks  are  governed  by  the  transverse  normal  tensile  stress. 

4.5.2  Impact  Energy  Threshold 

It  is  observed  from  the  previous  experiments  (Chapter  3)  that  there  exists  an  impact 
energy  threshold  beyond  which  impact  damage  occurs,  but  below  which  no  damage 
can  be  found.  The  same  phenomenon  is  also  reported  for  point-loading  impact  br¬ 
others  [25,34,47].  It  is  strongly  indicated  from  the  previous  experimental  observation 
that  the  impact  energy  threshold  is  related  to  the  energy  required  to  cause  the  first 
critical  matrix  crack. 

Thus,  numerical  calculations  were  performed  to  calculate  the  energy  required  to 
initiate  the  first  critical  matrix  crack  for  all  the  test  configurations.  For  a  given  mass 
of  the  impactor,  the  numerical  simulation  was  performed  at  various  velocities  until 
the  first  critical  matrix  crack  wac  predicted.  The  impact  energy  threshold  is  defined 
as  one  half  of  the  mass  multiplied  by  the  square  of  the  velocity  corresponding  to 
the  initial  matrix  crack  (1/2 MV2).  The  predicted  impact  energy  thresholds  for  vari¬ 
ous  ply  orientations  and  thickness  compared  to  the  test  data  are  presented  in  Table 
4.2.  Overall,  the  predictions  agree  with  the  test  data  very  well.  The  predictions  are, 
in  general,  consistently  conservative  compared  to  the  experiments.  The  predictions 
further  strongly  confirm  the  experimental  observation  that  the  impact  energy  thresh¬ 
old  is  associated  with  the  energy  required  to  initiate  the  first  critical  matrix  crack. 
Delamination  will  proceed  once  the  critical  matrix  is  produced. 
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Ply  Orientation 

Mass 

(kg/m) 

Test 

Velocity 

(m/s) 

Predicted 

Velocity 

(m/s) 

Predicted 

Energy 

(J/m) 

[07/902/07] 

1.14 

19±1.5 

20 

228 

[0^904/06] 

1.14 

14±1 

15.5 

137 

[0V±454/904/±454/Gi] 

1.57 

15±1 

15 

177 

[Q3/±454/906/±454/G3] 

1.57 

17±1 

14 

154 

[04/904/04/90V04] 

0.80 

25±2 

20 

157 

[03/903/03/903/03/903/03] 

4.40 

13±2 

11 

266 

Table  4.2:  Prediction  of  impact  velocity  threshold. 


Apparently,  by  examining  Table  4.2,  it  can  be  found  that  the  impact  energy  thresh¬ 
old  is  also  strongly  affected  by  ply  orientation  and  laminate  thickness.  Ply  orientation 
(especially  the  stacking  sequence)  seems  to  have  more  effect  on  impact  energy  thresh¬ 
old  than  thickness.  For  instance,  the  impact  energy  threshold  of  [04/904/04/904/04] 
laminates  is  about  15  perce;  t  higher  than  that  of  [06 /904 / 06]  laminates  with  the  lami¬ 
nate  thickness  increased  by  25  percent,  from  16  layers  to  20  layers.  However,  compar¬ 
ing  the  energy  thresholds  between  [03/903/ 03/903/ 03/903/03]  and  [04/904/04/904/04] 
laminates  with  only  one  layer  difference  in  thickness,  the  energy  threshold  increased 
by  40  percent  from  [04/904/04/904/04]  to  [O3/9O3/O3/9O3/O3/9O3/O3].  Accordingly, 
ply  orientation  and  thickness  are  also  very  important  factors,  especially  the  former, 
on  the  design  of  composite  structures  which  may  be  subjected  to  impact. 
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4.6.3  Delamination  and  Micro-cracks 


Once  a  critical  matrix  is  predicted,  the  material  stiffnesses  within  the  element  con¬ 
taining  the  crack  are  reduced  according  to  the  material  degradation  rules,  and  the 
stress  distributions  near  the  damaged  element  are  recalculated.  Figure  4.15  shows 
the  redistributed  out-of-plane  normal  stresses  along  the  ply  interfaces  of  a  [06/904/06] 
laminate  before  and  after  matrix  cracking.  Directly  before  the  occurrence  of  the 
matrix  crack,  the  out-of-plane  normal  stress  is  near  zero  around  the  damaged  loca¬ 
tion.  However,  two  peak  out-of-plane  normal  stresses  with  opposite  signs  are  found 
immediately  adjacent  to  the  damaged  element  along  each  upper  and  lower  interface 
between  the  0  and  90  degree  layers.  Furthermore,  the  distribution  of  interlaminar 
shear  stress  increases  in  the  neighborhood  of  the  damaged  area  and  reduces  sharply 
to  zero  w’ithin  the  damaged  area,  as  shown  in  Figure  4.16. 

Along  the  lower  interface,  the  concentrated  normal  stress  is  positive  leading  away 
from  the  impacted  area,  and  becomes  negative  toward  the  impacted  area.  How¬ 
ever,  along  the  upper  interface,  the  situation  of  the  stress  distribution  is  completely 
reversed.  This  result  indicates  that  two  delaminations  can  be  initiated  by  the  out- 
of-plane  normal  tensile  stresses  along  the  upper  and  lower  interfaces  resulting  from 
mode  I  fracture.  The  interlaminar  shear  stress  near  the  damage  area  consequently  can 
enhance  the  growth  of  the  delamination  once  it  propagates.  The  delamination  along 
the  lower  interface  propagates  away  from  the  impacted  region,  but  the  delamination 
along  the  upper  interface  moves  toward  the  center  of  the  specimen.  The  results  from 
the  prediction  can  be  illustrated  by  the  drawing  shown  in  Figure  4.17.  This  prediction 
is  very  consistent  with  the  experimental  findings  as  shown  in  Figure  3.5,  which  is  a 
side  view  of  a  lifesize  photograph  of  a  [06/904/0e]  specimen  after  impact. 

Similar  results  are  also  obtained  for  the  distribution  of  in-plane  transverse  tensile 
stresses  of  a  [06/904/0g]  laminate  before  and  after  the  occurrence  of  the  critical  matrix 
crack.  The  results  of  calculations  are  shown  in  Figure  4.18.  Again,  there  appears  to 
be  two  peak  in-plane  transverse  stresses  along  the  same  upper  and  lower  interfaces 
as  the  out-of-plane  normal  tensile  stresses  near  the  damaged  area  immediately  after 
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Figure  4.15:  Comparison  of  the  instantaneous  out-of-plane  normal  stress  distributions 
along  the  upper  and  lower  90/0  interfaces  in  the  90°  plies  of  a  [06/9004/06j  specimen 
before  and  after  impact. 
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Figure  4.16:  Comparison  of  the  instantaneous  interlaminar  shear  stress  distributions 
along  the  upper  and  lower  90/0  interfaces  in  the  90°  plies  of  a  [0e/9004 / 0e]  specimen 
before  and  after  impact. 
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Figure  4.17:  A  schematic  description  of  the  impact  damage  mechanism  resulting  from 
an  internal  matrix  crack. 


matrix  cracking.  This  can  imply  that  the  resulting  highly  concentrated  in-plane 
normal  tensile  stress  can  cause  micro-cracks  in  the  90  degree  layers  accompanying  the 
delamination  growth.  These  micro-cracks  along  with  delaminations  were  detected 
by  X-ray  from  the  experiments.  As  a  comparison,  a  typical  sample  of  a  radiograph 
of  these  micro-cracks  of  a  [06/904/06]  specimen  is  shown  in  Figure  3.5.  It  is  worth 
noting  that  the  extent  of  the  micro-cracks  is  quite  consistent  with  the  size  of  the 
delaminations. 

Similar  results  for  evaluating  the  initiation  of  delamination  and  micro-cracks  are 
also  obtained  for  [9O3/O3/9O3]  specimens.  The  distributions  of  the  out-of-plane  normal 
stress  and  transverse  tensile  stress  before  and  after  the  occurrence  of  the  critical 
matrix  crack  at  the  center  of  the  specimen  along  the  90  and  0  degree  ply  interface  are 
presented  in  Figure  4.19.  Again,  stress  concentrations  arise  near  the  crack  tip  along 
the  interface  and  trigger  the  initiation  of  delamination  and  micro-cracks.  Accordingly, 
based  on  the  predictions  and  the  experimental  observations  from  the  previous  chapter, 
the  physical  processes  of  impact  damage  mechanisms  in  laminated  composites  due  to 
low  velocity  impact  can  be  illustrated  schematically  in  sequential  steps  as  shown  in 
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Figure  4.18:  Comparison  of  the  instantaneous  in-plane  transverse  normal  stress  dis¬ 
tributions  along  the  upper  and  lower  90/0  interfaces  in  the  90°  plies  of  a  [Oe/900.,/06] 
specimen  before  and  after  impact. 


74 


Figure  4.20. 


4.6  CONCLUSION 

An  analytical  investigation  was  performed  to  study  impact  damage  in  laminated 
composites  as  a  result  of  line-loading  impact.  An  analysis  was  developed  for  predicting 
the  impact  damage  and  for  understanding  impact  damage  mechanisms  and  mechanics 
of  laminated  composites.  The  results  of  the  predictions  agreed  very  well  with  the 
experiments.  Based  on  the  study,  the  following  remarks  can  be  made: 

1.  Matrix  cracking  is  the  initial  failure  mode  of  impact  damage  of  laminated  com¬ 
posites. 

2.  Delamination  and  micro-cracks  due  to  impact  are  initiated  by  the  matrix  cracks 
(the  critical  matrix  cracks). 

3.  Impact  energy  threshold  is  associated  with  the  energy  required  to  initiate  the 
first  critical  matrix  crack. 

4.  Pre-existing  micro-cracks  can  significantly  reduce  the  impact  resistance  of  com¬ 
posites. 

5.  Interlaminar  shear  stresses/strength  and  in-plane  tensile  stress/strength  are  the 
dominant  factors  causing  the  critical  matrix  cracks. 

6.  Out-of-plane  normal  stress  (or  Mode  I)  and  interlaminar  shear  stress  (or  Mode 
II)  are  important  for  determination  of  delamination  growth  after  the  critical 
matrix  is  introduced. 

7.  In-plane  transverse  tensile  stress  and  interlaminar  shear  stress  near  the  critical 
matrix  cracks  produce  micro-cracks  as  delamination  propagates. 
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Figure  4.19:  Comparison  of  the  instantaneous  out-of- plane,  interlaminar  shear  stress, 
and  in-plane  normal  stress  distributions  along  the  lower  0/90  interfaces  in  the  90° 
plies  of  a  {9O3/O3/9O3]  specimen  before  and  after  impact. 
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Figure  4.20:  A  schematic  description  of  two  basic  impact  damage  growth  mechanisms 
of  laminated  composites. 


Chapter  5 

EFFECTS  OF  IMPACT  PARAMETERS 


5.1  INTRODUCTION 

The  effects  of  laminate  configurations,  such  as  stacking  sequence  and  thickness,  and 
mass  of  the  impactor  on  the  impact  damage  of  laminated  composites  due  to  low- 
velocity  line  load  impact  were  also  studied.  The  major  focus  of  the  study  concen¬ 
trated  on  the  initial  damage  and  the  corresponding  failure  mode.  Hence,  it  was  of 
particular  interest  to  determine  the  velocity  of  the  impactor  required  to  initiate  the 
impact  damage  as  a  function  of  the  stacking  sequence  and  the  thickness  of  the  com¬ 
posites  for  different  impactors1  masses.  In  order  to  understand  the  basic  mechanics 
and  mechanisms  of  the  test  results,  numerical  simulations  based  on  an  analytical 
model  (Chapter  4)  were  also  performed.  Based  on  the  study,  the  relationships  of  the 
initial  damage  with  respect  to  the  laminate  configuration  and  impactor’s  mass  were 
established. 

5.2  EXPERIMENTS 

In  order  to  achieve  the  objectives,  laminates  made  of  T300/976  Graphite/Epoxy 
prepregs  with  different  stacking  sequences  and  thicknesses  were  deliberately  selected 
for  the  tests.  For  stacking  sequence  effect,  cross-ply  composites,  containing  about  67 
percent  of  0  degree  plies  and  33  percent  of  90  degree  plies,  with  four  different  stacking 
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Ply  Orientation 

Thickness  (h) 
(mm) 

Span  length  (L) 
(cm) 

Width  (W) 
(cm) 

Impactor  Mass 
(kg/m) 

[06/906/06] 

2.58 

10.0 

23 

4.17 

6.67 

[OV9O3/OV9O3/04] 

2.58 

10.0 

23 

2.50 

4.17 

6.67 

[03/902 /W9Q2/Da/9O2/09] 

2.58 

10.0 

23 

2.50 

4.17 

6.67 

[0z/90/02/90/0/90/0]s 

2.58 

10.0 

2.3 

4.17 

[(03/903)2/0!] 

2.16 

10.0 

23 

2.50 

4.17 

6.67 

[(03/903)3/03] 

3.02 

10.0 

2.3 

2.50 

4.17 

6.67 

[(03/903  >/03] 

3.88 

10.0 

2.3 

2.50 

4.17 

6.67 

Table  5.1:  Ply  orientations  and  geometries  of  the  test  specimens. 


sequences  were  selected:  (06/906/06],  [O4/9O3/O4/9O3/O4],  [O3/9O2/O3/9O2/O3/9O2/O3] 
and  [02/90 /02/90/0/90/0]a,  all  of  which  have  a  constant  thickness.  To  evaluate  the 
thickness  effect,  the  [(03/903)„/03]  composites  were  chosen.  The  subscript  n  indicates 
the  number  of  the  repeat  of  (O3/9O3)  in  the  laminate.  Three  different  thicknesses  were 
considered.  In  addition,  three  different  masses  were  used  for  each  selected  configura¬ 
tion.  Table  5.1  lists  the  stacking  sequences,  geometries,  and  masses  of  the  impactors 
used  in  the  experiments. 

All  the  specimens  were  cured  under  the  same  cure  cycle  used  in  Chapter  3  (see 
Figure  3.3)  and  were  cut  into  the  same  size:  10.1  cm  in  length  and  2.54  cm  in  width. 
An  X-radiograph  was  taken  of  all  the  specimens  before  impact  to  examine  the  internal 
damage  resulting  from  manufacturing  or  cutting.  No  apparent  damage  was  found  in 
any  of  the  specimens. 

Each  specimen  was  firmly  clamped  on  two  parallel  edges,  and  the  other  edges  were 
left  free,  as  shown  in  Figure  5.1.  A  specially  designed  L-shaped  aluminum  tab  was 
adhesively  bonded  to  each  clamped  end  of  the  specimen  to  prevent  any  slippage  of 
the  specimen  from  the  fixture  during  impact.  A  schematic  of  the  test  configuration  of 
the  specimen  is  shown  in  Figure  5.1.  More  than  fifty  specimens  were  tested.  For  each 
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configuration,  four  to  five  replicants  were  tested.  After  every  test,  each  specimen  was 
unloaded  from  the  fixture  and  thoroughly  inspected  by  a  binocular  microscope  with 
20  times  of  magnification  to  determine  any  damage  from  the  sides  of  the  specimen. 

Initially,  X-radiographs  were  also  frequently  taken  to  confirm  the  eye  inspection. 
After  many  comparisons,  it  was  found  that  the  initial  damage  was  so  obviously  no¬ 
ticeable  by  the  naked  eye  that  a  dual  confirmation  by  X-radiograph  was  unnecessary. 
Therefore,  the  X-radiograph  was  only  used  afterward  when  the  damage  inspected  by 
the  binocular  microscope  was  in  question.  The  undamaged  specimens  were  reloaded 
and  tested  at  a  higher  velocity,  and  then  inspected  again  for  initial  damage.  The 
procedure  was  repeated  until  initial  damage  was  detected.  Most  specimens  were  re¬ 
peatedly  tested  less  than  three  times.  It  is  noted  that  the  attention  of  the  experiment 
is  focused  on  examining  the  initial  damage,  the  mode  of  failure,  and  determining 
the  corresponding  impactor’s  energy  or  velocity  required  to  cause  such  damage  in 
composites. 


5.3  RESULTS  AND  DISCUSSIONS 

The  results  of  the  measured  impact  damage  as  a  function  of  the  energy  of  the  im- 
pactor  at  the  instant  of  impact  for  various  specimen  configurations  are  summarized 
in  Figure  5.2.  Each  box  in  the  figure  represents  a  test  configuration.  The  test  data 
clearly  shows  that  an  impact  energy  threshold  exists  for  all  the  composites  tested, 
beyond  which  significant  damage  occurs,  but  below  which  no  damage,  including  de¬ 
laminations  and  microcracks,  is  found.  The  impact  energy  threshold  is  apparently 
the  minimum  energy  required  for  causing  impact  damage.  It  has  been  illustrated 
that  the  initial  impact  damage  consisted  of  several  critical  matrix  cracks  from  which 
the  subsequent  damage,  delamination  and  miciocracks,  is  induced.  Thus,  the  impact 
energy  threshold  corresponds  to  the  energy  that  produces  the  critical  matrix  cracks. 

Accordingly,  the  initial  impact  damage  and  the  corresponding  impact  energy  are 
very  important  for  characterizing  impact  resistance  of  laminated  composites,  because 
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Figure  5.1:  Description  of  the  geometry  and  the  boundary  conditions  of  the  test 
specimens. 
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Figure  5.2:  A  summary  of  the  test  results.  The  measured  impact  damage  size  as  a 
function  of  the  impact  energy  for  various  specimen  configurations. 


they  can  lead  to  extensive  damage  in  composites.  Therefore,  the  initial  impact  damage 
and  the  energy  threshold  are  evaluated  and  analyzed  thoroughly  in  the  foilwing. 

5.3.1  Stacking  Sequence  Effect 

The  test  results  of  cross-ply  laminated  composites,  consisting  of  three  different  stack¬ 
ing  sequences,  with  a  constant  thickness  and  volume  fraction  of  0  degree  and  90  degree 
plies  are  presented  in  Figure  5.3.  Three  different  masses  of  the  impactors  are  consid¬ 
ered.  The  results  presented  in  each  box  of  the  figure  correspond  to  a  different  mass.  In 
each  box,  the  ordinate  is  the  impact  velocity  threshold  corresponding  to  the  velocity 
of  the  impactor  of  a  given  mass  required  to  cause  initial  damage.  The  horizontal  axis 
indicates  the  three  ply  orientations.  At  a  given  mass,  the  test  results  clearly  show 
that  the  stacking  sequence  significantly  affects  the  impact  velocity  as  well  as  the  im¬ 
pact  energy  required  to  initiate  the  impact  damage.  The  impact  velocity  threshold 
Vc  can  be  enhanced  by  as  much  as  30  to  50  percent  by  reordering  the  sequence  of  the 
plies  within  the  laminates  from  l06/906/06]  to  [O3/9O2/O3/9O2/O3/9O2/O3],  depending 
upon  the  impactor’s  mass.  It  is  worth  noting  that  if  the  impact  energy  (1/2AH'2)  is 
used  instead  of  the  velocity  Vc  in  the  ordinate,  the  difference  in  the  impact  energy 
threshold  among  these  laminates  is  greatly  magnified. 

Apparently,  the  velocity  (energy)  of  an  impactor  must  be  increased  in  order  to 
produce  initial  impact  damage  as  the  number  of  layers  grouped  together  with  the 
same  ply  orientation  within  a  laminate  decreases.  In  other  words,  it  is  expected 
that  the  more  uniformly  dispersed  the  ply  sequence  is  in  a  laminate,  the  higher  the 
initial  impact  damage  resistance  will  be.  However,  this  statement  is  valid  only  under 
the  condition  that  the  initial  damage  is  governed  by  critical  matrix  cracks,  because 
sharp  nose  impactors  with  a  higher  impact  velocity  could  easily  cause  surface  damage 
such  as  fiber  breakage  on  the  surface  of  the  laminates  containing  thin  ply  groups  on 
outer  surfaces.  For  example,  a  few  tests  were  performed  on  [O2/9O/O2/9O/O/9O/O], 
composites,  and  a  photograph  of  a  life-size  tested  specimen  is  shown  in  Figure  5.4. 
The  specimen,  which  has  the  same  thickness  as  [O3/9O2/O3/9O2/O3/9O2/O3]  composites. 
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Figure  5.3:  The  effect  of  stacking  sequence  on  the  impact  velocity  threshold.  Com¬ 
parison  between  the  test  data  and  the  predictions. 


84 


Figure  5.4:  A  photograph  of  a  life-size  [02/90/02/90/0/90/0],  specimen  subjected  to 
line-loading  impact. 


broke  into  two  pieces  when  impacted.  Surface  damage  due  to  fiber  breakage  is  believed 
to  have  been  the  initial  damage  that  triggered  the  total  failure.  The  impact  velocity 
threshold  of  the  specimen  was  16  m/sec  at  an  impactor  mass  of  4.17  kg/m,  which  was 
nearly  the  same  as  that  of  the  [O3/9O2/O3/9O2/O3/9O2/O3]  composites.  Apparently,  if 
the  thickness  of  the  surface  ply  group  of  composites  is  too  thin  (i.e.,  less  than  three 
plies),  surface  damage  could  occur  earlier  than  internal  damage,  thus  resulting  in 
premature  failure.  Since  the  surface  damage  caused  by  fiber  breakage  was  not  the 
focus  of  the  study,  the  results  of  the  [02/90/02/90/0/90/0],  composites  were  not  used 
further  in  comparisons  with  others  which  failed  in  totally  different  modes. 

In  order  to  fundamentally  understand  the  phenomena  of  the  test  results,  numeri¬ 
cal  simulations  of  the  test  conditions  were  performed  based  on  the  model  (Chapter  4). 
The  predicted  impact  velocity  thresholds  for  the  tested  composites  with  the  three  dif¬ 
ferent  stacking  sequences  are  presented  in  Figure  5.3.  As  can  be  seen,  the  predictions 
correlated  with  the  data  very  well,  especially  for  the  lower  masses.  The  discrepancy 
between  the  predictions  and  the  data  increased  as  the  mass  of  the  impactor  increased. 


85 


This  can  be  primarily  attributed  to  the  basic  assumption  of  the  analysis  which  was 
based  on  a  small  deformation  theory.  It  was  observed  during  the  experiments  that 
the  amount  of  the  out-of-plane  deflection  of  the  specimens  increased  substantially  as 
the  mass  of  the  impactor  increased.  However,  the  predicted  velocities  overall  were 
conservative  and  consistently  lower  than  the  test  data  by  within  25  percent. 

The  predicted  locations  of  the  initial  critical  matrix  cracks  for  these  specimens  are 
shown  in  Figure  5.5  as  compared  to  the  test  data.  A  photograph  of  a  typical  life-size 
tested  specimen  corresponding  to  each  stacking  sequence  appears  in  Figure  5.6.  It 
is  worth  noting  that  the  locations  of  the  critical  matrix  cracks  shift  from  the  mid¬ 
plane  in  [06/90e/06]  composites  toward  the  back  surface  of  the  specimens  (away  from 
the  impacting  surface)  in  [O3/9O2/O3/9O2/O3/9O2/O3]  composites,  as  the  six  90  degree 
plies  located  at  the  mid-plane  of  the  laminates  are  dispersed  throughout  the  lami¬ 
nates.  For  [0e/906/06]  composites,  the  critical  cracks,  inclined  at  an  angle  of  nearly 
45  degrees,  are  located  in  the  90  degree  central  plies.  However,  for  [04/904/02 /904 /04] 
and  [O3/9O2/ O3/9O2 / O3/9O2 / O3]  composites,  the  critical  cracks  are  found  in  the  90  de¬ 
gree  ply  group  near  the  back  side  of  the  specimens,  and  the  direction  of  the  surface  of 
the  cracks  lean  more  and  more  toward  the  direction  vertical  to  the  loading  direction 
(about  90  degrees  measured  from  the  center  line). 

It  has  been  demonstrated  in  Chapter  4  that  both  interlaminar  shear  stress  cryz  and 
transverse  tensile  stress  crvv  control  the  initial  impact  damage  mechanism  due  to  line¬ 
loading  impact.  Excessive  interlaminar  shear  stress  can  cause  inclined  cracks,  and 
the  high  in-plane  transverse  tensile  stress  can  result  in  vertical  cracks.  For  transverse 
impact,  the  interlaminar  stress  distributions  across  the  thickness  reach  a  maximum 
near  the  mid-plane  of  the  laminate,  but  the  transverse  tensile  stress  is  minimal  near 
the  central  plane  and  increases  toward  the  outer  surfaces  due  to  bending.  According 
to  Eq.  (4.45),  the  combined  ratios  of  the  interlaminar  shear  stress  to  shear  strength 
and  transverse  tensile  stress  to  tensile  strength  determine  the  velocity  threshold  of 
the  laminate  and  the  locations  of  the  critical  cracks. 
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Figure  5.5:  The  effect  of  stacking  sequence  on  the  location  of  the  initial  matrix  crack. 
Comparison  between  the  test  data  and  the  predictions. 
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Figure  5.6:  Photographs  of  three  typical  life-size  tested  specimens  corresponding  to 
three  different  stacking  sequences. 


Figure  5.7  compares  the  calculated  stress/strength  ratio  of  each  stress  compo¬ 
nent  corresponding  to  the  impact  velocity  threshold  for  the  laminates  with  the  three 
different  stacking  sequences.  Note  that  the  sum  of  the  two  ratios  for  each  lam¬ 
inate  is  equal  to  unity.  For  [06/906/Oe]  composites,  initial  failure  is  dominated 
primarily  by  the  interlaminar  shear  stress/strength  ratio.  As  a  result,  the  criti¬ 
cal  cracks  are  inclined  by  nearly  45  degrees  in  the  90  degree  layers.  However,  the 
interlaminar  shear  stress/strength  ratio  is  reduced  for  the  [O4/9O3/O4/9O3/O4]  com¬ 
posites  because  the  90  degree  ply  group,  where  the  cracks  are  generated,  is  away 
from  the  mid-plane.  Hence,  failure  is  due  to  the  combination  of  interlaminar  shear 
stress  and  transverse  tensile  stress.  The  cracks  are  less  inclined  than  those  of  the 
[06/906/06]  composites.  As  the  90  degree  plies  are  further  dispersed  throughout 
the  laminates  in  [O3/9O2/O3/9O2/O3/9O2/O3]  composites,  the  locations  of  the  critical 
cracks  are  closer  to  the  back  of  the  outer  surface.  Accordingly,  the  transverse  tensile 
stress/strength  ratio  further  increases  and  become  a  dominant  factor  for  the  damage 
in  [O3/9O2/O3/9O2/O3/9O2/C3]  composites.  Based  on  the  calculations,  it  is  apparent 
that  both  interlaminar  shear  stress  and  transverse  tensile  stress  are  critically  impor¬ 
tant  for  causing  initial  impact  damage  in  laminated  composites.  By  reordering  the 
stacking  sequence,  the  effect  of  the  combined  stresses  on  the  location  of  the  initial 
failure  can  be  changed  significantly.  A  higher  impact  resistance  can  be  expected  for 
laminates  whose  initial  failure  occurs  away  from  the  mid-plane,  because  a  higher  ve¬ 
locity  is  required  to  make  up  the  reduction  of  the  interlaminar  shear  stress  near  the 
outer  surfaces. 


5.3.2  Thickness  Effect 

The  effect  of  laminate  thickness  on  the  initial  impact  damage  is  demonstrated  by  the 
results  of  the  test  data  on  [(03/903)n/03]  composites  (n  =  2  thru  4),  as  shown  in  Figure 
5.8.  As  n  increased,  the  repeating  ply  group  of  (03/903)  duplicate  n  times  within  the 
laminate;  only  the  thickness  is  changed,  and  the  ply  orientation  and  stacking  sequence 
remain  the  same.  Again,  each  box  in  the  figure  presents  the  data  obtained  from  a 
particular  mass  of  the  impactor.  By  comparing  the  test  data  for  various  laminate 
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Figure  5.7:  The  calculated  on-set  interlaminar  shear  stress/strength  and  transverse 
tensile  stress/strength  ratios  corresponding  to  the  impact  velocity  thresholds  for  lam¬ 
inates  with  three  different  stacking  sequences. 
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thicknesses,  it  is  clear  that  the  change  of  laminate  thickness  does  not  significantly 
alter  the  impact  velocity  threshold.  The  effect  of  thickness  is  not  as  clear  and  obvious 
as  the  stacking  sequence  on  the  initial  impact  damage  of  composites.  For  lower  masses 
(M  =  2.5  kg/m  and  4.17  kg/m),  the  impact  velocity  threshold  varies  within  ±10 
percent  as  the  laminate  thickness  increases  from  15  layers  to  27  layers.  A  consistently 
moderate  increase  of  the  impact  velocity  threshold  is  found  for  the  laminates  impacted 
under  the  heaviest  mass  (A1  =  6.67  kg/m). 

The  predictions  of  the  impact  velocity  thresholds  for  the  laminates  with  various 
thicknesses  based  on  the  analytical  model  are  also  presented  in  Figure  5.8.  Again, 
the  predicted  impact  velocity  thresholds  are  conservative  and  consistently  lower  than 
the  test  data.  The  difference  between  the  predictions  and  the  test  data  becomes 
significant  especially  for  the  thin  laminates  [O3/9O3/O3/9O3/O3]  with  15  layers  and 
for  the  laminates  impacted  under  the  heaviest  mass  (M  =  6.67  kg/m).  This  also 
can  be  primarily  attributed  to  the  large  out-of-plane  deformations  occurring  during 
impact  which  can  not  be  adequately  analyzed  by  the  model  based  on  the  present 
small  deformation  theory.  It  is  believed  that  a  better  prediction  can  be  obtained  by 
adopting  the  large  deformation  theory  in  the  model. 

The  predictions  indicate  that  the  impact  velocity  threshold  is  expected  to  increase 
slightly  as  the  thickness  of  the  composites  increases.  However,  the  rate  of  increase  of 
the  impact  velocity  threshold  with  thickness  is  considerably  smaller  than  that  with 
stacking  sequence.  Accordingly,  the  initial  impact  damage  and  the  subsequent  growth 
of  damage  in  composites  are  much  more  sensitive  to  the  change  of  stacking  sequence 
than  to  the  change  in  thickness. 

The  predicted  locations  of  the  critical  matrix  cracks  for  these  laminates  are  pre¬ 
sented  in  Figure  5.9  along  with  the  measured  test  data.  A  photograph  of  a  life- 
size  specimen  corresponding  to  each  configuration  is  shown  in  Figure  5.10.  Appar¬ 
ently,  the  predictions  correlate  with  the  test  data  very  well.  Figure  5.11  shows  the 
stress/strength  ratio  of  interlaminar  shear  and  transverse  tensile  stresses  in  the  90 
degree  layers  where  the  critical  matrix  cracks  are  predicted  for  the  laminates  with 
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Figure  5.8:  The  effect  of  laminate  thickness  on  impact  velocity  threshold.  Comparison 
between  the  test  data  and  the  predictions. 
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three  different  thicknesses.  It  seems  that  interlaminar  shear  stress  is  more  important 
than  the  transverse  tensile  stress  on  producing  the  initial  damage  for  thin  laminates, 
but  the  trend  is  gradually  reversed  when  the  thickness  increases.  For  [(03/903)4/03] 
composites,  the  contribution  of  the  transverse  tensile  stress  to  damage  is  greater  than 
that  of  the  interlaminar  shear  stress.  However,  the  difference  between  the  two  ratios 
for  each  of  these  three  composites  is  relatively  small  compared  to  that  for  the  compos¬ 
ites  with  various  stacking  sequences.  Accordingly,  this  is  the  reason  why  the  impact 
damage  is  much  less  sensitive  to  laminate  thickness  than  to  stacking  sequence.  It 
is  expected  that,  if  the  initial  impact  damage  is  dominated  by  the  transverse  tensile 
stress,  increasing  the  thickness  could  increase  the  impact  resistance  of  the  compos¬ 
ites.  However,  if  the  failure  mode  of  the  composites  is  strong,)'  associated  with  the 
interlaminar  shear  stress,  then  the  effect  of  thickness  on  impact  damage  is  expected 
to  be  at  a  minimum. 


5.3.3  Mass  Effect 

The  effect  of  the  mass  of  the  impactor  on  the  initial  impact  damage  can  be  best 
demonstrated  in  Figures  5.12  and  5.13  by  reorganizing  the  test  data.  The  ordinate  of 
each  box  in  the  figures  is  the  impact  velocity  threshold,  and  the  horizontal  axis  is  the 
impactor’s  mass.  Both  test  data  and  the  predictions  are  shown  in  the  figures.  Clearly, 
the  impact  velocity  threshold  strongly  depends  upon  the  impactor’s  mass;  the  heavier 
the  mass  is,  the  lower  impact  velocity  is  required  to  cause  damage.  Thus,  the  impact 
velocity  threshold  increases  if  the  mass  of  the  impactor  decreases.  However,  it  is  found 
from  all  the  composites  tested  that  the  percentage  of  increase  of  the  impact  velocity 
threshold  to  the  percentage  of  decrease  of  the  impactor’s  mass  is  out  of  proportion;  a 
nonlinear  relationship  exists  between  the  impact  velocity  threshold  and  the  impactor’s 
mass  for  a  given  laminate  configuration.  Overall,  an  increase  of  the  mass  by  three 
times  approximately  reduces  the  impact  velocity  threshold  by  as  much  as  50  percent. 

The  finding  of  the  test  results  from  Figures  5.12  and  5.13  is  very  important  and 
significant.  It  indicates  that  the  mass  of  the  impactor  should  be  considered  as  a  factor 
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Figure  5.9:  The  effect  of  laminate  thickness  on  the  location  of  the  critical  matrix- 
crack.  Comparison  between  the  test  data  and  the  predictions. 
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Figure  5.10:  Photographs  of  three  typical  life-size  tested  specimens  corresponding  to 
three  different  laminate  thicknesses.  _ _ 
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Stress/Stmgth  Ratio 


Figure  5.11:  The  calculated  on-set  interlaminar  shear  stress/strength  and  transverse 
tensile  stress/strength  ratios  corresponding  to  impact  velocity  threshold  for  laminates 
with  three  different  thicknesses. 
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Figure  5.12:  The  effect  of  impactor’s  mass  on  the  impact  velocity  threshold.  Com¬ 
parison  between  the  test  data  and  the  predictions. 
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Impact  Velocity  Threshold  (m/sec) 


Impactor  mass  (kg/m) 


Figure  5.13:  The  effect  of  impactor ’s  mass  on  the  impact  velocity  threshold.  Com¬ 
parison  between  the  test  data  and  the  predictions. 
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in  determining  impact  damage  in  laminated  composites.  According  to  the  figures, 
impact  damage  is  much  more  sensitive  to  the  impactor’s  velocity  than  to  its  mass. 
Given  the  same  impact  energy,  an  impactor  may  or  may  not  cause  damage,  strongly 
depending  upon  the  mass  of  the  impactor.  Therefore,  it  would  be  inadequate  to 
select  impact  energy  alone  as  an  parameter  to  characterize  impact  damage,  without 
specifying  the  mass  of  the  impactor.  Accordingly,  both  the  impactor’s  mass  and 
velocity  should  be  identified  in  order  to  characterize  impact  damage  in  composites. 

5.4  CONCLUSION 

An  investigation  was  performed  to  study  the  effect  of  laminate  configuration  and 
impactor’s  mass  on  impact  damage  in  laminated  composites.  Both  experimental  and 
analytical  work  were  conducted.  Based  on  the  study,  the  following  remarks  can  be 
made: 

1.  ply  orientation  and  stacking  sequence  can  significantly  affect  impact  damage. 

2.  laminates  with  uniformly  dispersed  ply  orientation  can  increase  impact  resis¬ 
tance  significantly. 

3.  impact  damage  is  more  sensitive  to  the  change  of  stacking  sequence  than  of 
thickness. 

4.  the  mass  of  the  impactor  considerably  affects  the  impact  velocity  threshold  as 
well  as  impact  damage. 

5.  interlaminar  shear  stress/strength  and  in-plane  tensile  stress/ strength  ratios 
dominate  the  initial  impact  damage  and  the  impact  velocity  threshold. 

6.  both  the  velocity  and  mass  of  the  impactor  are  required  in  order  to  characterize 
impact  damage. 

7.  the  use  of  impact  energy  alone  as  a  parameter  to  characterize  impact  damage 
is  inadequate. 
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Chapter  6 


3-DIMENSIONAL  POINT-LOADING  IMPACT  ANALYSIS  AND 

EXPERIMENTS 


6.1  INTRODUCTION 

A  model  was  developed,  based  on  the  previous  study  of  the  line-loading  impact, 
for  predicting  the  impact  damage  of  Graphite/Epoxy  laminated  composites  resulting 
from  point-nose  impact.  In  order  to  verify  the  model,  appropriate  experiments  were 
performed  during  the  investigation.  In  the  following  sections,  the  model  and  experi¬ 
ments  are  described,  and  the  comparisons  between  the  numerical  simulations  and  the 
test  data  from  the  experiments  are  presented. 


6.2  STATEMENT  OF  PROBLEM 

Consider  a  laminated  composite  panel  made  of  Graphite/Epoxy  fiber  reinforced  prepregs 
and  subjected  to  transverse  impact  by  a  low- velocity  spherical  nose  projectile  as  shown 
in  Figure  2.1  The  ply  orientation  of  the  laminate  can  be  arbitrary  but  must  be  sym¬ 
metric  with  respect  to  its  middle  plane.  For  a  given  mass  of  the  impactor,  it  was 
desired  to  determine  the  following: 

1.  The  velocity  of  the  impactor  required  to  initiate  the  impact  damage. 
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2.  The  extent  of  delaminations  inside  the  laminate. 


3.  The  effect  of  ply  orientation  and  laminate  thickness  on  the  impact  damage. 


6.3  ANALYTICAL  MODEL 

A  model  is  proposed  for  predicting  the  damage  in  fiber-reinforced,  epoxy  matrix  lam¬ 
inated  composites  resulting  from  point-loading  impact.  The  model  consists  of  a  stress 
analysis  for  determining  the  stress  distributions  inside  the  laminates  during  impact 
and  a  failure  analysis  for  predicting  the  initiation  and  the  extent  of  the  impact  dam¬ 
age.  The  information  of  the  transient  dynamic  stress  distributions  is  calculated  by 
using  the  three-dimensional  transient  dynamic  finite  element  code  previously  devel¬ 
oped  by  Wu  et  al  [13,16].  In  the  failure  analysis,  a  matrix  failure  criterion  and  a 
delamination  growth  criterion  are  proposed  for  predicting  the  initial  impact  damage 
and  the  extent  of  the  delaminations  due  to  impact,  respectively. 

6.3.1  Stress  Analysis 

The  finite  element  analysis  previously  developed  by  Wu  et  al  [13,16]  is  adopted  for 
calculating  the  stresses  and  strains  inside  the  composites  during  impact  resulting  from 
a  spherical  nose  impactor.  The  information  regarding  the  finite  element  procedures  is 
given  extensively  in  [13,16],  hence,  only  a  brief  description  of  the  analytical  approach 
is  given  as  follows: 

The  analysis  is  based  on  a  three-dimensional  linear  elasticity  theory.  The  ma¬ 
terials  in  each  layer  are  considered  homogeneous  and  orthotropic.  Accordingly,  the 
equilibrium  equations  at  instant  time  t  in  a  variational  form  can  be  expressed  as  [70] 

0=1  w,pui'ttdv+  /  e,jEtJkiikidv  —  /  w,(r,}Tijda  (6.1) 

JV  J  l» 
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where  are  the  stresses,  e*/  are  the  strains,  p  is  the  density,  uftt  are  the  accelerations 
(ui.tt  =  d2Uj/dt2),  Wi  are  the  arbitrary  variational  displacements,  e,7  are  the  strains 
from  the  arbitrary  variational  displacements,  v  is  the  entire  plate  volume,  sa  is  the 
surface  of  the  plate,  rij  is  the  outward  unit  normal  vector  on  the  plate  surface,  and 
Eijki  are  the  material  properties  of  the  laminate,  which  may  vary  from  layer  to  layer 
according  to  the  ply  orientation  of  the  composite. 

In  order  to  solve  Eq.  (6.1),  the  distribution  of  the  contact  force,  F(=  <7,jn;), 
between  the  impactor  and  the  impacted  laminate  must  first  be  known.  The  projectile 
is  modeled  as  an  elastic  body  with  a  spherical  nose.  The  contact  force  distribution 
during  impact  is  simulated  according  to  loading  and  unloading  processes. 

Upon  loading  (the  contact  force  was  increased),  the  contact  force  distribution  is 
determined  using  the  Hertzian  contact  law  [71].  Thus,  the  contact  force  F  can  be 
related  to  the  indentation  depth  a  (the  distance  between  the  center  of  the  projectile’s 
nose  and  the  mid-surface  of  the  plate)  by  the  expression  [71] 


F  =  kq1'5  (6.2) 

where  k  is  the  modified  constant  of  the  Hertz  contact  theory  proposed  by  Sun  et  al. 
[22]  and 


*  3^  [(1  -  +  l/E„)  (6'3) 

where  r,  v,,  and  E,  are  the  local  radius,  the  Poisson’s  ratio,  and  the  Young’s  modulus 
of  the  impactor,  respectively.  Eyy  is  the  transverse  modulus  normal  to  the  fiber 
direction  in  the  upper-most  composite  layer. 

Upon  unloading,  the  contact  force  is  simulated  by  the  following  relation  developed 
by  Sun  et  al.  [22] 
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2.5 


(6.4) 


F  = 


Fm 


a  -  Q0 

Otm  OCq 


where  Fm  is  the  maximum  contact  force  just  before  unloading,  am  is  the  maximum 
indentation  corresponding  to  Fm,  and  oc0  is  the  permanent  indentation  during  the 
loading/unloading  process.  The  permanent  indentation  can  be  determined  from  the 
following  expression  [22] 


<*o 

at0 


when  am  <  acr 
when  Qm  >  acr 


(6.5) 


where  acr  is  the  critical  indentation,  and  is  approximately  0.004  inches  for  glass/epoxy 
and  0.00316  inches  for  graphite/epoxy. 

An  eight-node  brick  element  incorporating  incompatible  modes  is  used  in  the  finite 
element  calculations,  and  a  direct  Gauss  quadrature  integration  scheme  is  adopted 
[70]  through  the  element  thickness  to  account  for  the  change  in  material  properties 
from  layer  to  layer  within  the  element.  Therefore,  plies  with  different  ply  orientations 
can  be  grouped  into  an  element,  resulting  in  a  significant  reduction  in  computational 
time  and  memory  space  for  the  three-dimensional  analysis. 

The  accuracy  of  the  computer  code  is  extensively  verified  by  comparison  with  the 
existing  analytical  and  numerical  solutions  [13,16].  To  demonstrate  the  capability  of 
the  finite  element  analysis,  Figure  6.1  shows  the  calculated  impact  force  and  central 
laminate  deflection  as  a  function  of  time  for  a  T300/976  [454/  —  45g/454]  composite 
plate.  The  finite  element  mesh  used  for  the  calculations  is  also  presented  in  Figure 
6.1.  The  laminate  is  discretized  into  a  total  of  N\  x  N2  x  N3  elements.  A  total  of 
four  elements  ( N3  =  4)  are  used  through  the  thickness  of  the  laminate.  Because 
of  low-velocity  impact,  the  maximum  impact  force  occurred  at  time  t  =  620/isec, 
corresponding  well  to  the  time  at  which  the  central  deflection  of  the  plate  reaches  the 
maximum.  After  the  rebound  of  the  impactor,  the  plate  vibrates  harmonically. 
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Finite  Element  Mesh 


N3 


Figure  6.1:  The  finite  element  calculations  of  the  impact  force  and  the  central  deflec¬ 
tion  of  a  [45«/  -  458/45«]  composite  plate  subjected  to  transverse  impact. 
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6.3.2  Failure  Analysis 


Once  the  stresses  are  obtained  from  the  finite  element  analysis,  two  failure  criteria  are 
proposed  to  predict  the  initiation  of  the  impact  damage  and  the  extent  of  the  delam¬ 
inations  in  the  composites  resulting  from  the  impact.  Since  impact  damage  is  a  very 
complicated  phenomenon,  predicting  such  damage  requires  a  thorough  understanding 
of  the  basic  damage  mechanisms  and  mechanics  governing  the  impact  damage  event. 
Based  on  the  previous  line-loading  impact  study  (Chapter  3,  4  and  5),  the  following 
are  concluded  for  line-loading  impact  damage: 

1.  Intraply  matrix  cracking  (referred  to  as  the  critical  matrix  cracks)  is  the  initial 
impact  damage  mode. 

2.  Delaminations  initiates  from  the  critical  matrix  cracks  which  propagates  into 
the  nearby  interface  with  the  dissimilar  materials. 

3.  If  a  critical  crack  is  located  within  the  inner  plies  of  the  laminates,  the  delami¬ 
nation  along  the  bottom  interface  of  the  cracked  ply  propagates  away  from  the 
location  of  the  impact  (see  Figure  6.2). 

4.  If  a  critical  crack  is  located  at  the  surface  ply  of  the  laminate,  a  delamination 
propagates  from  the  critical  crack  away  from  the  center  of  the  impact  along  the 
first  interface  of  the  crack  ply  (see  Figure  6.2). 

Although  the  results  of  the  study  are  obtained  from  the  line-loading  impact,  it 
is  believed  that  similar  phenomena  occur  in  point-loading  impact.  Therefore,  it  is 
postulated  that  the  damage  mechanisms  of  laminated  composites  resulting  from  point¬ 
loading  impact  follow  the  same  sequences  as  are  found  from  line-loading  impact  as 
follows: 

Transverse  impact  first  initiates  matrix  cracks  in  a  layer  within  the  laminate. 
Immediately,  delaminations  are  produced  along  the  bottom  or  upper  interface  of  the 
cracked  layer,  depending  on  the  position  of  the  cracked  layer  in  the  laminate.  As 
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the  duration  of  impact  proceeds,  additional  matrix  cracks  occur  subsequently  in  the 
other  layers  and  produce  additional  delaminations  along  the  other  interfaces. 

Critical  Matrix  Cracking  Criterion 

In  order  to  predict  the  occurrence  of  the  matrix  cracking,  the  matrix  failure  cri¬ 
terion  proposed  previously  for  the  line-loading  impact  is  adopted  here;  the  criterion 
can  be  expressed  as 


(n°yyV  (nTT yA2  2  f  eM  >  1  Failure 

\  nYt  )  \  nSi  )  M  \  eM  <  1  No  failure 

which  is  the  same  as  Eq.  (4.45)  in  Chapter  4. 


(6.6) 


Whenever  the  calculated  averaged  stresses  in  any  one  of  the  plies  in  the  laminate 
first  satisfy  the  criterion  (e\f  =  1)  during  impact,  initial  impact  damage  is  predicted. 
The  time  t  corresponding  to  the  initial  damage  is  designated  as  tM.  A  delamination 
could  immediately  follow  from  the  location  of  the  matrix  crack.  As  the  time  increases 
( t  >  <Af)  during  impact,  additional  matrix  cracking  could  be  produced  in  the  other 
layers.  Hence,  the  criterion  should  continuously  be  applied  at  the  other  layers  for 
determining  any  additional  matrix  failure.  It  is  noted  that,  if  no  additional  matrix 
cracking  is  found  at  any  other  layers  during  impact,  then  the  impactor’s  velocity  asso¬ 
ciated  with  the  only  matrix  cracking  is  defined  here  as  the  impact  velocity  threshold 
which  is  the  velocity  required  to  just  cause  the  initial  impact  damage  of  the  laminate. 

Delamination  Growth  Criterion 


Once  a  critical  matrix  crack  is  predicted  in  a  ply  within  the  laminate,  a  delamina¬ 
tion  can  be  initiated  from  the  crack.  To  accurately  simulate  delamination  propagation 
is  very  difficult  and  complicated,  involving  multiple  dynamic  crack  propagation  and 
delamination  surfaces  interaction.  Hence,  it  is  a  formidable  task  to  undertake.  There¬ 
fore,  in  the  investigation,  a  semi-empirical  model  is  proposed  to  estimate  the  extent 
of  delaminations  in  the  composites  after  impact.  Basically,  there  are  two  types  of 
critical  cracks  initiating  delamination  resulting  from  impact:  one  can  be  referred  to 


107 


as  the  shear  crack  generated  within  the  laminates  and  the  other  can  be  referred  to 
as  the  bending  crack  produced  from  the  bottom  surface  of  the  laminates  as  shown  in 
Figure  6.2.  In  order  to  effectively  estimate  the  extent  of  the  delamination  growth,  it 
is  necessary  to  include  in  the  model  the  basic  governing  parameters  controlling  the 
delamination  propagation,  once  it  is  initiated. 

It  is  well  reported  [14,30,31,33,46,53]  that  the  delamination  resulting  from  point- 
nose  loading  quasi-statically  or  dynamically  appears  generally  in  a  peanut  shape,  as 
shown  in  Figure  6.3.  The  longitudinal  axis  of  the  delamination  tends  to  orient  itself 
in  the  direction  parallel  to  the  fiber  direction  of  the  bottom  layer  below  the  interface. 
It  has  also  been  shown  that,  under  a  quasi-static  loading  condition,  the  growth  of 
a  delamination  induced  by  a  shear  crack  is  controlled  by  both  Mode  I  and  Mode  II 
fractures  due  to  excessive  interlaminar  shear  stresses  such  as  reported  in  [43,72.73]. 
The  growth  of  the  delamination  along  the  fiber  direction  of  the  bottom  layer  at  the 
interface  is  very  unstable  once  it  begins  to  propagate  [73].  From  the  previous  two- 
dimensional  line-loading  impact  study  (Chapters  3  and  4),  it  is  observed  empirically 
that  the  shear  cracks  and  interface  delamination  always  appear  simultaneously  and 
delamination  always  extends  fully  to  the  boundary  from  the  matrix  crack  tips  as 
shown  in  Figure  6.4. 

However,  several  investigators  [43,73,74]  have  also  demonstrated  recently  that  de¬ 
lamination  growth  induced  by  a  bending  crack  is  governed  primarily  by  Mode  I  tensile 
fracture  due  to  in-plane  bending  stress  (transverse  stress).  The  growth  of  the  delami¬ 
nation  induced  by  a  bending  crack  is  quite  stable  in  the  direction  normal  to  the  fiber 
direction  of  the  bottom  layer  beneath  the  interface  [73].  The  previous  two-dimensional 
study  also  reveals  a  stable  delamination  growth  due  to  a  bending  crack,  as  shown  in 
Figure  6.5.  Hence,  once  a  delamination  is  initiated  from  a  critical  matrix  crack,  it 
can  grow  much  more  extensively  along  the  fiber  direction  than  in  the  transverse  di¬ 
rection  of  the  bottom  layer  at  the  interface.  This  may  provide  an  explanation  why 
delamination  appears  to  be  in  a  peanut  shape  in  laminated  composites. 
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Figure  6.3:  A  typical  delamination  shape  in  laminated  composites  induced  by  a  point 
nose  impact. 
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Impact  Velocity  (m/s) 


Figure  6.4:  Measured  delamination  area  within  the  [Oe/906/06]  laminated  composite 
plates  resulting  from  2-dimensional  line  loading  impact  at  different  impact  velocities. 
Delamination  area  was  normalized  with  respect  the  corresponding  plate  area. 


Impact  Velocity  (m/s) 

Figure  6.5:  Measured  delamination  area  within  the  [904/0iO/904]  laminated  composite 
plates  resulting  from  2-dimensional  line  loading  impact  at  different  impact  velocities. 
Delamination  area  was  normalized  with  respect  to  the  corresponding  plate  area. 
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Therefore,  it  is  believed  that,  for  a  shear-crack  induced  delamination,  the  inter¬ 
laminar  longitudinal  shear  stress  (along  the  fiber  direction)  axz  in  the  layers  right 
below  the  interface  governs  the  major  delamination  growth,  while  the  interlaminar 
transverse  shear  stress  oyz  in  the  layers  right  above  the  interface  could  also  enhance 
the  delamination  propagation  depending  on  the  direction  of  its  ply  orientation.  A 
schematic  description  of  the  delamination  growth  mechanism  induced  by  a  shear 
crack  is  shown  in  Figure  6.6.  However,  for  a  bending  crack-induced  delamination, 
the  interlaminar  longitudinal  shear  stress  crxz  in  the  layers  right  below  the  interface 
still  controls  the  major  delamination  growth,  but  the  in-plane  bending  stress  eryy  in 
the  layers  right  below  the  interface  advances  the  delamination  propagation  in  the 
secondary  direction  (normal  to  the  fiber  direction).  A  schematic  description  of  the 
delamination  growth  mechanism  is  shown  in  Figure  6.7. 

Accordingly,  by  taking  both  failure  mechanisms  into  consideration,  it  is  considered 
that  the  distributions  of  the  interlaminar  longitudinal  shear  stress  axz  and  transverse 
in-plane  stress  cryy  throughout  the  thicknesses  of  the  bottom  layer  beneath  the  inter¬ 
lace  and  the  interlaminar  transverse  shear  stress  ayz  in  the  upper  layer,  contribute 
primarily  to  the  delamination  growth  resulting  from  point-nose  impact.  Therefore,  it 
is  postulated  that  delamination  growth  due  to  low-velocity  impact  occurs  only  when 
the  following  two  sequential  conditions  are  met: 

1.  one  of  the  ply  groups  intimately  above  or  below  the  concerned  interface  has 
failed  due  to  matrix  cracking  and 

2.  the  combined  stresses  governing  the  delamination  growth  mechanisms  through 
the  thicknesses  of  the  upper  and  lower  ply  groups  of  the  interface  reach  a  critical 
value. 

Based  on  the  hypothesis,  the  following  delamination  growth  criterion  for  low- 
velocity  impact  is  proposed 
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Impact  Damage  Growth  Mechanism 


Figure  6.6:  A  schematic  description  of  the  delamination  growth  mechanism  induced 
by  a  shear  crack  in  a  laminated  composite  subjected  to  point  nose  impact. 
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Impact  Damage  Growth  Mechanism 


Delamination  induced  by  Bending  Cracking 


Figure  6.7:  A  schematic  description  of  the  delamination  growth  mechanism  induced 
by  a  bending  crack  in  a  laminated  composite  subjected  to  point  nose  impact. 
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1  e£>  <  1 


Failure 
No  failure 


(6.7) 
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D 


where 

n+iy  =„+l  if  n+1eyy  >  0 

"+,y=n+,r£  ifn+1ayy<0 


(6.8) 


Da  is  an  empirical  constant  which  has  to  be  determined  from  the  experiments  de¬ 
scribed  in  the  Experiments  section.  Once  Da  is  chosen,  it  is  found  to  be  quite  insen¬ 
sitive  to  the  ply  orientation  and  thickness  of  the  laminates,  and  primarily  dependent 
only  on  the  material  system  used.  The  subscripts  x,  y ,  and  r  are  the  local  material 
coordinates  of  an  individual  ply  within  the  laminate,  and  the  superscripts  n  and  n  +  1 
correspond  to  the  upper  and  lower  ply  groups  of  the  n  —  th  interface,  respectively. 
ayz  and  ayy  are  the  averaged  interlaminar  and  in-plane  transverse  stresses  within  the 
n  —  th  and  n  +  1  —  th  ply,  respectively,  defined  in  Eq.  (4.46)  and  (4.47).  axz  is  the 
averaged  interlaminar  longitudinal  stresses  within  the  n  -f  1  —  th  ply  which  can  be 
expressed  as 


n+'*tz  = -j-  axzdz  (6.9) 

where  tn  and  tn_i  are  the  upper  and  lower  interfaces  of  the  n  —  th  ply  or  ply  group 
in  the  laminate  and  hn  is  the  thickness  of  the  ply  or  the  ply  group. 

Accordingly,  once  a  critical  matrix  crack  is  predicted  in  a  layer,  the  delamination 
growth  criterion  is  then  applied  to  estimate  the  extent  of  the  delamination  along  the 
interface  of  the  cracked  ply  in  the  laminate.  The  procedure  for  determining  the  extent 
of  the  impact  damage  can  be  described  as  follows: 

1.  Calculating  transient  dynamic  stresses  within  each  layer  as  a  function  of  time. 

2.  Applying  the  matrix  failure  criterion  for  predicting  the  critical  matrix  cracks  in 
each  layer. 
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3.  If  matrix  cracking  is  predicted  in  a  layer  of  the  laminate,  applying  the  delam¬ 
ination  criterion  subsequently  in  the  upper  and  bottom  layer  of  the  interface 
during  the  entire  period  of  impact. 


The  above  procedure  is  repeated  at  the  other  layers  during  impact  for  determining 
any  additional  matrix  cracking  and  delaminations.  The  final  size  of  each  delamination 
is  determined  by  the  area  within  which  the  stress  components  satisfy  the  delamination 
growth  criterion  during  the  entire  duration  of  impact.  No  material  degradation  was 
considered  in  the  model,  and  it  is  also  noted  that  the  model  does  not  take  into  account 
the  delamination  interaction  during  impact  which  may  be  important  for  multiple 
delaminations. 


6.4  EXPERIMENTS 

In  order  to  verify  the  proposed  model,  numerical  predictions  from  the  model  were 
compared  with  the  data  generated  by  the  tests  performed  during  the  investigation. 
The  same  impact  testing  facility  described  in  Chapter  3  was  used  for  the  experiments. 
A  spherical-nosed  impactor  was  selected  for  the  study  (see  Figure  3.2).  The  radius 
of  the  spherical  nose  head  made  of  steel  was  0.635  cm.  The  specimens  were  firmly 
clamped  along  two  parallel  edges  as  shown  in  Figure  3.1. 

T300/976  composites  were  selected  for  the  study.  Numerous  different  ply  orien¬ 
tations  were  chosen.  The  dimensions  of  each  specimen  were  10  cm  long  and  6.6  cm 
wide.  Table  6.1  lists  the  ply  orientations  and  the  configurations  of  the  specimens 
used  in  the  tests.  All  the  specimens  were  cured  under  the  cure  cycle  previously  men¬ 
tioned  in  Chapter  3  which  was  found  to  produce  finished  composites  of  fairly  good 
quality,  without  thermal-induced  pre-matrix  cracks.  All  the  specimens  were  cut  by 
a  diamond  coated  saw  and  X-rayed  after  the  cutting  to  inspect  any  internal  damage 
due  to  manufacturing  and  cutting.  No  apparent  damage  was  found  after  cutting  the 
specimens. 
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Ply  Orientation 

Thickness  (h) 
(mm) 

Span  length  (L) 
(cm) 

Width  (W) 
(cm) 

[454/-45S/454] 

2.30 

10.0 

7.6 

[03/903/03/903/03] 

2.16 

10.0 

7.6 

[04/452/-454/452/04] 

2.30 

10.0 

7.6 

[04/454/-454/904/-454/454/04] 

4.03 

10.0 

7.6 

Clamped  Area 
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Fifteen  replicants  were  prepared  for  the  three  ply  orientations:  [454/  —  458/454], 
[O3/9O3/O3/9O3/O3],  and  [04/452/  —  454/452/04].  Each  specimen  was  impacted  only 
once  at  a  selected  impact  velocity.  Specimens  were  then  X-rayed  using  enhanced  dye- 
penetrate.  For  some  specimens  without  any  visual  surface  damage,  dye  might  not 
have  penetrated  effectively  into  the  specimens  to  produce  maximum  effect.  In  these 
cases,  a  small  hole  was  drilled  at  the  impacted  location  and  dye- penetrate  was  applied 
through  the  hole.  An  X-Radiograph  was  then  taken.  For  a  given  configuration,  the 
velocity  of  the  impactor  was  started  at  a  very  low  level  and  gradually  increased  until 
a  considerable  amount  of  damage  was  visualized  from  X- Radiographs. 

The  velocity  just  causing  the  initial  matrix  cracking  was  recorded  and  considered 
as  the  impact  velocity  threshold  for  the  test  specimen  under  the  given  boundary  con¬ 
dition.  A  peanut  shape  as  given  in  Figure  6.4  was  found  to  be  the  typical  shape  for  the 
delaminations  from  the  X-Radiographs.  This  is  also  consistent  with  the  findings  ob¬ 
served  by  other  investigators  [14,30,31,33,46,53].  The  sizes  of  the  delaminations  shown 
in  the  X-radiographs  were  then  measured  and  recorded  in  terms  of  their  length  (lon¬ 
gitudinal  direction)  and  width  (transverse  direction)  as  a  function  of  the  impactor  s 
velocity.  Although  the  delamination  size  would  be  estimated  from  the  X- Radiographs 
indicated  by  the  white  color  area,  the  actual  size  of  the  delamination  was  still  very 
difficult  to  determine,  especially  when  the  delamination  size  was  small.  In  order  to 
confirm  the  direct  measurement  of  the  delamination  size  from  the  X- Radiographs, 
selected  specimens  were  also  sliced  along  the  longitudinal  or  tranverse  axis  of  the  rel¬ 
evant  delamination,  and  the  sliced  cross-section  of  the  specimens  was  then  X-rayed. 
The  length  and  width  were  then  remeasured  to  compare  with  the  previous  measure¬ 
ments.  Overall,  both  types  of  measurements  were  quite  consistent  with  each  other. 
However,  an  error  of  about  ±10  percent  in  the  measurement  from  the  actual  size  could 
possibly  have  been  introduced.  The  measured  delamination  sizes  are  summarized  in 
the  Appendix  C. 
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6.5  SENSITIVITY  STUDY 


Determining  the  value  of  Da  in  the  delamination  growth  criterion  requires  the  use  of 
impact  test  results.  The  measured  delamination  sizes  and  shapes  from  X-radiographs 
were  used  for  determining  the  value.  The  test  results  of  [O3/9O3/O3/9O3/O3]  composites 
were  first  selected  and  compared  with  numerical  simulations  based  on  the  model.  A 
value  of  l.S  was  chosen  for  Da  because  the  prediction  of  the  delamination  sizes  would 
best  fit  the  measured  test  data. 

In  order  to  determine  the  sensitivity  of  the  empirical  constant  Da  to  lamination 
configuration,  the  delamination  lengths  of  the  laminates  with  three  tested  ply  orien¬ 
tations  were  calculated  from  the  values  of  Da  varying  in  a  wide  range.  The  calculated 
lengths  were  then  compared  with  the  corresponding  data  from  the  experiments.  A 
formula  [75]  was  selected  to  estimate  the  discrepencies  between  the  predictions  and 
measurements  based  on  various  of  Da.  The  error  formula  is  expressed  as 


where 


Lc  =  Calculated  Delamination  Length 
Lt  =  Tested  Delamination  Length 
L  =  Specimen  Length  (10  cm) 

N(  =  Number  of  Test  points 


(6.10) 


(6.11) 


The  comparison  of  the  study  is  presented  in  Figure  6.8.  As  can  be  seen  clearly, 
the  value  of  Da  which  best  matched  to  the  data  varied  from  1.2  to  2.0  for  the  ply 
orientations  studied.  An  average  value  of  1.8  would  be  a  fair  selection  for  T300/976 
composites.  The  effect  of  the  variation  of  Da  values  on  the  prediction  of  the  delami¬ 
nation  size  can  also  be  seen  from  the  Figure  6.8.  Depending  upon  the  location  of  the 
delamination  in  laminates,  Figure  6.8  shows  that  a  change  of  the  parameter  by  up  to 
±50%  from  the  averaged  value  does  not  affect  significantly  the  predicted  delamina¬ 
tion  size.  Therefore,  it  seems  that  the  value  of  Da  is  insensitive  to  the  change  of  ply 
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Figure  6.8:  The  sensitivity  evaluation  of  the  effect  of  Da  value  on  the  predicted 
delamination  size  as  compared  to  the  experiments. 


orientation  and  thickness,  but  may  depend  on  the  material  system  chosen.  However, 
a  study  may  be  needed  to  further  evaluate  the  characteristic  of  the  parameter.  In  this 
investigation,  once  the  value  of  Da  was  determined,  it  was  then  used  as  a  constant  in 
the  model  for  the  T300/976  composites.  The  material  properties  of  T300/976  com¬ 
posites  used  in  the  calculations  are  listed  in  Table  4.1.  The  results  of  the  predictions 
compared  with  experiments  are  presented  in  the  next  section. 

In  order  to  evaluate  the  effect  of  finite  element  mesh  size  on  the  stability  of  the 
numerical  predictions  based  on  the  model,  numerical  calculations  are  obtained  based 
on  different  meshes,  ranging  from  very  fine  to  rather  coarse.  Figure  6.9  shows  the 
comparisons  of  the  calculations  in  which  the  meshes  were  generated  by  proportionally 
decreasing  or  increasing  the  number  of  elements  along  length  (Ari)  and  width  (A2) 
direction  of  the  laminate  while  keeping  the  number  of  total  elements  (A3)  through  the 
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thickness  fixed.  It  can  be  seen  that  the  predicted  delamination  sizes  quickly  converged 
as  the  number  of  the  elements  was  more  than  192  {Ni(=  16)  x  jV2(=  12)),  beyond 
which  the  calculated  maximum  impact  force  and  maximum  deflection  of  the  plate 
were  also  fairly  stabilized.  The  results  indicate  that  the  proposed  failure  analysis 
does  not  require  the  use  of  an  extensively  fine  mesh.  Therefore,  a  regular  mesh  of 
768  (N\{=  16)  x  A2(=  12)  x  N3(=  4))  was  used  in  the  numerical  calculations  for 
generating  the  results  in  comparison  with  the  test  data. 


6.6  RESULTS  AND  DISCUSSIONS 
[454/-  458/454] 

Figure  6.10  shows  three  X-Radiographs  of  [454 /  —  458/454]  composites  in  series 
after  being  hit  ea  ?,  a  different  velocity  by  an  impactor.  The  top  X-radiograph  was 
taken  from  the  'aininate  subjected  to  impact  at  a  velocity  slightly  above  the  impact 
velocity  threshold  of  the  laminate.  Two  distinct,  parallel,  short  matrix  cracks  which 
are  located  near  the  impacted  area  within  the  middle  —45  degree  ply  group  of  the 
laminate  are  clearly  shown  in  Figure  6.10.  A  long  matrix  crack  aligned  in  the  45 
degree  fiber  direction  and  located  in  the  bottom  45  degree  ply  group  is  also  clearly 
indicated.  These  cracks  are  most  likely  associated  with  the  critical  matrix  cracks 
from  which  delamination  is  initiated.  A  white  color  area  in  a  peanut  shape  oriented 
along  the  45  degree  fiber  direction  shows  that  a  delamination  exists  at  the  interface 
between  the  middle  -45  degree  ply  group  and  the  bottom  45  degree  ply  group. 

As  the  impactor’s  velocity  increased,  the  first  delamination  found  at  the  last  in¬ 
terface  between  —45/45  ply  groups  was  substantially  increased,  as  shown  in  the  X- 
Radiograph  located  at  the  middle  of  Figure  6.10.  Meanwhile,  a  relatively  smaller 
delamination  at  the  first  interface  between  45/  -  45  ply  groups  was  also  spotted, 
in  addition  to  a  considerable  number  of  microcracks  in  the  45  and  —45  degree  ply 
groups.  At  a  much  higher  velocity,  as  given  in  the  X- Radiograph  at  the  bottom  of 
Figure  6.10,  the  first  delamination  at  the  last  interface  always  governed  the  major 
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Figure  6.9:  The  sensitivity  evaluation  of  the  effect  of  finite  element  mesh  size  on  the 
calculated  delamination  sizes  based  on  the  proposed  model. 
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T300/976  [454/458/454] 


Mass  =  0. 1 6  kg  Velocity  =  4.15  m/s 


Mass  =  0.16  kg  Velocity  =  5.89  m/s 


Mass  =  0.16  kg  Velocity  =  9.02  m/s 


Figure  6.10:  Three  X-Radiographs  of  [454 /  —  458/454]  composites  in  s 
to  impact  by  a  spherical  nosed  projectile  at  three  different  velocities. 


ies  subjected 


122 


delamination  size,  which  was  considerably  larger  than  the  second  delamination. 

A  summary  of  the  measured  delamination  sizes  in  [45 4/  —  458/454]  composites  as 
a  function  of  the  velocity  of  the  impactor  with  0.16  kg  is  shown  in  Figure  6.11.  Solid 
circular  and  rectangular  symbols  represent  the  measured  delamination  length  and 
width  in  the  longitudinal  and  transverse  directions  as  a  function  of  the  impactor’s 
velocity,  respectively.  The  upper  figure  corresponds  to  the  first  delamination  along  the 
—45/45  degree  plies  interface,  and  the  bottom  one  relates  to  the  second  delamination 
which  occurred  at  the  upper  45/  —  45  degree  plies  interface. 

Clearly,  there  exists  an  impact  velocity  threshold  for  the  laminate.  The  first 
delamination  apparently  appears  much  earlier  than  the  second  delamination.  There 
also  seems  to  be  a  velocity  threshold  for  the  occurrence  of  each  delamination.  It  is 
worth  noting  that  no  matrix  cracks  or  delaminations  were  found  in  the  laminate  in 
the  X-Radiographs  when  the  impactor’s  velocity  was  less  than  3.3  m/s. 

The  predicted  delamination  size  at  each  interface  as  a  function  of  the  impactor's 
velocity  is  presented  by  the  solid  and  dashed  lines  in  Figure  6.1 1.  The  impact  velocity 
threshold  is  slightly  underestimated  but  still  agrees  with  the  data  reasonably  well. 
The  first  critical  matrix  cracking  is  predicted  at  the  middle  —45  degree  layer  of  the 
laminate.  The  predicted  delamination  length  and  width  also  correlate  well  with  the 
data.  Not  only  does  the  prediction  match  with  the  size  of  the  first  delamination,  but 
also  it  correlates  very  well  with  the  initiation  and  the  size  of  the  second  delamination. 

Although  the  test  data  is  quite  scattered,  a  relationship  which  is  quite  consistent 
with  the  predictions  seems  to  exist  between  the  delamination  size  and  the  velocity  of 
the  impactor.  The  longitudinal  length  of  the  delamination  seems  to  be  more  sensitive 
than  the  transverse  width  to  the  increase  in  the  impactor’s  velocity.  The  longitudinal 
length  of  the  delamination  is  always  oriented  along  the  fiber  direction  of  the  bottom 
ply  below  the  delamination  interface. 

As  a  comparison,  the  numerical  simulations  of  the  delamination  size  of  the  com¬ 
posites  subjected  to  impact  at  the  velocities  corresponding  to  those  given  in  Figure 
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Velocity  (m/s) 

Figure  6.11:  Comparison  of  the  delamination  sizes  of  T300/976  [454/  —  45g/454] 
composites  between  the  measurements  from  the  experiments  and  the  calculations 
based  on  the  model. 
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6.10  are  generated,  and  the  results  of  the  predictions  are  plotted  in  Figure  6.12.  The 
delamination  contours  shown  in  Figure  6.12  bound  the  area  of  the  predicted  delami¬ 
nation.  Corresponding  to  each  velocity  given  in  Figure  6.10,  the  predicted  shape  and 
size  of  each  delamination  shown  in  Figure  6.12  are  quite  consistent  with  the  results 
of  the  X-Radiographs  given  in  Figure  6.10. 

[O3/9O3/O3/9O3/O3] 

The  X-Radiographs  of  three  [O3/9O3/O3/9O3/O3]  composites  impacted  by  an  im- 
pactor  at  three  different  velocities  are  shown  in  Figure  6.13.  Again,  the  impact 
damage  shown  in  the  top  X-Radiograph  given  in  the  figure  corresponds  to  the  impact 
velocity  threshold.  Differing  from  [464/  —  45g/ 454]  composites,  only  a  single  matrix 
crack  located  along  the  center  line  of  the  bottom  0  degree  plies  was  found.  A  small 
peanut-shaped  delamination  was  also  found  located  at  the  interface  between  the  last 
90  degree  ply  group  and  the  cracked  outermost  0  degree  ply  group. 

Multiple  delaminations  occured  in  the  laminates  as  shown  in  the  middle  and  bot¬ 
tom  X-Radiographs  of  Figure  6.13,  as  the  impactor’s  velocity  increased.  Again,  each 
delamination  oriented  itself  along  the  fiber  direction  of  the  bottom  ply  of  the  delam¬ 
inated  interface.  Three  delaminations  were  found  in  the  X-Radiographs  which  were 
also  confirmed  by  a  X-Radiograph  taken  from  a  cross-secion  of  a  sliced  specimen.  The 
major  delamination  size  was  still  governed  by  the  first  delamination  which  occurred 
at  the  bottom  interface. 

The  measured  length  and  width  of  each  delamination  in  [O3/9O3/O3/9O3/O3]  com¬ 
posites  as  a  function  of  impactor’s  velocity  are  presented  in  Figure  6.14.  Since  only 
three  delaminations  occurred  in  the  laminates,  no  delamination  was  found  at  the 
first  O3/9O3  interface.  The  first  delamination  which  occurred  at  the  last  interface  ap¬ 
peared  earlier  than  the  others.  The  results  of  the  test  indicate  that  the  size  of  the 
first  delamination  is  much  more  sensitive  to  the  impactor’s  velocity  than  are  the  other 
delaminations. 

The  predicted  delamination  length  and  width  for  each  delamination  as  a  function 


125 


T300/976  [45V-458/454] 


Mass  =  0.16  kg  Velocity  =  5.89  m/s 


Figure  6.12:  The  predicted  delamination  sizes  of  [454/  —  458/454]  composites  corre¬ 
sponding  to  the  impact  condition  given  in  Figure  6.10. 
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T300/976  [03/903/03/903/03] 


<y 


Mass  =  0.16  kg  Velocity  =  3.22  m/s 


Mass  =  0. 1 6  kg  Velocity  =  4.0  m/s 


Mass  =  0.16  kg  Velocity  =  6.70  m/s 


Figure  6.13:  Three  X-Radiographs  of  [O3 / 903 / 03 / 9O3 / O3]  composites  in  series  sub¬ 
jected  to  impact  by  a  spherical  nosed  projectile  at  three  different  velocities. 
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Velocity  (m/s) 

Figure  6.14:  Comparison  of  the  delamination  sizes  of  T300/976  [O3/9O3/O3/9O3/O3] 
composites  between  the  measurements  from  the  experiments  and  the  calculations 
based  on  the  model. 
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of  the  impactor’s  velocity  are  also  indicated  in  Figure  6.14  by  solid  and  dashed  lines. 
It  is  worth  noting  that  the  model  predicts  matrix  cracking  at  the  bottom  0  degree  plies 
as  the  initial  failure  mode,  which  correlates  very  well  with  the  X- Radiograph’s  finding 
,  shown  in  Figure  6.13.  The  predicted  impact  velocity  threshold  is  also  very  close  to 
the  measured  one.  A  delamination  is  predicted  to  be  generated  at  the  last  interface 
between  9O3/O3  ply  groups  immediately  after  the  occurrence  of  matrix  cracking. 

Matrix  cracking  is  also  predicted  afterward  in  the  inner  90  and  0  degree  plies 
which  initiates  second  and  third  delaminations  as  shown  in  Figure  6.14.  Overall, 
the  predicted  delaminations  are  quite  consistent  with  the  measured  ones.  No  matrix 
cracking  is  predicted  in  the  top  surface  0  degree  ply  group,  hence  there  is  no  de¬ 
lamination  predicted  to  occur  at  the  first  interface  between  the  surface  0  degree  ply 
group  and  the  second  90  degree  ply  group.  This  prediction  is  also  confirmed  by  the 
experiment. 

The  experimental  results  shown  by  the  X-Radiographs  given  in  Figure  6.13  are 
numerically  simulated  by  the  model.  The  numerical  simulations  of  the  delamination 
sizes  corresponding  to  the  test  conditions  are  shown  in  Figure  6.15.  Compared  to 
Figure  6.13,  the  prediction  slightly  overestimates  the  initial  size  of  the  first  delami¬ 
nation.  However,  the  overall  correlations  between  the  predictions  and  the  results  of 
the  X-Radiographs  are  fairly  good,  especially  at  higher  velocities. 

[04/452/  -  454/452/04] 

Figure  6.16  shows  the  X-Radiographs  of  [04/452/  —  454/452/04]  composites  result¬ 
ing  from  impact  by  an  impactor  with  a  mass  of  0.16  kg  at  three  different  velocities. 
The  lowest  velocity  corresponded  to  the  impact  velocity  threshold  of  the  compos¬ 
ites.  Only  a  long  0  degree  matrix  crack  in  the  bottom  ply  group  was  found  from  the 
X-Radiograph  as  well  as  by  the  eye  inspection  of  the  surface  of  the  specimen.  Ap¬ 
parently,  this  matrix  crack  was  the  initial  failure  mode.  No  delamination  was  found 
from  the  X-Radiograph  of  this  particular  specimen,  but  an  X-Radiograph  taken  from 
a  specimen  impacted  at  slightly  higher  velocity  shows  a  sign  of  delamination  near 
the  impacted  area  in  Figure  6.16.  Except  for  the  first  interface  between  the  0  degree 
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T300/976  [03/903/03/903/03] 


Mass  =  0. 1 6  kg  Velocity  =  6.70  m/s 


Figure  6.15:  The  predicted  delamination  sizes  of  [O3/9O3/O3/9O3/O3]  composites  cor¬ 
responding  to  the  impact  condition  given  in  Figure  6.13. 
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surface  ply  group  and  the  second  45  degree  ply  group,  delaminations  were  found  at 
the  rest  of  the  interfaces. 

The  measured  delamination  lengths  and  widths  as  a  function  of  the  impactor 
velocity  are  presented  in  Figure  6.17.  The  first  delamination  was  initiated  at  the 
last  interface  between  the  452/04  ply  groups  and  appeared  earlier  than  the  other 
delaminations.  Delaminations  with  relatively  smaller  sizes  were  also  found  at  the 
second  and  third  interfaces  of  the  laminates  at  higher  impactor  velocities.  Again, 
it  seems  that  the  first  delamination  governed  the  overall  delamination  size  and  was 
much  more  sensitive  to  the  impactor’s  velocity  than  the  other  delaminations.  No 
delamination  was  found  at  the  first  interface  between  the  surface  0  degree  ply  group 
and  the  second  45  degree  ply  group. 

The  predictions  of  the  delamination  length  and  width  of  the  composites  as  a 
function  of  the  impactor’s  velocity  are  also  shown  in  Figure  6.17  along  with  the  test 
data.  According  to  the  prediction,  matrix  cracking  at  the  bottom  0  degree  layer  group 
is  the  initial  failure  mode  which  initiated  the  first  delamination  at  the  last  interface 
between  452/04  ply  groups.  This  prediction  is  consistent  with  the  finding  from  the 
X-Radiograph  shown  in  Figure  6.16.  The  model  slightly  underestimated  the  impact 
velocity  threshold  of  the  laminate  and  the  threshold  velocity  corresponding  to  the 
occurrence  of  the  second  delamination,  but  the  predictions  are  correlated  fairly  well 
with  the  measured  first  and  second  delamination  sizes. 

However,  the  model  considerably  overestimated  the  velocity  threshold  correspond¬ 
ing  to  the  occurrence  of  the  third  delamination.  The  test  data  shows  that  both  the 
second  and  third  delaminations  initiate  at  roughly  the  same  time,  but  the  prediction 
indicates  that  the  third  delamination  will  occur  at  a  much  later  time  than  the  second. 
In  order  to  initiate  the  third  delamination,  the  model  predicts  that  matrix  cracking 
in  the  45  degree  ply  group  would  occur  much  later  than  the  matrix  cracking  in  the 
—45  ply  group  which  initiates  the  second  delamination. 

The  early  occurrence  of  the  matrix  cracking  in  the  45  degree  ply  group  from  the 
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T300/976  [04/452/-454/452/04] 


Mass  =  0.16  kg  Velocity  =  2.51  m/s 


Mass  =  0.16  kg  Velocity  =  3.90  m/s 


Figure  6.16:  Three  X-Radiographs  of  [04 /452 /  —  454/452/O4]  composites  in  series 
subjected  to  impact  by  a  spherical  nosed  projectile  at  three  different  velocities. 
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Velocity  (m/s) 


Figure  6.17:  Comparison  of  the  delamination  sizes  of  T300/976  [04/452/  —  454/452/04] 
composites  between  the  measurements  from  the  experiments  and  the  calculations 
based  on  the  model. 
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experiment  may  have  been  caused  by  the  loss  of  material  properties  in  the  neighbor¬ 
hood  of  the  impacted  area  due  to  the  combined  damage  of  the  cracking  in  the  —45 
degree  ply  groups  and  the  first  and  second  delaminations.  Such  interaction  among 
delaminations  and  matrix  cracks,  apparently  significant  in  [04/ 452/ — 454/452/ 04]  com¬ 
posites,  is  not  considered  in  the  model.  Therefore  a  progressive  failure  analysis  taking 
into  account  material  degradation  due  to  damage  would  be  necessary  to  accurately 
model  the  damage  interaction  during  impact. 

In  Figure  C.1S,  the  numerical  simulations  of  the  delaminations  of  the  test  results 
given  in  Figure  6.16  are  presented.  Overall,  the  predicted  delamination  sizes  are  quite 
consistent  with  the  results  of  X-Radiographs  shown  in  Figure  6.16. 

[04/454/  -  454/902]s 

Quasi-isotropic  laminates  were  also  tested  at  two  selected  velocities,  and  the  re¬ 
sults  of  the  test  are  shown  by  X-Radiographs  in  Figure  6.19.  Delaminations  were 
found  at  each  interface,  and  each  delamination  also  oriented  itself  along  the  fiber 
direction  of  the  bottom  ply  group  below  the  delaminated  interface.  The  predicted 
delamination  sizes  of  the  composites  corresponding  to  the  test  condition  are  also  pre¬ 
sented  in  Figure  6.19.  The  model  predicts  a  delamination  occurring  at  each  interface, 
consistent  with  the  test  results.  The  predicted  overall  size  and  shape  of  the  delami¬ 
nations  compared  to  the  X-Radiograph  pictures  both  given  in  Figure  6.19  are  fairly 
accurate.  The  predicted  delamination  size  and  shape  at  each  interface  within  the 
laminate  subjected  to  impact  at  velocity  of  7.8  m/s  are  presented  in  Figure  6.20. 


6.7  CONCLUSION 

An  investigation  was  performed  to  study  impact  damage  in  Graphite/Epoxy  lami¬ 
nated  composites  resulting  from  point-nose  impact.  A  model  was  developed  for  pre¬ 
dicting  the  initiation  of  the  damage  and  the  extent  of  the  delaminations  in  Graphite/Epoxy 
laminated  composites  resulting  from  the  impact.  Experiments  were  also  performed  to 
verify  the  model  and  the  computer  simulations.  Based  on  the  model,  a  user-friendly 
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T300/976  [04/45 2/ -454/452/04] 


Figure  6.18:  The  predicted  delamination  sizes  of  [04 / 452 /  —  454 / 4^2/ 04]  composites 
corresponding  to  the  impact  condition  given  in  Figure  6.16. 


135 


T300/976 

[04/454/-454/904/-454/454/04] 


Mass  =  0.16  kg  Mass  =  0.16  kg 

Velocity  =  7.8  m/s  Velocity  =  5.2  m/s 


Figure  6.19:  Comparison  of  delamination  sizes  and  shapes  of  [04/454/  —  454/902]s 
composites  between  X- Radiographs  taken  from  the  tested  specimen  and  predictions. 
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T300/976 
Mass  =  0.16  kg 
Velocity  =  7.8  m/s 
[04/454/-454/904/-454/454/04] 

It  II 

12  3  4 


Figure  6.20:  Predicted  delamination  sizes  and  shapes  of  [04/454/  —  454/902]s  com¬ 
posites  at  different  interfaces. 
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computer  code  was  developed.  Overall,  the  predictions  based  on  the  model  agreed 
with  the  test  data  very  consistently.  Based  on  the  study,  the  following  remarks  can 
be  made: 

1.  An  impact  velocity  threshold  exists  for  Graphite/Epoxy  laminated  composites 
below  which  no  delamination  occurs,  but  above  which  significant  damage  can 
be  produced. 

2.  Matrix  cracking  is  the  initial  failure  mode. 

3.  Matrix  cracking  in  composites  can  lead  to  interface  delaminations. 

4.  The  occurrence  of  the  critical  matrix  cracks  is  primarily  due  to  the  interlaminar 
tr  nsverse  shear  stress  and  transverse  in-plane  stress  of  the  concerned  laver. 

5.  Delamination  growth  is  governed  by  the  interlaminar  longitudinal  shear  stress 
and  transverse  in-plane  stress  in  the  layer  below  the  delaminated  interface  and 
the  interlaminar  transverse  shear  stress  in  the  upper  layers  directly  above  the 
interface. 
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Chapter  7 


COMPUTER  CODE 


During  the  investigation,  a  computer  code,  designated  as  “3DIMPACT,”  was  devel¬ 
oped  based  on  the  proposed  analysis.  The  user-friendly  code  can  be  used  to  provide 
the  following  information: 

1.  Stresses  and  strains  as  a  function  of  time  up  to  the  initial  damage  during  impact. 

2.  The  contact  force  as  a  function  of  time. 

3.  The  impact  velocity  threshold  corresponding  to  the  initial  impact  damage. 

4.  The  extent  of  final  delaminations  after  impact. 

The  material  properties,  laminate  ply  orientation,  thickness,  and  the  mass  of  the 
impactor  can  be  input  arbitrarily.  In  Appendix  D,  the  user-friendly  INPUT-OUTPUT 
of  the  “3DIMPACT”  code  is  presented. 
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Appendix  A 


THE  STIFFNESS  MATRIX 


On-axis  stress-strain  relationship  for  a  lamina,  is  given  by: 


where 
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The  off-axis  material  property  matrix  is  then  given  by: 
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Applying  plane  strain  condition  such  that 

£2  =  <4  =  <6  =  0 
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Appendix  B 


DETERMINATION  OF  THE  DISPLACEMENT  OF  IMPACTOR 


From  the  known  displacement  of  impactor  d\  at  time  t,  the  displacement  d2  at  time 
t  -+  At  can  be  calculated  as 


d7 


dr 


(B.l) 


where  the  v(t)  is  velocity  of  impactor.  From  the  known  velocity  t>i  at  time  t,  the 
velocity  at  time  r  is  calculated  as 


v(t)  =  v ,  ~  —  /  /(r)dr  (B.2) 

m  J  i 

where  the  /(r)  is  contact  force  between  impactor  and  plate.  The  contact  force  at 
time  t  is  given  and  equal  to  f\.  Assuming  the  contact  force  is  increasing  or  decreasing 
linearly  between  time  t  and  t  +  At,  the  contact  force  f2  at  time  t  +  At  can  be  expressed 
as 


/(t)  =  f\  +  ■  (r  -  0  :  t  <  t  <  t  +  At  (B.3) 

Eqs.  (B.2)  and  (B.3)  give 
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(B.4) 
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At  time  t  +  At,  the  velocity  u2  becomes 
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(B.5) 


By  combining  Eqs.  (B.l)  and  (B.4),  the  displacement  ofimpactor  <f2  at  time  t  +  At 
is  obtained  as  follows; 
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Appendix  C 


DATA  OF  POINT-LOADING  IMPACT  EXPERIMENTS 
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T300/976  [4547-458/454] 


2  1 


Mass 

(kg) 

mm 

■ 

De  lamination 
at  the  interface  1 
(cm) 

Delamination 
at  the  interface  2 
(cm) 

Length 

Width 

Length 

Width 

0.16 

2.31 

0.0 

0.0 

0.0 

0.0 

0.16 

3.04 

0.0 

0.0 

0.0 

0.0 

0.16 

3.15 

0.0 

0.0 

0.0 

0.0 

0.16 

3.49 

1.6 

0.5 

0.0 

0.0 

0.16 

3.63 

1.6 

0.5 

0.0 

0.0 

0.16 

4.06 

1.7 

0.5 

0.0 

0.0 

0.16 

4.15 

1.9 

0.5 

0.0 

0.0 

0.16 

4.24 

1.7 

0.5 

0.0 

0.0 

0.16 

4.87 

2.4 

0.9 

0.0 

0.0 

0.16 

4.88 

3.0 

1.5 

0.0 

0.0 

0.16 

5.89 

3.5 

1.6 

0.0 

0.0 

0.16 

7.26 

5.0 

2.4 

1.9 

1.0 

0.16 

9.02 

5.7 

2.5 

2.0 

1.1 

Table  C.l:  Measured  dalmination  sizes  of  T300/976  [45 4/  —  45g/454]. 
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T300/976  [03/903/03/903/03] 


H  t 


3  2  1 


Mass 

(kg) 

v  .  .  '  .  / 

Delamination 
at  the  interface  1 
(cm) 

Delamination 
at  the  interface  2 
(cm) 

Delamination 
at  the  interface  3 
(cm) 

Length 

Width 

Length 

Width 

Length 

Width 

0.16 

2.73 

EH 

0.0 

0.0 

0.0 

0.0 

0.0 

0.16 

3.16 

EH 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

3.22 

E9 

0.3 

0.0 

0.0 

0.0 

0.0 

0.16 

3.34 

1.0 

1.3 

0.0 

0.0 

0.0 

0.0 

0.16 

3.59 

0.8 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

3.60 

1.1 

0.5 

0.0 

0.0 

0.0 

0.0 

0.16 

3.68 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

3.93 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

4.00 

0.8 

1.6 

1.2 

Eh 

0.16 

4.03 

0.6 

1.1 

0.8 

na 

0.16 

4.12 

0.8 

1.5 

1.1 

0.5 

0.16 

4.18 

1.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.16 

4.72 

2.2 

0.8 

1.4 

0.8 

1.1 

0.5 

0.16 

5.35 

3.4 

1.1 

1.9 

1.1 

1.5 

0.7 

0.16 

6.70 

3.6 

1.1 

1.9 

1.2 

1.6 

0.9 

Table  C.2:  Measured  dalmination  sizes  of  T300/976  [O3/9O3/O3/9O3/O3]. 
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T300/976  [04/452/-454/452/04] 

t  t  t 

3  2  1 


Mass 

(kg) 

Test 

Velocity 

(m/s) 

Delamination 
at  the  interface  1 
(cm) 

Delamination 
at  the  interface  2 
(cm) 

Delamination 
at  the  interface  3 
(cm) 

Length 

Width 

Length 

Width 

Length 

Width 

0.16 

1.78 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.16 

2.25 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.16 

2.51 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.16 

2.63 

1.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.16 

2.64 

1.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.16 

2.86 

1.8 

0.2 

0.0 

0.0 

0.0 

0.0 

0.16 

2.92 

1.5 

0.2 

0.0 

0.0 

0.0 

0.0 

0.16 

3.33 

2.2 

0.3 

0.0 

0.0 

0.0 

0.0 

0.16 

3.36 

2.0 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

3.37 

2.5 

0.5 

0.0 

0.0 

0.0 

0.0 

0.16 

3.69 

2.5 

0.4 

0.0 

0.0 

0.0 

0.0 

0.16 

3.90 

3.1 

0.5 

1.5 

0.4 

1.5 

0.16 

4.34 

3.3 

0.8 

1.5 

0.5 

1.5 

0.16 

5.90 

4.3 

1.1 

2.3 

0.8 

2.1 

0.7  | 

Table  C.3:  Measured  dalmination  sizr*  of  T300/976  [04/452/  —454/452/04]. 
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Appendix  D 


USER-FRIENDLY  INPUT-OUTPUT  OF  “3DIMPACT”  CODE 


158 


•  * 
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D 

II 

MEM 

NM4 

PP  PP  AAAAA 

CC 

rr 
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D 

D 

II 

mm  m 
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PP  AAAAAAAAA  CC 

TT 
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II 
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PP  AA  AA 
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TT 

•  • 
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ODO 

mi 

EMM 

NMM  PPPP  AA  AA 

cccc 

IT 

»•  -  «  IMPACT  DAMAGE  Of  COtBOSITE  PLATES  »  - 

•  • 

••  THE  PURPOSE  Of  THIS  PROGRAM  IS  TO  PREDICT  THE 
••  DAMAGE  IN  COWOSITE  PLATES  SUBJECTED  TO  LOW-VLEOCITY 

••  FOREIGN  OBJECT  INPACT. 

*•  HYUNG  YUN  CHOI  4  FU-KUO  CHANG 

••  DEPARTMENT  Of  AERONAUTICS  AM)  ASTRONAUTICS 

••  STANFORD  UNIVERSITY 

••  STANFORD,  CA  94305 

••  (415)  123-3466 

•* 

••  CREATED:  JUttt.  1989 

••  MODIFIED:  NOVEPBER.  1990 


THIS  PROGRAM  FOLLOWS  THE  NUMERICAL  STRUCTURE 
OF  THE  'IMPACT*  BY  H-Y  T.  WU  t  G.  S.  SPRINGER 

«  COPYRIGHT  1990,  H.Y.  CHOI  4  F.K.  CHANG 


*«  IN  ORDER  TO  MINE  THE  INPUT  EASIER,  YOU  MAY  SAVE  THE 

••  INPUT  PARAMETERS  IN  THE  DATA  FILE  “INPUT.DAT*. 

"  THE  PROGRAM  HILL  INPUT  THE  NUBERS  AUTOMATICALLY . 

•  • 

••  IF  YOU  ARE  THE  FIRST  TINE  USER.  YOU  SHOULD  USE  THE 

"  INTERACTIVE  PROCEDURE  TO  INPUT  YOUR  DATA. 


••  WHEN  YOU  ARE  US  IMS  THE  INTERACTIVE  PROCEDURE,  THE  PROGRAM 

"  WILL  GENERATE  A  DATA  FILE  MNfS  •INPUT.DAT*  ACCORDING 

••  TO  YOUR  INTERACTIVE  :UOI. 

•  • 

••  ONCE  TW  IWUT  DATA  FILE  *:WUT.DAT*  HAS  BEEN  CREATED, 

••  YOU  CRY  INPUT  THE  DATA  CITHER  BY  THE  INTERACTIVE 

••  PROCEDURE  OR  DIRECTLY  FROM  THE  DATA  FILE. 


♦00  YOU  WANT  TO  INPUT  DATA  TOOM  THE  DATA  FILE  OR  DO  YOU  WANT  TO 
INPUT  DATA  INTERACTIVELY  ? 

ENTER  0  FOR  DATA  FILE 
ENTER  1  FOR  INTERACTIVELY 
0 


PROBLEM  STATQBNT: 

THE  PROBLEM  IS  SHOW  IN  THE  FOLLOWING  FIGURE: 
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0 


IPPACTOR 


( 

I 

1/ 

V 


BOUNDARY  2 


BOUNDARY  1  / 

/ 


X3  X2  /| 

I  /  II 

1/  II 

-  XI  /  /  BOUNDARY  3 

/  / 

/  / 

/  I 


I  / 


BOCACARY  4 


IN  THIS  PROGRAM,  THE  FOLLOWING  PAP  INTERS  ARE  SPECIFIED 
BY  THE  USER: 

(1)  I  VP  ACTOR: 

<A)  mss 

(B)  RADIUS 

(C)  DENSITY 

(D)  YOUNG'S  MODULUS 

(E)  POISSON'S  RATIO 

(F)  VELOCITY 


(2)  COWOSITE  PLATE: 

(A)  GECMETRY 

(I)  LENGTH 

(II)  WIDTH 
(lit)  THICKNESS 

(B)  PLY-ORIENTATIONS 

(C)  MATERIAL  PROPERTIES 

(D)  BOUNDARY  CONDITIONS  ON  BOUK2ARIE5  1,2.3,  AW  4 

(I)  SDVLY  SUPPORTED 

(II)  CLAMPED 

(III)  FREE 


OBJECTIVES: 

(1)  TO  FIND  THE  FOLLOWING  INTOWOMION: 

(A)  Q®  ACTOR 

(I)  VELOCITY  AS  A  rJCTION  OF  TIP* 

(II)  DISPLACXMBrT  AS  A  rUNCTICN  OF  TIP* 

(B)  C OPPOSITE  PLATE 

(I)  contact  force  as  a  function  or  tdc 

(ID  DISPLACOBarr  AS  A  FUNCTION  OF  POSITION  t  TOC 

(III)  STRESSES  AS  A  FUNCTION  OF  POSITION  4  TIME 

(IV)  STRAINS  AS  A  FUNCTION  OF  POSITION  4  TOC 


(2)  TO  FIM9  THE  SIZES  MC  LOCATIONS  OF  DAMAGES  FROM 
THE  FOLLOWING  FAILURE  31  ITER LA: 

(A)  miRIX  CRACK  CRITERION  FRW  CHOI  4  CHANG (1911) 

(B)  DE LAMINATION  C3IOWT*  CRITERION  FROM  CHOI  4  CHANG (1990) 


ANALYSIS: 

THE  STRESS  HISTORY  IS  CALCULATED  USING  A  TRANSIENT 
FINITE  ELSCNT  (CTHOD  TOMOLATED  ON  THE  BASIS  OF  THREE 
DIMENSIONAL  ELEtCNT  PESHES. 


near  instructions 


WHENEVER  YOU  SEE  THE  QUESTION  MUW  *  ?  *, 

YOU  MUST  INPUT  THE  DESIRED  VALUE. 

•ENTER  THE  NUOER  IDENTIFYING  THIS  RUN  (A  THREE  DIGITS  INTEGER)  ? 
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•THE  MESH  or  THE  PLATE  IS  16  BY  12  BY  4  (LENGTH,  WIDTH,  THICKNESS)  . 
DO  YOU  WANT  TO  CHANGE  THE  MESH  ? 

ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
0 


I  •FACTOR 


- THE  PROPERTIES  OF  THE  FOLLOWING  I  NT  ACTOR  MATERIALS 

ARE  BUILT  INTO  THE  PROGRAM: 


NAME 

MASS 

RADIUS 

DENSITY 

YOUNG'S  MODULUS 
POISSON'S  RATIO 


I  MAT.  NO.  1: 

I  STEEL  BALL 
I  0.842E-2 
I  0.635E-2 
I  7.8SE3 

I  2.1E11 

I  0.3 


MAI.  NO.  2:  I 

ALUMINUM  BALLI 
0.289E-2  I 
0. 635E-2  I 
2.70E3  I 
7.3E10  I 
0.3  I 


•DO  YOU  WANT  TO  USE  ONE  OF  THE  ABOVE  MATERIALS 
AS  I  (FACTOR  ? 

ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
0 

•ENTER  THE  MUSS  (KG)  7 
0.1610000E+00 

•ENTER  THE  RADIUS  (M>  7 
0.6350000E-02 

•ENTER  THE  DENSITY  (KG/CUBIC  M)  7 
0 . 7850000E+04 

•  ENTER  THE  YOUNG'S  MODULUS  (PA)  7 
0.21Q0000E+12 

•ENTER  THE  POISSON'S  RATIO  7 
0 . 3000000E+00 

- 1  (FACTOR'S  VELOCITY: 

•ENTER  THE  VELOCITY  (M/SECCND)  7 
0.5000000E+01 


COMPOSITE  PLATE 


-  BOCTCARY  COtOITIONS: 


•ENTER  THE  BOUtCARY  CONDITION  ON  BOUtCARY  1  7 
TYPE  0  POR  SIMPLY  SUPPORTED 
TYPE  1  FOR  CLMFEB 
TYPE  2  FOR  FRIZ 

1 

•ENTER  THE  BOUtCARY  C08CITI0N  ON  BOUtCARY  2  7 
TYPE  0  FOR  SDFLY  SUPPORTED 
TYPE  1  FOR  CLAIFED 
TYPE  2  FOR  FRIZ 

2 

•ENTER  THE  BOUtCARY  COCITION  ON  BOUtCARY  3  7 
TYPE  0  FOR  SDFLY  SUPPORTED 
TYPE  1  FOR  CUtVtD 
TYPE  2  FOR  FREE 
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1 


♦ENTER  THE  BOUWARY  CONDITION  ON  BOUWARY  4  ? 
TYPE  0  FOR  SIMPLY  SUPPORTED 
TYPE  1  FOR  CLAWED 
TYPE  2  FOR  FREE 
2 


REFERRING  TO  THE  FIGURE  SHOWN  ABOVE,  WE  DEFINE  THE 
FOLLOWING: 

(1)  THE  LENGTH  OF  THE  PLATE  IS  IN  THE  Xl-OIRECTION 

(2)  THE  WIDTH  OF  THE  PLATE  IS  IN  THE  X2 -DIRECTION 
!3>  THE  FIRST  PLY  IS  THE  ONE  CLOSEST  TO  THE  ItCACTOR 
(4)  THE  PLY  ORIENTATION  IS  DEFINED  TO  BE  POSITIVE  IN 

THE  COUNTERCLOCKWISE  DIRECTION  FROM  THE  Xl-AXIS 
WHO!  LOOKING  DOWNWARD  FROM  THE  DIRECTION  OF  THE 
IMPACT 


-  GEOMETRY: 

•ENTER  THE  LENGTH  OF  THE  PLATE  (M>  ? 
0.1016000E+00 

•ENTER  THE  WIDTH  OF  THE  PLATE  (M)  7 
0.7«20000E-01 


COORDINATES  OF  THE  (CSH  OF  THE  PLATE 


Xl-tCSH(l)-  0.0 

/ 

/  XI -MESH (  9)-  0.50SE-01 

I  /  / 

I  /  / 

-  X2-(C5H<  7).  0.3S1E-01 


— *  X2-4CSHU)-  0.0 


- THE  COORDINATES  OF  THE  (CSH  OF  THE  PLATE 

(XI  AW  X2  DIRECTION)  ARE  INPUT  BY  THE  USER. 

- PLEASE  XtffUT  THE  VALUES  OF  Xl-tCSH  AW  X2-4CSB. 

• . * . .  XI  COORDINATE  . . . . 

Xl-MSB(l)-  0.0;  ;  ;  ;  X1-4CSH(  9)-  LENGTH/ 2  »  0.50IE-01 

Xl-tCSB(l)-  0.0  (CENTER  OF  THE  PLATE) 

•ENTER  THE  VALOR  OF  Xl-tCSH  (  2)  ?  (POSITIVE  REAL  NUMBER) 
0.31730001-02 

•ENTER  THE  VALOR  OF  Xl-tCSH  (  3)  7  (POSITIVE  REAL  NUMER) 
O.C3SOOOOE-02 

•ENTER  THE  VALOR  OF  Xl-tCSH  (  4)  7  (POSITIVE  REAL  NIMER) 
0.93250001-02 

•ENTER  TOE  VALOR  OT  Xl-tCSH (  5)  7  (POSITIVE  REAL  NUMER) 
0.1270000E-01 

•ENTER  THE  VALOR  0 T  Xl-tCSH  (  «)  7  (POSITIVE  REAL  (AMBER) 
0.19030001-01 

•ENTER  THE  VALOR  OT  X1-+CSR<  7)  7  (POSITIVE  REAL  NUMBER) 
0.2S40000E-01 

•ENTER  THE  VALOR  OF  Xl-tCJH  (  ()  7  (POSITIVE  REAL  NIMSER) 
0.3110000E-01 

X1-4C3HI  9)-  O.SOK-Ol  (RIGHT  EDGE  OF  THE  PLATE) 

. . . . .  «  COORDINATE  * . . . 

X2-(C3i(l)-  0.0;  .*  ;  ;  X2-(CSH(  7).  WIDTH/2  -  0.3(1E-01 

X2-4CSH11)*  0.0  (MI DOLE  OF  THE  PLATE) 

•ENTER  THE  VALOR  OF  X2-4CSHI  2)  7  (POSITIVE  REAL  NUMER) 

0 . 31 73000E-02 
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IEMRR  THE  VALUE  OF  X2-*CSH(  3)  ?  (POSITIVE  REAL  NUMBER) 
0.63500001-02 

♦ENTER  THE  VALUE  OF  X2-tCSH(  4)  ?  (POSITIVE  REAL  NUMBER) 
0.9525000E-02 

♦ENTER  THE  VALUE  OF  X2-*CSH(  5)  ?  (POSITIVE  REAL  NUMBER) 
0 . 1270000E-01 

♦ENTER  THE  VALUE  OF  X2-tCSH<  6)  ?  (POSITIVE  REAL  NUMBER) 
0.2S40000E-01 

X2-MESHI  7)-  0.3HE-01  (TOP  EDGE  OF  THE  PLATE) 


DEFINITION  OF  PLY  PROPERTIES: 

EXX  -  LONGITUDINAL  YOUNG'S  MODULUS  (PA) 

EYY  -  TRANSVERSE  YOUNG'S  MODULUS  (PA) 

VXY  -  POISSON'S  RATIO  IN  THE  X-Y  DIRECTION 

VYZ  -  POISSON'S  RATIO  IN  THE  Y-Z  DIRECTION 

GXY  -  SHEAR  MOOULUS  IN  X-Y  DIRECTION  (PA) 

H  -  PLY  THICKNESS  <M) 

RHO  -  DENSITY  (KG/CUBIC  M) 

ACR  -  CRITICAL  IWJENTATION  (M) 

LT  -  LONGITUDINAL  TENSILE  STRENGTH  (PA) 

LC  -  LONGITUDINAL  CO»RESSIVE  STRENGTH  (PA) 

TT  -  TRANSVERSE  TENSILE  STRENGTH  (PA) 

TC  -  TRANSVERSE  COWRESSIVE  STRENGTH  (PA) 

IS  -  IOSIPESCU  SHEAR  STRENGTH  (PA) 


THE  PROPERTIES  OF  THE  FOLLOWING  MATERIALS  ARE  BUILT  INTO 
THE  PROGRAM  (PA  OR  M  OR  KG/CUBIC  M)  : 


MATERIAL  NO.  1: 
MATERIAL  NO.  2: 
MATERIAL  NO.  3: 
MATERIAL  NO.  4: 
MATERIAL  NO.  S: 
MATERIAL  NO.  6: 
MATERIAL  NO.  7; 
MATERIAL  NO.  I: 


CFRP  T300/976 
CFRP  T300/934 
CFRP  T300/NS20I 
BFRP  B(4)  /NS50S 
CFRP  AS/3S01 
GFRP  E -GLASS /EPOXY 
KFRP  KEV  49/EPOXY 
CTRTP  AS4/PEEX 


THE  PROPERTIES  ARE: 


1 

MAT.  NO.  1 

1  MAT.  NO.  2 

1  MAT.  NO.  3 

1  MM.  NO.  4 

1 

1  EXX 

I.S6E1I 

1  1.4SE11 

1  l.tlEU 

i  2.04E11 

1 

1  EYY 

9.09E09 

1  9.99E09 

1  I.03E10 

1  1.85E09 

1 

1  VXY 

0.221 

1  0.30 

1  0.21 

1  0.23 

1 

1  VYZ 

0.400 

1  0.30 

1  0.21 

1  0.23 

t 

1  GXY 

6.96Z09 

1  S.61E09 

1  7.17E09 

1  5.S9E09 

1 

1  H 

1.44E-4 

1  1 . S9E-4 

1  I.25E-4 

1  1.25E-4 

! 

1  ABO 

1.S4E03 

1  1.54E03 

1  1.60E03 

1  2.00E03 

1 

1  ACR 

6.03E-5 

1  ♦ . 03E-5 

1  1.03E-5 

1  I.03E-5 

1 

I*  LT 

1.S2E09 

i  1.7IE09 

1  I.S0E09 

I  1 .26X09 

1 

1  LC 

1.S9E09 

1  1.73E09 

I  1.50E09 

1  2.S0E09 

! 

1  TT 

4.4SE07 

1  5.51E07 

1  4.00E07 

1  6.10E07 

I 

1  TC 

2.S2E0I 

1  2.94E0* 

1  2.46E0I 

1  2.02E01 

1 

1  IS 

3.S6E07 

1  3.29E07 

1  6.90E07 

1  6.70E07 

1 
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1  1 

MAI.  NO.  S 

1  MAI.  NO.  6 

1  MAT.  NO.  7 

1  MkT.  NO.  8  1 

1  EXX 

1.31E11 

1  3.86E10 

1  7.59E10 

I  1.34E11  1 

1  EYY 

S.96E09 

1  8.26E09 

1  S.S0E09 

1  B.B9E09  1 

1  VXY  1 

0.30 

1  0.26 

1  0.34 

1  0.28  1 

1  VYZ  1 

0.30 

1  0.26 

1  0.34 

1  0.28  1 

1  GXY  1 

7.10E09 

1  4.14E09 

1  2.30E09 

1  S.I0E09  1 

1  H  1 

1.2SE-4 

t  1.2SE-4 

1  1.2SE-4 

1  1.2SE-4  1 

1  RHO  | 

1.60E03 

1  1.80E03 

1  1.46E03 

I  I.6IE03  1 

i  ACR  | 

B.03E-5 

1  1.02E-4 

1  8.03E-S 

1  8.03E-5  1 

i  LT  1 

1.4SE09 

1  I.06E09 

1  1.40E09 

1  2.13E09  1 

1  LC  1 

1.45E09 

1  6.10E08 

1  2.3SE08 

1  I.10E09  1 

1  TT  1 

S.20E07 

1  3.10E07 

1  1.20E07 

1  7.99E09  i 

1  TC 

2.06EOI 

1  1.17E08 

1  S.30E07 

1  2.OOEO0  1 

1  IS  1 

9.29E07 

1  7.19E07 

1  3.40E07 

1  I.60E08  1 

-  YOU  MAY  USE  THE  PROPERTIES  OF  ANY  OTHER  COMPOSITE. 

»00  YOU  WANT  TO  ENTER  THE  PROPERTIES  OF  ADDITIONAL  MATERIALS  ? 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

0 


PLIES  »Y  BE  ENTERED  INDIVIDUALLY  OR  IN  GROUPS. 

THE  FIRST  PLY  GROUP  IS  THE  ONE  CLOSEST  TO  THE  IM>ACTOR. 
THE  PLY  ANGLE  IS  MEASURED  FROM  THE  Xl-AXIS  (SEE  FIGURE 
ABOVE)  IN  THE  COUNTERCLOCKWISE  DIRECTION 
(LOOKING  DOWNWARDS  FRCM  THE  DIRECTION  OF  THE  IMPACT)  . 


- SPECIFY  THE  PLY  ORIENTATION  OF  THE  PLATE. 

FOR  EXAMPLE,  IS  LAYERS,  4  INTERFACES,  5  PLY  GROUPS: 
(0(3)/ 90  (3)/0(3)/ 90(31/0(3)  ) 

[THETX1  (Nil  /THETA2  (N2) / . . . . ] 


♦  ENTER  THE  NUMBER  OF  LAYERS  OF  THE  PLATE  ( INTEGER)  ? 

IS 

♦ENTER  THE  NUMBER  OF  PLY  GROUPS  OF  THE  PLATE  (INTEGER)  ? 

5 

♦ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  NO.  1  ( INTEGER)  ? 
3 

♦ENTER  THE  ORIENTATION  IN  PLY  GROUP  NO.  1  (DEGREE)  ? 

o.oooooooe*oo 

♦ENTER  THE  NOOER  OF  PLIES  IN  PLY  GROUP  NO.  2  (INTEGER)  ? 
3 

♦ENTER  THE  ORIENTATION  IN  PLY  GROUP  NO.  2  (DEGREE)  7 
0 . 90000008*02 

♦ENTER  THE  MMKR  OF  PLIES  IN  PLY  GROUP  NO.  3  (INTEGER)  7 
3 

♦ENTER  THE  ORIENTATION  IN  PLY  GROUP  NO.  3  (DEGREE)  7 
0 .  OOOOOOOE+OO 

♦ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  NO.  4  (INTEGER)  7 
3 

♦ENTER  THE  ORIENTATION  IN  PLY  GROUP  NO.  4  (DEOTEE)  7 
0.9000000E+02 

♦ENTER  THE  NUMER  CT  PLIES  IN  PLY  GROUP  NO.  S  (INTEGER)  7 
3 

♦ENTER  THE  ORIENTATION  IN  PLY  GROUP  NO.  S  (DEGREE)  7 
O.OOOOOOOE*OQ 


- SINCE  THE  WSH  OF  THE  PLATE  IS  1«  BY  12  BY  4 

(LENGTH,  WIDTH,  THICKNESS).  THERE  ARE  4  GROUPS 
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THROUGH  THE  THIOOCSS  OF  THE  PLATE. 

♦DO  YOC  WANT  ALL  THE  PLIES  IN  EACH  GROUP 
HAVE  DIFFERENT  ORIENTATIONS  7 
ENTER  0  FOR  NO 
ENTER  1  FOR  TES 
1 

♦ENTER  THE  MATERIAL  NUMBER  GIVEN  ABOVE  TO  SPECIFY  THE 
MATERIAL  USED  IN  GROUP  NO.  1  ( INTEGER)  7 
1 

♦ENTER  THE  NUMBER  OF  PLIES  IN  GROUP  NO.  1  (INTEGER)  ’ 

3 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  1 
IN  GROUP  NO.  I  (DEGREES)  7 

O.OOOOOOOE+OO 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  2 
IN  GROUP  NO.  1  (DEGREES)  7 

o.oooooooe*oo 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  3 
IN  GROUP  NO.  1  (DEGREES)  7 

o.oooooooe*oo 

♦ENTER  THE  MATERIAL  NUMBER  GIVEN  ABOVE  TO  SPECIFY  THE 
MATERIAL  USED  IN  GROUP  NO.  2  (INTEGER)  7 
1 

♦ENTER  THE  NUMBER  OF  PLIES  IN  GROUP  NO.  2  (INTEGER)  ’ 

4 


♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  1 
IN  GROUP  NO.  2  (DEGREES)  7 
0 . 9000000E+02 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  2 
IN  GROUP  NO.  2  (DEGREES)  7 
0.9000000E*02 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  3 
IN  GROUP  NO.  2  (DEGREES)  7 
0.9000000E-02 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  4 
IN  GROUP  NO.  2  (DEGREES)  7 
0.0000000E+00 

♦ENTER  THE  MATERIAL  NUMBER  GIVEN  ABOVE  TO  SPECIFY  THE 
MATERIAL  USES  IN  GROUP  NO.  3  ( INTEGER)  7 
1 

♦ENTER  THE  NIMER  OF  PLIES  IN  GROW  NO.  3  1  INTEGER)  7 

4 

♦ENTER  THE  PLY-ORIENTATION  OF  TW  PLY  NO.  1 
IN  GROUP  NO.  3  (DEGREES)  7 
0.0000000£*00 

♦ENTER  THE  PLY-ORIENTATION  OF  Tm  PL T  NO.  2 
IN  GROUP  NO.  3  (DEGREES)  7 
0. 00000001*00 

♦ENTER  THE  PLY-ORIENTATION  OT  T7«  PLT  NO.  3 
IN  GROUP  NO.  3  (DEGREES)  7 
0.90000002*02 

♦ENTER  THE  PLY -ORIENTATION  CF  DK  PLT  NO.  4 
IN  StOUP  NO.  3  (DEGREES)  7 
0.9000000E*02 

♦ENTER  THE  MATERIAL  NUMBER  GIVEN  ABOVE  TO  SPECIFY  THE 
MATERIAL  USES  IN  GROUP  NO.  4  I  INTEGER)  7 
1 

♦ENTER  THE  NOOER  OF  PLIES  IN  GROUP  NO.  4  (INTEGER)  7 

4 

♦ENTER  THE  PLY -ORIENTATION  OT  THE  PLY  NO.  1 
IN  3WUP  NO.  4  (DEGREES)  7 
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0.9000000E+02 

*  ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  2 
IN  GROUP  NO.  4  (DEGREES)  ? 

0.0000000E+00 

♦ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  3 
IN  GROUP  NO.  4  (DEGREES)  ? 

o.oooooooe*oo 

♦  ENTER  THE  PLY-ORIENTATION  OF  THE  PLY  NO.  4 
IN  GROUP  NO.  4  (DEGREES)  ? 

o.oooooooe*oo 


SUMMARY  OF  THE  INPUT  DATA 


Iff  ACTOR: 

MASS 

RADIUS 

DENSITY 

YOUNG'S 

POISSON 


MODULUS 
S  RATIO 


VELOCITY 


0.16100E+00  KG 

0 . 63500E-02  M 

0 .  78500E*04  KG/ (CUBIC  M) 

0.210001*12  PA 

0 . 30000E+00 

0 . SC000E+01  M/SECOND 


COMPOSITE  PLATE: 
BOUNDARY  1  IS  CLAMPED 
BOUM3ARY  2  IS  FREE 
BOUMJARY  3  IS  CLAMPED 
BOUWAAY  4  IS  FREE 


-  GECPCTRICAL  CONFIGURATION: 

LENGTH  -  0.101S0E*0Q  M 

WIDTH  -  0.76200E-OI  M 


-  GROUP  NO.  1  : 

MATERIAL  NUKBER  • 
PLIES  NUMBER 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 

-  GROUP  NO.  2  : 

MATERIAL  NUMBER  - 
PLIES  NUMBER 

ply-orientation  of 

PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 

-  GROOP  NO.  3  : 

MATERIAL  NUKBER  - 
PLIES  NUMBER 
PLY-ORIENTATION  of 
PLY-ORIENTATION  OF 

ply-orientation  of 

PLY-ORIENTATION  OF 

-  GROW  NO.  4  : 

MATERIAL  NUKBER  - 
PUBS  NQMBER 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 
PLY-ORIENTATION  OF 


1 

3 


NO. 

1 

PLY  ■ 

0 . OOOOOE-OO 

DEGREES 

NO. 

2 

PLr  • 

O.OOOOOE *00 

DEGREES 

NO. 

3 

PLY  - 

O.OOOOOE-OO 

DEGREES 

1 

4 


NO. 

1 

PLY  « 

0.90000E-02 

degrees 

NO. 

2 

PLY  » 

0 .900002*02 

DEGREES 

NO. 

3 

PLY  - 

0.90000E-02 

DEGREES 

NO. 

4 

PLY  - 

O.OOOOOE+OO 

DEGREES 

I 

4 


NO. 

1 

PLY  - 

Q.OOOOOE+OO 

DEQUES 

NO. 

2 

PLY  - 

0 . OOOOOE-OO 

DEQUES 

NO. 

3 

PLY  • 

0 . 90000E-02 

DEGREES 

NO. 

4 

PLY  - 

0.90000E«-02 

DEGREES 

I 

4 


NO. 

1 

PLY  - 

0 . 90000E-02 

OECREES 

NO. 

2 

PLY  - 

0 . OOOOOE+OO 

DEQUES 

NO. 

3 

PLY  - 

O.OOOOOE *00 

DEGREES 

NO. 

4 

PLY  - 

0.00000E*00 

DEGREES 

INTEGRATION  TOC  STEP 


- THE  APPROX  MUZ  TO*  TO  REACH  THE  MAXIMUM  CONTACT  FORCE  IS: 

Til*  -  O.CE-3  SECONDS 

IMPACT  COICITION  IS: 

UfACTOR  MASS:  0.1(  KG 
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IteACTOR  VELOCITY:  5  M/SEC 
PLATE  THICKNESS:  0.0022  M 


IT  IS  RECCPtCWED  TO  CHOOSE  120  TIJC  STEPS  OF  INTEGRATION 
FOR  THE  END  OF  CONTACT. 

AS  IMPACTOR  MASS  INCREASED  OR  PLATE  BECOfrE  LESS  STIFF ER, 
THE  TIME  TO  REACH  THE  MAXIMUM  CONTACT  FORCE  WOULD 
BE  LONGER. 


•TIME  STEPS  OF  INTEGRATION: 


IF  NO, 

♦ENTER  THE  TOTAL  INTEGRATION  TI(C  TO  BE  USED  (INTEGER)  ? 

0 . 1200000E-02 

♦ENTER  THE  NOSER  OF  TIME  STEPS  TO  BE  USED 
(INTEGER)  ? 

100 

- FOR  THE  PRESENT  PROBLEM,  WE  HAVE: 

FINAL  TIME  OF  INTEGRATION  -  0.1200000E-02  SECONDS 

ONE  TK  STEP  -  0 .1200000E-04  SECONDS 


OUTPUT 


- THE  FOLLOWING  OUTPUTS  CAN  BE  GENERATED  BY  THIS  CODE. 

OUTPUT  GROUP  41: 

MATRIX  CRACK  IN  EACH  PLY  AT  ANY  TIME 

DELAMINATION  AREA  AMD  SHAPE  ALONG  THE  ANY  INTERFACE  AT  ANY  TIMS 
OUTPUT  GROUP  #2: 

CONTACT  FORCE  BETWEEN  PLATE  AND  IMPACTOR  AT  ANY  TIME 
DISPLACEMENT  OF  THE  IMPACTOR  AT  ANY  TIME 
VELOCITY  OF  THE  IMPACTOR  AT  ANY  TIME 


OUTPUT  GROUP  »3: 

OISPLACnCNT  OF  ANY  POINT  AT  ANY  TIME 
VELOCITY  OF  ANY  POINT  AT  ANY  TIJC 
STRAINS  OF  ANY  POINT  AT  ANY  TIME 
STRESSES  OF  ANY  POINT  AT  ANY  TIW 


NOTE: 

1)  THE  OUTPUT  OF  GROUP  #3  MIGHT  OCCUPY  LARGE  SPACES  OT 
STORAGE.  IT  YOU  DO  NOT  REALLY  MS)  THIS 
IWCRMttlON,  THEN  DO  NOT  ASK  FOR  IT  TO  BE  CALCULATED. 

2)  WHEN  INPUTTING  THE  PLY  MISERS.  THE  FIRST  PLY  IS  THE  ONE 
THAT  IS  CLOSEST  T3  THE  Dff ACTOR. 


—OUTPUT  GROUP  II: 

- MATRIX  CRACK  AM)  DELAMINATION  SHAPE/AREA  FROM  THE 

FAILOm  CRITERIA  PROPOSED  BY  CHOI  AM)  CHANG 


THE  OUTPUT  or  DAMAGE  SHAPE/AREA  IS  STORED  AND  CAN  BE  PRINTED 
FROM  THE  FILE  * DAMVZ.DAT*  AFTER  THE  RUNNING  OF  THE  PROGRAM. 

♦oo  you  twrr  to  calculate  group  no.  i  ? 

ENTER  o  for  no 
ENTER  1  FOR  YES 

1 


MATRIX  CRACK 
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■THE  LAMINATED  PLATE  OF  THIS  PROBLEM  HAS  5  PLY  GROUPS. 
THUS,  MATRIX  CRACX  CRITERION  WOULD  BE  APPLIED  AT  THOSE 
S  DIFFERBiT  PLY  GRIOUPS  AT  EACH  THC  STEP. 


DELAMINATION 


- DELAMINATION  IS  ASSCPED  TO  OCCUR  ALONG  THE  INTERFACE  BETWEEN 

TWO  ADJACENT  PLIES  WHEN  A  MATRIX  CRACX  OCCURS  WITHIN  ONE  OR  BOTH 
ADJACENT  PLIES. 

THE  LAMINATED  PLATE  OF  THIS  PROBLEM  HAS  4  INTERFACES. 

THUS,  DELAMINATION  GROWTH  CRITERION  WOULD  BE  APPLIED  AT  THE 
INTERFACE  WHERE  THE  MATRIX  CRACK  CRITERION  IS  SATISFIED  FOR  ANY 
ADJACENT  PLIES. 

- ©EXAMINATION  IS  ILLUSTRATED  BY  *•'  AM}  UNDAMAGE  AREA  IS  REMAINED 

WITH  V  IN  THE  OUTPUT  FILE  DAMAGE.DAT 

- IF  THE  PLATE  IS  SYWCTRIC  WITH  RESPECT  TO  THE  XI  OR  X2  AXIS. 

THE  DELAMINATION  AREA  AND  SHAPE  WOULD  BE  ALSO  SYMMETRIC. 

THEREFORE,  ONLY  HALF  OR  QUARTER  PORTION  OF  THE  PLATE  IS  NEEDED  TO 
CALCULATE  THE  OELAMI NATION  AREA  AM}  SHAPE. 

CHOOSE  ONE  OF  FOLLOWING  NUTBERS  FOR  DELAMINATION  SHAPE  AND  AREA. 

(1)  FULL  SIZE 

(2)  HALF  SIZE 

(3)  QUATER  SIZE 

•DO  YOU  WANT  TO  HAVE  1  OR  2  OR  3  7 
1 

- FOR  MOST  OF  CAESES,  DELAMI NATIONS  ARE  CONFINED  WITHIN  THE  CENTER 

REGION  OF  THE  PLATE.  FOR  THE  CLOSER  VIEW  OF  THE  PREDICTED  D EXAMINATIONS 
FORM  THE  SCREEN,  YOU  KAY  ZOOM  IN  THE  IS  ONLY  NEEDED  TO  APPLY 
DELAMINATION  GROWTH  CRITERION  LOCAL  AREA  BY  CHOOSING  SCALING 
PARAMETER,  R.  FROM  0  TO  1 .  IF  R-l,  THE  ORIGINAL  SCALE. 

•ENTER  THE  SCALE  FOR  DELAMI ANTICN  AREA  CALCULATION  <0<  R  <-1.0>  ? 

0 .7000000E+00 

- DELAMINATION  AREA  AM)  SHAPE  WILL  PRINTED  DURING  INTERMEDIATE  TIME. 

•AFTER  HOW  MANY  TOC  STEPS  00  YOU  WANT  THE  RESULTS 
TO  BE  PRINTED  7 
ENTER  AN  INTEGER 
1 

- AN  EMPIRICAL  CONSTANT  (DO  IS  NEEDED  FOR  DELAMINATION  GROWTH 

CRITERION. 

RECOM3DZZ}  CONSTANT  D«  FOR  GRAPHITE/EPOXY  T300/97S  MATERIAL 
SYSTEM  IS  1 .1. 

CONSTANTS  FOR  THE  OTHER  MATERIAL  SYSTEMS  ARE  NOT  CURRENTLY  AVAILABLE. 
FOR  FUTHER  INFOfMATION,  CONTACT  DR.  CHANG  AX  THE  ADDRESS  GIVEN  ON 
THE  FIRST  PAGE. 

•ENTER  THE  CONSTANT  FOR  DELAMINATION  GROWTH  CRITERION  7 
0 .  ItOOOOOE*01 


- OUTPUT  GROUP  *2: 

—CONTACT  FORCE,  II*  ACTOR  DISPLATOCNT/VELOCITY 

THE  OUTPUT  or  CONTACT  FORCE,  IKPACTOR  DISPLACEMSNT/VELOCITY 
IS  STORES  AM)  CAN  BE  PRINTS)  FROM  THE  FILE 
"3DGP2.DAT"  AFTER  THE  ROWING  OF  THE  PROGRAM. 


•DO  YOU  WANT  TO  CALCULATE  GROUP  NO.  2  7 
ENTER  0  FOR  NO 
ENTER  I  FOR  YES 

1 


IF  YES, 

•DO  YOU  WANT  TO  CALCULATE  THE  CONTACT  FORCE  BETWEEN  THE  PLATE 
AM)  THE  IWACTCR  DURING  THE  CONTACT  ? 

ENTER  0  FOR  NO 
ENTER  l  FOR  YES 

1 
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ido  you  want  to  calculate  the  displacd«nt  of  the  impactor 

DURING  THE  CONTACT  7 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
1 

IDO  YOU  NANT  TO  CALCULATE  THE  VELOCITY  OF  THE  IMPACTOR 
DURING  THE  CONTACT  7 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
1 


- OUTPUT  GROUP  13: 

- DISPLACEMENT,  VELOCITY,  STRAIN,  AND  STRESS  OF  THE  PLATE 

- REMINDER:  DO  NUT  CALCULATE  GROUP  13  UNLESS  IT  IS 

NECESSARY  BECAUSE  THE  OUTPUT  MIGHT  OCCUPY 
MANY  SPACES  TO  STORE. 

IDO  YOU  WANT  TO  CALCULATE  GROUP  NO.  3  7 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
1 


IF  YES, 

. .  DISPLACEMENT  AND  VELOCITY 


- THE  COORDINATES  OF  THE  PLATE  ARE: 

-O.SOtE-OI  (BOUNDARY  1)  <  XI  <  0.508E-01  (BOUNDARY  3) 

-0. 31  IE -01  (BOUNDARY  4)  <  X2  <  0.381E-01  (BOUNDARY  2) 

-0. 10IE-02  (BOTTOM)  <  X3  <  0.10IE-02  (TOP) 

PLEASE  ENTER  COORDINATES  MITHIN  THESE  LIMITS. 

THE  OUTPUT  OF  DISPLACEMENT /VELOCITY  OF  THE  PLATE  IS  STORED 
AMD  CAN  BE  PRINTED  FROM  THE  FILE  '3DOV.DAT* 

AFTER  THE  RUNNING  OF  THE  PROGRAM. 

IDO  YOU  KANT  TO  CALCULATE  THE  DISPLACQdTS  AM)  VELOCITIES 
IN  THE  PLATE  7 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 
1 


IF  YES. 

IAT  HO*  MANY  POINTS  DO  YOU  WANT  THE  DISPLACEJCNTS  AND 
VELOCITIES  TO  BE  CALCULATED  ? 

ENTER  AN  INTECXR 
I 

IWOT  THE  COORDINATES  OT  EACH  POINT  OM  BY  ONE. 

•ENTER  THE  XT-COORDINATE  OF  POINT  NO.  1  7 
0 . 1000000E-03 

•ENTER  THE  X2 -COORDINATE  OT  POINT  NO.  1  7 
0 . 10000001-03 

•ENTER  THE  XXOOHD INATE  OF  POINT  NO.  1  7 
0 . 10000001-05 

the  Dispucaerrs  aw  velocities  *ill  be  printed 
DQR2NG  INIEMDIAXE  TIME. 

•AFTER  BOR  MANY  TO«  STEPS  X  TOO  RANT  THE  RESULTS 
TO  BE  PRINTED  7 
ENTER  AN  INTK» 

1 


STRAINS  A7C  STRESSES 


- THE  COORDINATES  OF  THE  PLATE  ARE: 

-0. SO IE -01  (BOUNDARY  1)  <  XI  <  0.S0IE-01  (BOUNDARY  3) 

-0.3I1E-01  (BOOWARY  4)  <  X2  <  0.3I1E-01  (BOUKDARY  2) 

-0. 10IE-42  (BOTTOM)  <  X3  <  0.10IE-02  (TOP) 

SINCE  THE  STRAINS  AW  STRESSES  ARE  NOT  ALWAYS  CONTINUOUS 
THROOa  THE  THIOOCSS,  THE  PLY  NLMKR  (COUNTED  FROM  TOP) 


169 


>«JST  BE  IWUT  TOGETHER  WITH  THE  X3-COORDINATE. 

PLEASE  ENTER  COORDINATES  WITHIN  THESE  LIMITS. 

THE  OUTPUT  OF  THE  STRAINS  OF  THE  PLATE  IS  'RED  AND  CAN 
BE  PRINTED  FROM  THE  FILE  *  3D  STRAIN .  DAT*  AF  '  .  THE  RUNNING 
OF  THE  PROGRAM. 

THE  OUTPUT  OF  THE  STRESSES  OF  THE  PLATE  IS  STORED  AM)  CAN 
BE  PRINTED  FROM  THE  FILE  ■  3DSTRESS . DAI*  AFTER  THE  RUNNING 
OF  THE  PROGRAM. 


STRAINS 


♦DO  YOU  WANT  TO  CALCULATE  THE  STRAINS  ? 

ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

1 

IF  YES, 

♦AT  HOW  MANY  POINTS  DO  YOU  WANT  THE  STRAINS  TO  BE  CALCULATED  ? 

1 

INPUT  THE  COORDINATES  OF  EACH  POINT  ONE  BY  ONE. 

♦ENTER  THE  Xl-COOSINAIE  OF  POINT  NO.  1  ? 

0.1000000E-03 

♦ENTER  THE  X2 -COORDINATE  OF  POINT  NO.  I  7 
0 . I000000E-03 

♦ENTER  THE  X 3-COORDINATE  OF  POINT  NO.  1  ? 

0 . 1000000E-05 

♦ENTER  THE  PLY  NUMBER  IN  WITCH  THIS  POINT  IS  LOCATED  ? 

t 

THE  STRAINS  WILL  BE  PRINTED  DURINS  INTERMEDIATE  TIME. 

♦AFTER  HON  MANY  TIME  STEPS  DO  YOU  WANT  THE  RESULTS 
TO  BE  PRINTED  7 
ENTER  AN  INTEGER 
1 


STRESSES 


♦00  YOU  WANT  TO  CALCULATE  THE  STRESSES  T 
ENTER  0  FOR  NO 
ENTER  1  FCR  YES 
1 


IF  YES. 

♦AT  HON  (ANY  POINTS  DO  YOU  WANT  THE  STRESSES  TO  BE  CALCULATED  7 
ENTER  AN  INXECZR 
I 


INPUT  THE  COORDINATES  OF  EACH  POINT  0(C  BY  ONE. 


♦ENTER  THE  Xl-COOHDINAIE  OF  POINT  NO.  I 
0 . 1000000E-03 

♦ENTER  THE  X2-C00H0 INAIZ  OF  POINT  NO.  1 
0 . 1000000E-03 

♦ENTER  THE  X3-COOHDINATE  OF  POINT  NO.  I 
0 . lOOOOOOE-OS 

•ENTER  THE  PLY  NCMER  IN  WITCH  THIS  POINT 


7 

7 


IS  LOCATED  7 


THE  STRESSES  WILL  BE  PRINTS  DURING  INTER7C0IAIE  TO* . 

♦AFTER  ROW  MANY  TO*  STEPS  00  YOU  WANT  THE  RESULTS 
TO  BE  PRINTS  7 
ENTER  AN  TNIEOR 

I 


QC  OF  IWH7T  AND  OUTPUT  INSTRUCTION 


♦  ♦♦♦»♦«  PROGRAM  STARTS  RUNNING 


Til*  STEP- 
Til*  STEP- 
TII*  STEP- 


1; 

1; 

2; 


TOe- 

TO*- 

TD*« 


0.120000E-04 
0.120000E-04 
0 • 240000E-04 


PROGRAM  IS  RUNNING 
PKX3LAM  IS  RUNNING 
PROGRAM  IS  RUM1TNG 
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TItC 

STEP- 

92; 

THE- 

0.110400E-02 

PROGRAM 

IS 

RUNNING 

THE 

STEP- 

93; 

TIME- 

0.111600E-02 

PROGRAM 

IS 

RUNNING 

TIME 

STEP- 

94; 

TIME- 

0.112807E-02 

PROGRAM 

IS 

RUNNING 

THE 

STEP- 

95; 

TIME- 

0.114000E-02 

PROGRAM 

IS 

RUNNING 

TItC 

STEP- 

96; 

TIME- 

0. 115200E-02 

PROGRAM 

IS 

RUNNING 

THE 

STEP- 

97; 

TIME- 

0. 116400E-02 

PROGRAM 

IS 

RUNNING 

THE 

STEP- 

98; 

TIME- 

0.117600E-02 

PROGRAM 

IS 

RUNNING 

TIME 

STEP- 

99; 

TIME- 

0.UI800E-02 

PROGRAM 

IS 

RUNNING 

THE 

STEP- 

100; 

THE- 

0.120000E-02 

PROGRAM 

IS 

RUNNING 

ttttl 

III  PROGRAM  FINISHED 

RUNNING  *»***( 

II 

>  the 

OUTPUT 

OF  ' 

THE  CONTACT  FORCE,  IMPACT 

BODY 

0 1 SPLACEMENT/VELOC ITY  IS  STOREDTN  FILE  -3DGP2.DAT* 


THE  OUTPUT  OF  THE  D I SPLACEJCNTS/ VELOCITIES  OF  THE  PLATE 
IS  STORED  IN  FILE  ■3DDV.DAX" 


**  THE  OUTPUT  OF  THE 
IS  STORED  IN  FILE 


STRAINS  OF  EVERT  POINTS 
■3DSTHAIN.DAT* 


••  THE  OUTPUT  OF  THE 
IS  STORED  IN  FILE 


STRESSES  OF  EVERY  POINTS 
•3DSTRESS.DAT* 


THE  DATA  STORED  IN  THESE  FILES  CAN  BE  MANIPULATED 
(PRINTED  .PLOTTED.  *tc. )  USING  THE  CCMAND  APPROPRIATE 
TO  THE  OPERATING  SYSTEM  OF  YOUR  COMPUTER 

.  END  . . . . . 
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••  AT  TI*  STEP  1 

MATRIX  CRACK 

PLY  GROUP  •  1  ( 
PLY  GROUP  I  2  < 
PLY  GROUP  I  3  ( 
PLY  GROUP  »  4  ( 
PLY  GROUP  »  5  ( 

DELAMINATION 


(  0.1200E-04)  •• 


O.OOOOE+OO) 

NO 

0 . 9000E+02) 

NO 

0.0000E+00) 

NO 

0.9000E+02) 

NO 

0.0000E+00) 

NO 

INTERFACE 

1 

1 

(BETWEEN  PLY 

GROUP 

» 

1 

AW) 

* 

2) 

NO 

INTERFACE 

# 

2 

(BETWEEN  PLY 

GROUP 

1 

2 

AW) 

1 

3) 

NO 

INTERFACE 

♦ 

3 

(BETWEEN  PLY 

(SOUP 

* 

3 

AW) 

» 

4) 

NO 

INTERFACE 

» 

4 

(BETWEEN  PLY 

GROUP 

1 

4 

AW) 

* 

5) 

NO 

*•  AT  TKC  STEP  100  (  0.1200E-02>  •* 

MATRIX  CRACK 


PLY 

CROUP 

» 

1  ( 

0 . 0000E+00) 

NO 

PLY 

(SOUP 

« 

2  ( 

0.9000E+02) 

YES 

PLY 

GROUP 

* 

3  ( 

0 . OOOOE+OO) 

YES 

PLY 

GROUP 

1 

4  ( 

0 . 9000E+02) 

YES 

PLY 

GROUP 

# 

5  ( 

O.OOOOE»OOl 

YES 

DELAMINATION 


INTERFACE  f 
INTERFACE  t 
INTERFACE  « 
INTERFACE  # 


1  (BETWEEN 

2  (BETWEEN 

3  (BETWEEN 

4  (BETWEEN 


PLY  GROUP  » 
PLY  GROUP  » 
PLY  GROUP  I 
PLY  GROUP  » 


1  AW}  I  2) 

2  NO  *3) 

3  NO  14) 

4  NO  t  5) 


NO 

YES 

YES 

YES 


DELAMINATICN  AREA  AW3  SHAPE  AT  THE  INTERFACE  *  2 
FULL  SIZE  <  XI  *  X2  *  0.7U2E-01  •  0.5334E-01) 

GRID  SIZE  -  0.191CE-02 


•  •  4  #  • 

•  •  •  3  •  •  • 


XI 
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X2 


DELAMnftTIOJ  AREA  AW  SHARE  AT  THE  INTERFACE  *  3 
FULL  SIZE  <  XI  •  X2  -  0.7112E-01  *  0.5334E-01) 

GRI1-  SIZE  -  0.1976E-02 


XI 


OELAKINATICN  AREA  AK)  SHAPE  AT  THE  INTERFACE  *  4 
FULL  SIZE  (  XI  ■  X2  -  0.71I2E-01  •  Q.5334E-01) 

GRID  SIZE  -  0.197CE-02 


0  •  •  •  • 


XI 
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